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Abstract 


This  study  examines  a  scenario  for  bolstering  the 
operational  control  exercised  over  the  U.S.  satellite  fleet. 
An  Extended  Space  Transportation  System  (STS)  composed  of  a 
Shuttle,  Space  Station,  and  Orbital  Transfer  Vehicle  (OTV) 
is  analysed  using  a  nonlinear  optimization  technique.  The 
OTV  deploys  a  postulated  fleet  of  military  satellites  across 
the  entire  gamut  of  inclinations  and  altitudes.  The  Shuttle 
payload  mass  and  the  total  Station,  OTV  and  Shuttle 
propellant  mass  required  in  orbit  are  calculated  as  a 
function  of  inclination  and  altitude,  and  used  to  minimize 
annual  Shuttle  launches.  Launch  rates  for  one,  two  and 
three  Station  scenarios  are  compared  with  the  corresponding 
rates  for  current  STS  operations.  'The  use  of  both 
chemically  and  ion  propelled  OTVs  are  evaluated.  Applying  a 
vector  optimization  process  to  the  latter  simultaneously 
minimizes  both  the  average  OTV  mission  duration  and  annual 
Shuttle  launches.  The  resulting  efficient  operating 
frontier  specifies  a  series  of  optimal  inclinations, 
altitudes  and  OTV  sizes  at  which  the  system  should  be 
operated.  Total  Shuttle  launch  rates  for  the  Extended  STS 
are  significantly  less  than  for  direct  orbital  insertion  of 
satellites  with  the  Shuttle.  Equally  important,  the  ion 
propelled  OTV  satellite  deployment  times  are  probably  fast 
enough  to  satisfy  the  military  requirement  of  rapid 
deployment. 


OPTIMIZING  THE  SPACE  TRANSPORTATION  SYSTEM 


1 .0  Introduction 

In  the  past  two  decades  the  United  States  space  program 
grew  from  a  fledgling  experiment  to  a  fleet  of  operational 
satellites  crucial  to  the  civil  and  defense  needs  of  the 
nation.  The  growing  ability  of  military  space  systems  to 
greatly  enhance  the  effectiveness  and  reduce  the  cost  of 
land,  sea  and  air  operations  places  them  in  the  vanguard  of 
United  States  operational  forces.  Space  systems  will  be  the 
first  to  identify  an  enemy  attack  should  it  occur,  and  they 
will  be  used  to  locate,  deploy,  command  and  control  counter¬ 
attacking  forces.  The  recent  introduction  of  the  Space 
Shuttle  and  the  establishment  of  a  USAF  Space  Command  are 
important  interim  steps  toward  consolidating  control  over 
military  satellite  systems. 

An  operational  Space  Transportation  System  (STS)  will 
greatly  enhance  the  effectiveness  of  military  satellites. 
Some  potential  benefits  include:  proliferating  space 
systems  thereby  making  them  more  difficult  to  negate,  adding 
redundant  subsystems  to  satellites  to  make  them  more 
reliable,  adding  shielding  and  other  materials  to  make 
satellites  less  vulnerable  to  radiation  and  other  forms  of 
attack,  and  adding  propellants  to  make  satellites  more 
maneuverable  and  to  increase  their  mission  duration  (Ref.  1) 


1 


A  top-level  defense  study  chaired  by  former  Secretary 
of  the  Air  Force,  Dr.  Hans  Mark,  examined  the  utility  of 
manned  operations  in  space.  The  study  envisioned  the 
Shuttle  being  used  to  assemble  large  structures  in  space,  to 
test  military  systems,  to  repair  valuable  spacecraft,  to  act 
as  a  command  post  during  contingencies,  and  to  be  used  for  a 
variety  of  evolving  missions.  With  the  easy  access  of 
military  personnel  to  orbit,  the  Air  Force  can  develop  new, 
less  complex,  and  less  expensive  satellite  systems  specifi¬ 
cally  designed  for  orbital  repair,  refurbishment,  service 
and  retrieval  (Ref  1). 

Orbital  service  and  test  operations  may  eventually 
reduce  factory  and  launch  base  test  requirements,  and  give 
rise  to  a  sequel  to  the  Factory-to-Pad  test  concept  employed 
with  military  satellites;  namely,  Factory-to-Orbit  testing 
(Ref  2,  p.  10).  A  NASA  study  published  in  1974  examined 
the  performance  of  57  different  spacecraft  during  the  first 
month  in  orbit.  A  total  of  154  malfunctions  were  noted  with 
over  5^i;  occurring  on  the  first  operational  day.  During  the 
subsequent  29  days,  there  was  a  dramatic  decrease  in  the 
number  of  daily  failures  (Ref  3).  Manned  systems  like  the 
Space  Shuttle,  capable  of  returning  or  repairing  a  satellite 
prior  to  final  orbital  insertion,  should  help  to  relieve  the 
serious  problem  of  spacecraft  early  on-orbit  failures. 

Although  the  Space  Shuttle  has  great  potential,  it  also 
has  some  inherent  limitations.  Unable  to  reach  higher 
altitudes,  its  use  for  repair,  refurbishment,  service  and 


retrieval  operations  is  limited  to  Low  Earth  Orbit  (LEO) 
satellites.  The  limited  flight  duration  and  cumbersome 
working  environment  further  qualify  the  Shuttle's  utility  as 
an  orbital  work  platform.  The  Shuttle's  chief  virtue  is 
exactly  what  it  was  designed  for — repetitively  hauling  large 
massive  cargoes  into  LEO.  Unfortunately,  as  currently 
employed  the  Shuttle  cannot  take  full  advantage  of  its 
cargo-hauling  capabilities.  Indeed,  the  Shuttle  is  seldom 
loaded  to  anywhere  near  capacity.  A  logical  extension  over¬ 
coming  the  Shuttle's  drawbacks  is  the  construction  of  a 
manned  Space  Operations  Center  (SOC)  and  Orbital  Transfer 
Vehicle  (OTV)  to  augment  the  STS.  The  Shuttle  launched  to  a 
SOC  could  be  r  sistently  loaded  with  its  maximum  payload  of 
cargo  and  OTV  fuel. 

As  currently  envisioned  the  Extended  STS  composed  of 
the  SOC,  OTV  and  Shuttle  would  be  the  nexus  binding  together 
our  national  space  assets.  From  the  military  perspective 
the  SOC  would  serve  as  a  central  depot  from  which  the  satel¬ 
lite  fleet  could  be  managed.  Some  missions  accomplished  by 
the  Extended  STS  might  include: 

1.  To  provide  a  depot  for  satellite  checkout 
prior  to  final  orbital  deployment. 

2.  To  provide  a  depot  for  satellite  deployment, 
retrieval,  test,  service,  refurbishment  and 
repair . 

3.  To  provide  a  base  for  a  manned/unmanned 
reconnaissance  fleet. 


4.  To  provide  a  base  allowing  for  rapid  replen¬ 
ishment  of  pretested  operational  spacecraft. 

5.  To  provide  a  base  for  an  Anti-Ballistic  Mis¬ 
sile  (ABM)/Anti-Satellite  (ASAT)  system. 

These  cardinal  missions  consolidate  and  bolster  the 
operational  control  exercised  over  allied  and  enemy 
satellites.  Although  not  exhaustive,  these  missions  lend 
themselves  to  quantitative  analysis. 

2.0  Problem  Statement 

The  Extended  STS  envisioned  presents  a  transportation 
and  logistics  problem  similar  in  principle  to  the 
determination  of  an  optimal  resupply  scenario  for  remote 
forces.  The  Shuttle  would  deploy  satellites,  personnel, 
fuel  and  other  consumables  directly  to  the  SOC.  An  OTV 
would  then  be  used  to  deploy  satellites  to  their  final 
orbits.  The  Satellite  Traffic  Model  which  the  OTV  services 
is  a  hypothetical  model  which  crudely  reflects  possible 
military  traffic  requirements  between  1985  and  2000.  Data 
within  the  model  includes  inclinations,  altitudes, 
eccentricities,  masses  and  launch  frequencies  for  a  fleet  of 
postulated  satellites  (see  Table  4.2).  More  realistic 
military  models,  compiled  by  Rockwell  International 
(Ref  4),  are  available  for  authorized  users.  Since  WWII 
such  transportation  problems  have  been  routinely  solved 
using  linear  programming  techniques.  The  Extended  STS 


differs  from  the  traditional  transportation  problem  in  that 
the  orbital  equations  used  are  nonlinear. 

By  utilizing  nonlinear  optimization  techniques  it  is 
possible  to  determine  an  optimal  orbit  and  resupply  scenario 
for  the  SOC.  The  performance  index  used  minimizes  the  mass 
of  fuel  consumed  by  OTV  deployments,  SOC  orbit  maintenance 
and  the  Shuttle  Orbital  Maneuvering  System  (OMS).  This 
total  fuel  mass  is  calculated  as  a  nonlinear  function  of 
inclination  and  altitude.  Alternatively,  the  number  of 
Shuttle  launches  can  be  minimized  by  dividing  the  total  fuel 
mass  by  the  payload  mass  of  the  Shuttle  similarly  calculated 
as  a  function  of  inclination  and  altitude.  The  SOC  altitude 
and  inclination  are  variables  which  the  optimization  proce¬ 
dure  specifies.  Constraints  on  the  optimization  process 
restrict  the  altitude  and  inclination  to  realizable  values 
resulting  in  a  positive  payload  mass  and  an  OMS  fuel  tank 
capacity  less  than  10830  Kg.  Both  models  are  based  on 
accurate  assumptions,  and  both  can  be  further  refined. 

The  launch  and  fuel  minimization  models  directly  assess 
the  launch  and  fuel  costs  associated  with  each  component  of 
the  Extended  STS.  Of  special  significance  is  their  ability 
to  identify  the  costs  associated  with  placing  one  SOC  at  a 
near  equatorial  inclination  versus  two  SOCs,  one  each  at 
polar  and  equatorial  inclinations.  The  sensitivity  or 
robustness  of  the  model  is  analysed  by  varying  parameters 
like  the  SOC  ballistic  coefficient,  OTV  mass  and  the  engine 
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specific  impulses.  The  resulting  impact  on  launch  and  fuel 
costs  gives  the  engineer  cost  saving  insight  into  efficient 
designs  for  an  Extended  STS.  Minor  changes  to  the  models 
allow  for  an  accurate  assessment  of  the  impact  of  increasing 
the  specific  impulses  of  the  SOC  and  OTV  engines  into  the 
range  of  the  ion  motor.  With  the  completion  of  the  sensiti¬ 
vity  analysis  several  operational  scenarios  are  evaluated. 

One  of  the  operational  scenarios  utilizes  a  vector 
optimization  technique  to  minimize  two  performance  indices: 
the  number  of  Shuttle  launches  and  the  average  flight  time 
for  an  ion  propelled  OTV.  The  trade  off  between  the  number 
of  launches  versus  the  flight  times  provides  an  indication 
of  whether  current  ion  engine  technology  is  sufficient  to 
satisfy  the  military  requirement  of  rapid  deployment.  In 
addition,  a  simplified  ABM/ASAT  system  is  postulated  where 
the  flight  time  to  a  coplanar  target  is  minimized  subject  to 
a  given  characteristic  velocity  (  aV).  The  latter  problem 
is  treated  as  an  extension  of  the  basic  thesis,  and  not 
evaluated  in  detail. 

The  remainder  of  this  thesis  includes  five  sections 
explaining  the  solution  methodology,  results,  summary, 
conclusion  and  recommendations.  The  former  derives  the  six 
mathematical  models  used  in  the  optimization  program. 
Results  of  computer  runs  done  for  the  different  models  are 
then  tabulated  and  evaluated.  A  summary  of  the  results 
compares  the  calculated  Shuttle  launch  rates  to  the 
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corresponding  rates  without  an  Extended  STS.  The  concluding 
sections  review  the  key  findings,  recommend  additional 
research  and  recommend  a  tentative  national  space  policy. 
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3.0  Solution  Methodology 
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Minimization  of  the  annual  Shuttle  launches  or  fuel 
consumption  requires  that  the  total  mass  transported  to 
orbit  be  known.  The  total  mass  of  the  SOC,  satellites  and 
consumables  is  easily  calculable  from  the  Satellite  Traffic 
Model  and  a  few  initial  definitions.  Similarly,  the  fuel 
mass  consumed  by  the  Shuttle  OMS,  SOC  orbit  maintenance  and 
OTV  engines  is  calculable  with  only  a  little  extra  effort. 
Combined,  the  mass  of  the  SOC,  satellites,  consumables  and 
propellant  represents  the  total  mass  required  to  deploy, 
maintain  and  operate  the  Extended  STS  in  orbit. 

The  mass  of  fuel  used  by  the  OMS  and  OTV  engines  can  be 
calculated  with  the  Rocket  Equation: 

Mi 

AV  =  Ig  X  ln(-i-)  (Eq.  3.1) 

”0 

where  Mj^  and  Mq  are  the  masses  before  and  after  the  aV 
maneuver.  Specific  impulse  and  gravitational  acceleration 
are  indicated  by  I  and  g  respectively.  Subtracting  the 

final  from  the  initial  mass  and  substituting  Mq  and  from 
Eq.  3.1  gives  the  fuel  mass  (MF)  consumed  by  the  maneuver: 

MF  =  Mi  -  Mo 

AV/Ig 

MF  =  MqCs  -  1 )  (Eqs.  3.2) 

-AV/Ig 

MF=  Mi  (1  -  e  ) 

Consequently,  the  mass  of  propellant  consumed  is  a  function 
of  I,  g,  aV  and  either  the  final  or  initial  mass  of  the 
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rocket.  Alternatively,  the  annual  SOC  orbit  maintenance 
fuel  mass  is  calculated  from  the  definition  of  specific 
impulse : 


Thrust 

(dMF/dt)  X  g 


dMF/dt 


M  X  A 

Ig 


(Eqs.  3.3) 


where  M  and  A  are  the  average  mass  and  acceleration 
respectively.  Although  the  units  for  the  above  equations 
are  arbitrary,  all  of  the  calculations  done  in  this  study 
are  in  MKS  units. 

The  major  simplifying  assumptions  within  the  derived 
mass  models  include: 

1.  The  SOC  orbit  is  circular  between  120  and 
800  km. 

2.  The  atmosphere  is  modeled  as  a  rotating 
sphere  whose  density  decreases  exponentially 
with  altitude. 

3.  Velocity  changing  maneuvers  are  modeled  as 
impulsive  Hohraann  Transfers  for  the 
chemically  propelled  OTV. 

4.  The  chemical  OTV  accomplishes  a  two  impulse 
aV  maneuver  with  the  initial  transfer  being 

an  altitude  change  to  the  destination  orbit. 

A  combined  plane  and  altitude  change  then 
inserts  the  OTV  into  its  destination  orbit. 

To  return  to  the  SOC  the  same  transfers  are 
accomplished  in  the  reverse  order. 

The  outgoing  and  return  aVs  are  assumed 
equivalent  (actually,  this  study  calculates 
the  outgoing  ^V  which  is  slightly  greater 
than  the  return). 


5. 


6. 


Each  satellite  within  the  Satellite  Traffic 
Model  is  deployed  individually  by  the  OTV  (a 
rigorously  realistic  scenario  would  allow  for 
the  simultaneous  deployment  of  satellites 
whenever  possible). 

7 .  Realistic  values  are  assumed  for  design  para¬ 
meters  such  as  specific  impulses,  the  SOC 
ballistic  coefficient  and  the  OTV  mass. 

8.  The  atmospheric  scale  height  is  assumed  con¬ 
stant  at  30  km  (this  gives  an  average  atmos¬ 
pheric  drag  approximately  equivalent  to  that 
experienced  by  Skylab— see  Appendix  B). 

9.  The  generic  Orbiter  used  is  OV-99 
(Challenger)  with  mass  configuration  and 
flight  profile  assumptions  listed  in  Tables 
3.1  and  3.2  respectively.  Two  assumptions 
inherent  within  the  tables  include  not  using 
any  Reaction  Control  System  (RCS)  fuel  during 
flight  and  the  Orbiters  return  to  Earth 
without  any  cargo  on-board. 

10.  The  gravitational  acceleration  (g)  at  an 
altitude  of  390  km  above  the  equator  is  used 
throughout  this  analysis,  8.7  m/sec2. 

11.  The  Orbiter  ascent  aV,  deorbit  aV  and  Main 
Engine  Cut  Off  (MECO)  mass  are  modeled  by 
reqression  analysis  utilizing  flight  planning 
data . 

12.  A  military  Satellite  Traffic  Model  for 
deploying  new  satellites  is  postulated,  while 
the  satellite  retrieval  mission  for  repair 
and  refurbishment  is  ignored. 


All  of  the  assumptions  used  to  derive  the  mass  models  are 
extremely  good  approximations,  and  usually  reflect  a 
slightly  high  estimate  of  the  mass  needed  in  orbit. 
Although  it  would  be  superfluous  to  the  purposes  of  this 
study,  many  of  the  assumptions  could  be  further  refined  or 
eliminated  with  more  elaborate  models  and  a  more  accurate 
assessment  of  the  current  STS  capabilities. 


3.1  OMS  Propellant  Consumption 

In  deriving  an  expression  for  the  OMS  propellant 
consumption  it  is  necessary  to  calculate  the  inserted  mass 
at  MECO.  The  masses  for  various  inclinations  are  estimated 
with  NASA  flight  planning  data  (Ref  5).  Derived  MECO  masses 
are  tabulated  in  Appendix  A  and  modeled  as  a  function  of 
inclination  degrees)  by  linear  regression 

analysis.  The  maximum  3-sigma  error  associated  with  each 
MECO  mass  calculation  is  estimated  to  148  Kg  and  245  Kg  for 
the  Eastern  and  Western  Test  Ranges  (ETR  &  WTR) 
respectively.  The  corresponding  correlation  coefficients 
are  greater  than  0.99999  for  both  cases.  The  resulting 
equations  measuring  mass  in  kilograms  are: 

”^etr  =  159275.829  +  3.3575  x  l^oc 

2  3 

-  3.70994  X  Igoc  +  0.01339  x  Igoc 

MMwtr  =  164628.209  -  79.17494  x  Igoc  (EQS.  3.4) 

2  3 

-  3.01734  X  Igoc  +  0.01256  x  Igoc 

Met  =  38399.0 

The  variables  MM^^,^  and  MM„tr  designate  the  inserted  MECO 
mass  from  the  ETR  and  WTR,  while  Met  designates  the  mass  of 
the  Space  Shuttle  External  Tank  (ET),  residual  fuel  and 
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propellant  required  for  ET  Separation  (ETS).  Table  3.1 
delineates  the  mass  configuration  of  the  generic  Orbiter 
used  to  resupply  the  SOC.  The  STS  configuration  and  assump¬ 
tions  made  in  calculating  the  MECO  masses  are  shown  in  Table 
3.2.  Care  was  taken  to  insure  that  all  assumptions  were 
conservative  estimates  for  the  postulated  STS  configuration. 
Consequently,  Eqs.  3.4  represent  conservative  estimates  of 
the  Shuttle  fleet's  capability  to  nominal  MECO  when 
outfitted  with  HPM  WTR  Motors  and  32945  Kg  (inert  mass) 
Solid  Rocket  Motors  (SRM).  Challenger  was  chosen  as  the 
generic  Orbiter  because  of  its  median  mass  with  respect  to 
the  other  vehicles. 

Two  other  quantities  modeled  by  linear  regression 
analysis  are  the  characteristic  velocities  required  for 
Orbiter  ascent  from  MECO  (aVA)  and  Deorbit  (aVD).  The  data 
for  the  regression  analysis  was  generated  by  the  ascent  and 
deorbit  aV  equations  in  Appendix  A.  The  Appendix  A 
equations  are  used  by  NASA  JSC  personnel  for  mission 
planning  (Ref  6)  and  include  a  velocity  reserve  for 
contingencies . 
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Table  3-1  Mass  Configuration  of  Generic  Orbiter 
(based  on  STS-7  mission) 


ITEM 

Capability  to  Nominal  MECO 

QMS  Propellant  (Ascent) 

QMS  Propellant  (Deorbit) 

Orbiter  OV-99  Inert 
Space  Shuttle  Main  Engines  (SSME) 
Non  Propulsive  Consumables 
STS  mass  charges  to  operator 
Personnel  (2  Men/6  Days) 

RCS  Propellant 

STS  Operations  Reserve 

SUBTOTAL  (reference  mass-M^,) 

RCS  Propellant  for  ET  Separation 
MPS  Unusable  Fluids 
MPS  Flight  Performance  Reserve 
External  Tank  (Block  II) 

SUBTOTAL  (Mass  at  ETS-Met) 


MASS  (Kg) 
+  MM 

MFoa 
-  Ml^od 

68346 

9461 

2166 

2588 

813 

3306 

1361 

88041 

100 

3991 

2518 

31790 

38399 


Payload  Mass 


Table  3-2  Flight  Profile  Assumptions  for 


MECO  Mass  Calculation 


ETR  Launch  WTR  Launch 


-Std  MECO  Conditions 
-Altitude  =  57  NM 
-Fit  Path  Angle=0.65° 
-Iner.  Vel  =  25680  fps 
-102?  SSME  Power  Level 
-2.75  Sec  delay  till 
SRM  Ignition 
-3o  Flight 
Performance  Reserve 
-15  July  Launch 
-680  PSF  Pressure 
-HPM  WTR  Motor 
-Light  Weight  SRM 
-Light  Weight  ET 
-OV-99  (see  Table  3-1) 


-Std  MECO  Conditions 
-Same 
-Same 

-Iner.  Vel  =  25374  fps 
-Same 
-Same 

-Same 

-15  March  Launch 

-650  PSF  Pressure 

-Same 

-Same 

-Same 

-Same 


AVAetr  =  (-2.023703587  x  10-5)  ^  R^oc 

+  0.8474467398  x  Rgoc  -  4623.658232 


AVAwtr  =  ^VAetr  *  85.65  m/s 


(Eqs.  3.5) 


aVD  =  (1  .974610171  X  10-^)  X  R. 


-  2.50173195  X  +  7994.436744 

The  correlation  coefficients  calculated  in  Appendix  A  are 
0.99999  for  ascent  and  0.99696  for  deorbit.  It  is  worth 
noting  that  the  ascent  aVs  could  be  calculated  exactly  by 
using  the  ETR  and  WTR  standard  MECO  conditions  shown  in 
Appendix  A.  Similarly,  a  good  estimate  for  the  deorbit  aV 
could  be  calculated  by  targeting  the  maneuver's  perigee  for 
a  constant  altitude  within  the  atmosphere. 

Utilizing  the  curve-fitted  data  and  Eqs.  3.2  the  Shuttle 
OMS  fuel  consumption  for  ETR  launches  is: 


,,,  (-^VAetr/Ig)^ 

MFoa  =  (MMetr  -  ) 


( AVD/Ig) 

MFod  =  (MMetr  “  Met  “  MFoa  ’  Mpi)(1-e  ) 


(Eqs.  3.6) 


where  Mpj^^  MFoa  MFod  the  Orbiter  payload  mass,  the 
propellant  mass  required  for  Orbiter  ascent  and  the 
propellant  mass  required  for  deorbit.  The  payload  mass  of 
the  Orbiter  is  calculable  if  the  Orbiter  reference  mass 
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calculated  in  Table  3.1,  is  known  (Mj.  is  approximately 
equal  to  the  Orbiter  reentry  mass): 

”pl  =  MMetr  “  ^et  "  MFoa  -  MFod  -  (Eq.  3.7) 

Solving  Eqs.  3.4,  3.6  and  3.7  simultaneously  for  the  payload 
mass  yields: 


(-^^^etr/Ig) 

Mpi  =  —  Mgt)® 


( AVD/Ig),_ 

Mr  X  e  (Eq. 


3.8) 


The  WTR  equations  corresponding  to  Eqs.  3.6,  3.7  and  3.8  are 

easily  derived  by  substituting  MM„tr  and  ^VAv^tr  MMetr 
and  AVAg^j..  All  three  equations  are  predicated  on  not 
using  any  RCS  propellant  during  flight.  This  assumption  is 
reasonable  with  careful  mission  planning,  and  even  when 
ignored  results  only  in  a  conservatively  high  estimate  of 
the  OMS  propellant  mass  consumed  in  orbit. 
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3.2  SOC  Orbit  Maintenance  Propellant  Consumption 


The  mass  of  SOC  orbit  maintenance  propellant  can  be 
estimated  by  first  calculating  the  acceleration  due  to  air 
drag.  From  the  definition  of  impulse: 

D  X  dt  =  MgQ^  X  dVsoc 
or 

■^soc  ~  ^/^soc  =  P  X  F  X  S  X  C(j  X  VpQ/2MgQ2  (Eq,  3.9) 

where  Agoc  acceleration  of  the  SOC  due  to  air  drag 

(D).  The  individual  drag  terms  are  defined  in  King-Hele 
(Ref  7f  p70p.  20-26)  as: 

~  2.2  =  drag  coefficient 
S  =  average  cross  sectional  area 

^soc  =  mass  of  SOC 

B  =  ^CciS)/Msoc  =  ballistic  coefficient  (m^/kg) 

F  =  ^(Rpoj^soc^  ~  1  “Rpo^^e  ®  ^  ^  SOC ) ) /^  po  ^  ^ 

^  =  2  -  fCRpo  ,  180°  -  Isoc)  90°<Isoc<180 

F  =  0.00175  X  Igoc  0.84004  250<Isoc<‘' 

^po  =  Rsoc  =  perigee/SOC  radius 

(Eqs.  3.10) 

^soc  =  80C  inclination 

^po  =  ^O/Rgoc)^'^^  =  perigee  velocity 

We  =  7.292115856  x  10~5  rad/sec  =  earth's  angular 

velocity 

U  =  398601 .2  km3/sec2  =  gravitational  parameter 
R  ”  Rq  X  “  Rsoc)/^) 

Rq  =  120  km  =  reference  altitude 


120  km  =  reference  altitude 


Pq  =  24.9  Kg/km^  =  atmospheric  density  at 
H  =  30.0  km  =  scale  height 


The  equation  specified  by  F  represents  the  effect  of 
atmospheric  rotation  on  drag,  and  varies  between  about  0.9 
and  1.1  over  the  inclination  range.  This  rotating  atmos¬ 
phere  factor  is  shown  for  inclinations  below  and  above  90°, 
and  is  also  modeled  by  regression  analysis  for  the 
inclination  range  of  interest.  The  extreme  sensitivity  of 
the  atmospheric  model  resulted  in  unrealistically  high  esti¬ 
mates  of  fuel  consumption  when  a  linear  fit  to  the  Scale 
Height  was  attempted.  Consequently,  the  constant  value  of 
Scale  Height  was  chosen  which  resulted  in  an  atmospheric 
drag  approximately  equal  to  that  experienced  by  the  Skylab 
mission  (see  Appendix  B).  Substituting  the  terms  of 
Eqs.  3.10  into  Eq.  3-9  gives: 

Asoc  =  0  X  Po  x  B  X  (0.00175  x  Igoc  +  0.84004) 

(Pq  “  Rsoc^/H) 


'O  ~  ^SOC‘ 

2Rsoc 


(Eq.  3.11) 


and  utilizing  Eq.  3-3: 

^^SOC  =(^300  X  A3QC  )  31557600/Ig 


(Eq.  3.12) 


where  MF30C  average  annual  mass  of  propellant 

required  to  keep  the  SOC  in  orbit. 
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3.3  OTV  Propellant  Consumption 


The  OTV  propellant  mass  required  to  deploy  or  retrieve 
a  satellite  is  derived  by  first  calculating  the  required 
aV.  Figure  3.1  depicts  the  two-impulse  Hohraann  Transfer 
used  by  the  maneuver.  The  first  impulse  accomplishes  an 
altitude  change  to  the  mission  orbit's  apogee,  then  a 
combined  inclination  and  altitude  change  places  the  OTV  in 
its  destination  orbit. 


Destination 

Satellite 

Orbit 


Fig  3.1  Hohmann  Transfer 


The  equations  for  the  first  impulse  are: 


Vsoc  “  (U/Rsoc^ 

^to  “  (2U(1/Rsoc  "^/^^asat  ^soc^^^^^ 

(Eqs.  3.13) 

^asat  =  ^sat^^  ■*■  ®sat^ 

=  Vto  -  Vsoc 

and  for  the  second  impulse: 

^tf  "  (2U(1/R333t  “  ^/^^asat  ■*■  l^soc^^^^^ 

Vsat  =  (2U(1/R3sat  -  1/2Asat))l/2  (Eqs.  3.14) 

AV2  =  (Vtf2  +  ^sat^  -  Vtf  X  Vsat 

X  cosdsoc  -  Isat^^ 

Recognizing  that  the  aV  required  to  travel  from  the  SOC  to 
the  satellite  mission  orbit  is  approximately  equal  to  the 
return  aV  (  ”  ^Vj.goc^  »  ^nd  combining  Eqs.  3.13 

and  3.14: 

^^soc-i  =  4Vi  +  AV2  (Eq  3.15) 

Solving  Eqs.  3.2  and  3.15  simultaneously  gives  the  total  mass 
of  fuel  required  to  service  a  satellite  orbit  with  inclina¬ 
tion  dsat)  »  semi-major  axis  (Asat^  »  eccentricity  (esat^ 
and  satellite  mass  ; 


Mfrtn  =  Motv^® 


^^^soc-i/Ig) 


-  1) 


(Eqs.  3.16) 


( AVsoc-i'^^®^ 

MFdpl  =  ^^sat  ■*■  ^otv 

where  is  the  OTV  mass,  MF^pi  is  the  propellant  mass 

required  to  deploy  a  satellite  and  MF^tn  propellant 

mass  required  to  return  to  the  SOC.  Summing  the  two 
equations : 


(Msat  ^otv  ^otv  *  ® 


^  ^^soc-i/^g^ 


X  (e  -  1 ) 


(Eq.  3.17) 


where  MF^sn  mass  of  propellant  required  per  OTV 

mission  depolying  a  payload  and  returning  to  the  SOC. 


3.4  Calculation  of  Launch  Windows 

The  Hohmann  Transfer  specifies  all  of  the  orbital 
elements  except  the  argument  of  perigee  (id)  and  the 
longitude  of  right  ascension  (n).  Both  parameters  can  be 
set  by  choosing  an  appropriate  launch  window  for  the  OTV. 
Alternatively,  the  frequency  at  which  launch  windows  occur 
can  be  calculated  if  the  SOC  and  satellite  orbits  are  known. 
The  window  occurs  once  during  each  synodic  period  of  the 
satellite  and  SOC  argument  of  perigee.  Similarly,  the 
window  occurs  twice  during  each  synodic  period  of  the 


satellite  and  SOC  longitude  of  right  ascension.  The  u 
window  occurs  much  more  frequently  than  the  corresponding 
n  window  which  is  calculated  from  the  definition  of  synodic 
period  (Tg); 


^•^sat/*^^  “  dA5QQ/dt 

-3  X  J2  X  cosdsat)  x  (U/Agat^  3-18) 

dflsat^^^  =  - 

2  . 2 

2(1  -  esat^ 

7  , 

ddgoc^^^  s  «3  X  J2  X  cos(IgQQ)  X  (U/RgQQ/ 

J2  =  0.001082642  (Ref  8,  p.  422) 

where  J2  is  a  constant,  and  dnsat/*^^  ^nd  dfigoc^*^^ 
time  rates  of  change  of  the  longitude  of  right  ascension  for 
the  satellite  and  SOC  respectively.  Similarly,  Eqs.  3.18 
could  be  used  to  calculate  the  synodic  period  between  the 
satellite  and  SOC  argument  of  perigee  by  replacing  dfigat'^*^^ 
and  dfigQ^/dt  with  the  satellite  and  SOC  orbital  rates. 
However,  since  the  propellant  cost  associated  with  orbital 
transfers  to  a  less  than  optimal  u*  is  negligible  in  compari¬ 
son  with  similar  n  transfers  they  are  ignored  in  this 
analysis . 

Besides  being  useful  for  deriving  launch  windows,  the 
time  variation  of  right  ascension  might  also  be  useful  for 
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deploying  spare  orbital  satellites  or  even  an  ABM/ASAT 
system.  Objects  deployed  at  the  SOC  inclination  but  with  a 
different  altitude  have  different  nodal  regression  rates 
and,  consequently,  different  ground  tracks  over  the  surface 
of  the  Earth.  Spare  satellites  deployed  in  this  manner  are 
geographically  dispersed  over  the  surface  of  the  Earth,  and 
yet  easily  serviced  at  periodic  intervals  with  only  a  small 
coplanar  Hohmann  Transfer  from  the  SOC.  An  ASAT  system 
deployed  similarly  has  two  coplanar  intercept  opportunities 
during  each  orbit  of  a  target  satellite.  The  corresponding 
ABM  scenario  is  more  complex  in  that  there  is  only  one 
launch  opportunity,  and  consequently  many  more  interceptors 
are  required.  The  geographic  dispersal  of  all  the  systems 
makes  them  difficult  targets  which  must  be  destroyed  one  at 
a  time. 


3.5  Single  Objective  Optimization  Problems 

Several  nonlinear  optimization  problems  are  posed  and 
solved  utilizing  a  numerical  routine,  the  Sequential 
Unconstrained  Minimization  Technique  (SUMT)  computer  program 
(Ref  9).  This  program  uses  a  penalty  function  approach  to 
find  the  minimum  of  a  single  multivariable  nonlinear 
function  (called  an  objective  function  or  performance  index) 
subject  to  inequality  and  equality  constraints: 


MIN:  F(Xi,  Xp,  ...,  Xn) 
SUBJECT  TO  CONSTRAINTS: 


Gk(Xi,  Xp,  ...»  Xn)  >  0 


Hk(Xi,  X2 . Xn)  =0 

k=M+T,  M+2,  M+MZ 

The  procedure,  developed  by  Fiacco  and  McCormick  (Ref  10), 
uses  the  problem  constraints  and  the  original  objective 
function  to  form  an  unconstrained  objective  function  which 
is  minimized  by  any  appropriate  unconstrained  multivariable 
technique.  A  requirement  on  the  F,  G  and  H  functions  is 
that  they  be  continuous  and  twice  differentiable.  Fig  3.2 
is  a  brief  flowchart  of  the  SUMT  algorithm.  The  performance 
indices  individually  minimized  by  SUMT  include  the  total 
f;.>el  consumption  and  number  of  Shuttle  launches  for  an 
Extended  STS. 

cr 

3.5.1  Model  I — ETR/WTR  Fuel  Minimization  Model 

With  all  of  the  propellant  mass  terms  defined,  the 
minimization  problem  is: 


MIN:  MF^q^  =  ^orb  ^  (MFq3  +  MFqj)  +  MF50C 


iJi  ^^sat  ^ 


SUBJECT  TO  CONSTRAINTS: 


28.5° 

<  Isoc 

<  570 

560  < 

^soc  ^ 

104° 

6500 

km  <  Rsoc  ^  7200 

MF 

“^oa 

+  MFod  < 

10830  Kg 

Mpl  > 

0  Kg 

^SOC  t 

Rsoc 

(Eqs.  3.19) 

ETR  Inclination  constraint 
WTR  Inclination  constraint 
Altitude  constraints 

OMS  fuel  constraint 
Orbiter  payload  constraint 

Problem  variables 
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Fiacco  and  McCormick  (SUMT  Algorithm) 
Logic  Diagram 


where  all  terms  have  been  previously  defined  except  the 
annual  satellite  launch  frequency  (Fgat)*  the  total  fuel 
consumption  (MFtot^’  number  of  satellite  missions  (n) 
and  the  number  of  Shuttle  launches  The  latter  can 
either  be  estimated  or  calculated  more  rigorously  by: 


F orb  ~ 


i  =  1  t^sat  ^  (Mgat+Mfmsn)  3  ^^soc  +  Me 


(Eq.  3.20) 


where  Me  includes  the  mass  of  the  SOC,  OTVs  and  consumables 
(food,  oxygen,  men,  equipment,  etc.)  required  to  build  and 
operate  the  station.  Some  estimates  of  the  required  masses 
are:  (Ref  11,  p.  189) : 

Msoc  ^  Motv  100000  Kg 


'food  “  210  Kg/yr/man 


Mwater  ~  KG/yr/man 


^air  ~  330  KG/yr/man 


^misc 


5920  Kg/yr 


Averaging  the  above  for  a  six  ’-an  station  over  a  16  year 
period  gives  an  of  20000  Kg/yr. 


3.5.2  Model  II--ETR/WTR  Launch  Minimization  Model 

As  an  alternative  to  the  fuel  minimization  problem  the 
number  of  Shuttle  launches  are  also  minimized. 
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MIN:  Forb 


i^l  '■^sat^Msat  MFmsn)]  +  MFsoc  Me 


SUBJECT  TO  CONSTRAINTS: 


(Eqs.  3.22) 


28.5°  <  <  57° 


ETR  constraint 


56°  <  Isoc  < 


WTR  constraint 


6500  km  <  R  <  7200  km 


Altitude  constraint 


MFoa  +  MFod  <  10830  Kg 


OMS  fuel  constraint 


Mpi  >  0  kg 


Payload  constraint 


Isoc>  ^soc 


Problem  Variables 


3.5.3  Model  III--ETR/WTR  Launch  Minimization  Model  with 
OTV/Shuttle  Rendezvous  in  LEO 
One  possible  method  of  further  reducing  the  number  of 
Shuttle  Launches  is  to  deploy  the  OTV  to  a  low  altitude 
where  it  would  rendezvous  with  the  Shuttle.  After 
rendezvous  the  OTV  would  pick  up  the  Shuttle  payload  and 
sufficient  fuel  for  the  return  journey  to  the  SOC.  Although 
manned  and  unmanned  cargo  modules  would  be  needed,  they 
should  not  be  difficult  to  build  and  could  be  reused. 
Alternatively,  the  necessity  for  manned  modules  could  be 
eliminated  by  occassionally  deploying  the  Shuttle  all  the 
way  to  the  SOC  for  personnel  changes. 

To  implement  the  changes  to  Model  II  a  Hohmann  Transfer 
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is  calculated  from  the  rendezvous  inclination  and  altitude 


to  the  SOC.  The  rendezvous  inclination  is  left  as  a  vari¬ 
able,  while  the  altitude  is  set  at  the  maximum  altitude 
attained  by  the  Shuttle  during  ballistic  flight  after  MECO 
(i.e.  the  apogee  altitude  for  the  standard  MECO  conditions). 
The  Shuttle  rendezvous  radius,  and  the  corresponding  ascent 
and  descent  aVs  are: 

^ren  =  653d. 1  km  (Altitude  =  160  km) 

Iren  =  variable 

^^^etr  =  51.7  m/sec  (Eqs.  3.23) 

AVA^tr  =  137.4  m/sec 
aVD  =  84.1  m/sec 


where  and  I^en  3'"®  rendezvous  altitude  and 
inclination.  The  values  of  aVA  and  aVD  were  calculated 
utilizing  Eqs.  3.5.  Substituting  into  Eqs.  3.13,  3.14  and 

3*15  Rpen  for  Rsoo  ^soc  ^asat  ^sat’  Hohmann 
Transfer  characteristic  velocity  is  calculable. 

Utilizing  Eqs.  3.2,  3-23  and  AVren  mass  required 
for  the  OTV  to  descend  to  (and  ascend  from)  the  rendezvous 
orbit  is  calculated  as  follows; 


,  ^ ^^ren/^S) 

MFdsc  “  ^otv  X  ^®  ”1) 

(Eqs.  3.24) 

^^asc  ~  ^^otv  ^pl ^  X  (1  -  e  ) 


MFdsc  and  MFasc 


where  all  of  the  variables  except  ^ren> 
are  defined  as  in  Eqs.  3.16.  The  last  minor  change  to  the 
model  requires  the  addition  of  MF^gj,  into  the  numerator,  and 
the  subtraction  of  MFgsc  fnom  the  denominator  of  Eq.  3.20 
The  resulting  Orbiter  launch  frequency  is: 


^  ^^sat^^sat+f^^rnsn^  ]  +  MFgoc  +  Mq  +  MFdsc^^f'orb 
i=  1 


(Eq.  3.25) 


Simplified  and  restated  the  minimization  problem 
becomes : 


MIN :  f orb 


n 

1=1  C^^sat^^sat+^^rnsn^  ^  *  ^^soc  ^c 


(Mpi  -  MFdsc  “  MFasc) 


SUBJECT  TO  CONSTRAINTS:  (Eqs.  3.26) 

(same  as  in  Eqs.  3.22  except  Ipg^  is  also  a  variable) 


3.5.4  Model  IV— ETR/WTR  Launch  Minimization  Model  with  Ion 
Propelled  OTV 

Due  to  the  high  fuel  cost  of  deploying  satellites  with 
the  OTV  more  efficient  engine  designs  may  be  desirable.  A 
likely  candidate  is  the  ion  engine  which  can  deploy  payloads 
with  vastly  less  fuel  than  required  by  a  chemical  vehicle. 
The  major  drawback  of  the  ion  system  is  its  very  low  thrust 
which  results  in  lengthy  missions.  Unfortunately,  the 
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Hohmann  Transfer  yields  an  optimistically  low  estimate  of 
the  velocity  change  required  for  constant  low  thrust  orbit 
changes . 

A  better  model  partially  developed  at  the  Air  Force 
Institute  of  Technology  requires  the  deletion  of  assumptions 
3,  4  and  5  of  Section  3.0  and  the  addition  of  two  new 
assumptions : 


1 .  Propulsion  system  thrust  is  low  requiring  an 
infinite  number  of  revolutions  to  reach  the 
final  orbit. 

2.  The  radius  change  is  accomplished  first 
followed  by  the  inclination  change. 

3.  The  ion  OTV  mass  is  approximated  with  Eq.  3.30. 


Applying  these  assumptions  Alfano  and  Wiesel  derived  minimum 
flight  times  for  radius  and  plane  changes  accomplished 
independently  (Ref  12): 


-1/2  -1/2 
So  -  af 


X  t  X  U 


-1/2 


(Eqs.  3.27) 


lo  ’  If  = 


2t  X  t  X  (a/U 


J/2 


The  semimajor  axis  and  inclination  are  designated  by  a  and 
I,  while  T  and  t  represent  the  thrust  to  mass  ratio  and  the 
mission  duration.  The  gravitational  parameter  is  represen¬ 
ted  as  before  by  U.  The  same  investigators  also 
successfully  evaluated  combined  plane  and  radius  changes 
which  were  slightly  more  efficient,  but  their  implementation 
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is  much  more  complex. 


i. 

i 

I 


Applying  Eqs  3.27  to  the  OTV  orbit  change  gives  AVs 
for  the  outgoing  and  return  trips: 

-1/2  -1/2 


aVI  = 

^sat 

“  Rsoc 

TT 

AV2  =■ 

Isat  “  Isoc 

^^tot 

=  AVI  +  aV2 

X 


-1/2 

Asat 


(Eqs.  3.28) 


where  represents  the  total  velocity  change  required  in 
canonical  units.  These  aV  equations  can  be  used  as  in  Eq. 
3.16  to  define  the  fuel  mass  required  to  deploy  and  retrieve 
satellites : 


I 

I 
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MFrtn  =  ^otv  ^ 


(e 


(Eqs.  3.29) 


MFdpl  =  tMotv  ^sat  .  ^  ® 


where  ail  of  the  variables  except  aV^-^j.  are  defined  as  in 
Eq.  3.16.  An  interesting  alternative  to  using  a  constant 
OTV  mass  is  to  define  the  mass  in  terms  of  the  number  of  ion 
thruster  used  (Ref  13,  p.  82): 


^  ^otv  =  N  X  (Mthr  ^sp  *  ^wr  ^ft^  ^str  (Eq.  3.30) 

where  N,  Mthr»  Msp,  Pwr»  Mft  and  Mstr  number  of 

thrusters,  thruster  mass,  specific  mass,  thruster  power, 
I  fuel  tank  mass  and  the  mass  of  the  OTV  structure  and  control 

systems  respectively. 
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Deployment  and  return  times  for  the  OTV  can  be 
calculated  by  integrating  the  mass  flow  rate,  Eq.  3.3: 

dMF/dt  =  (N  X  T)/Ig 
or 

TIM  =  (MF  X  Ig)/N  X  T  (Eq.  3-31) 

The  flight  time,  TIM,  is  thus  easily  calculable  once  the 
required  fuel  mass  and  individual  motor  thrusts,  T,  are 

known.  Applyng  Eqs.  3.31  and  3.29  allows  for  the 

calculation  of  the  OTV  deployment  and  return  times  for  each 
mission  within  the  Satellite  Traffic  Model: 

™rtn  =  MFrtn  x  Ig/(N  x  T) 

-  MFdpi  X  Ig/(N  X  T)  (Eq.  3.32) 

where  TIM^tn  and  TIM^pi  represent  the  return  and  deployment 
times  respectively. 

The  resulting  problem  is  identical  to  the  ETR/WTR 
Launch  Minimization  problem  in  section  3.5.2,  except  that 
the  A  Vs  required  for  OTV  orbital  changes  are  defined 
by  Eqs.  3*28.  In  addition,  the  OTV  flight  times  are 

defined,  and  could  be  minimized  in  place  of  the  number  of 

STS  launches.  Utilizing  Eq.  3.30  and  varying  the  ion  engine 
design  parameters  also  provide  insight  into  the  feasibility 
of  an  ion  propelled  OTV  for  rapid  satellite  deployment. 
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3.5.5  Model  V--ETR/WTR  Launch  Minimization  Model  with 
Chemical  OTV/STS  Rendezvous  in  LEO  and  Ion  OTV  for  Satellite 
Deployment 

The  final  '•operational"  scenario  considered  combines 
the  methods  of  sections  3.5.3  and  3.5.4.  A  chemical  OTV  is 
again  used  to  rendezvous  with  the  Shuttle  and  pick  up  pay- 
load  and  fuel  for  the  return  journey.  Concurrently,  ion- 
propelled  OTVs  deploy  satellites  to  their  operational 
orbits.  The  final  minimization  problem  is  identical  to  that 
in  section  3.5.3,  except  the  ion-propelled  OTV  mass,  delta- 
vees  and  mission  durations  are  calculated  by  Eqs.  3.30,  3.28 
and  3.32  respectively.  Also,  the  OTV  rendezvous  inclination 
is  set  equal  to  the  SOC  inclination  to  reduce  the  number  of 
variables.  The  resulting  two  variable  problem  is 
considerably  easier  for  the  numerical  technique  SUMT  to 
optimize . 
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3.6  Multiple  Objective  Optimization  Problems 


In  addition  to  problems  containing  only  one  objective 
function,  similar  problems  can  be  formulated  with  a  vector 
of  several  objective  fuj.ctions.  A  vector  or  multiple 
objective  optimization  technique  optimizes  the  vector  of 
objective  function  with  a  comuter  program  such  as  PROCES. 
The  procedure  (Ref.  14,15,21  &  22)  sequentially  minimizes 
each  objective  function  subject  to  the  original  problem  con¬ 
straints  and  internally  imposed  equality  constraints  due  to 
the  remaining  objective  functions.  The  sequential  minimiza¬ 
tion  problem  is  accomplished  by  using  the  SUMT  optimization 
program  as  a  subroutine  within  PROCES. 

The  solution  of  the  problem  is  a  Non-Dominated  Solution 
Set  (NDSS)  known  as  an  efficient  frontier.  Specific  values 
for  the  problem  variables  (for  example  SOC  inclination, 
radius,  etc.)  are  associated  with  each  point  of  the  NDSS. 
Figure  3.3  shows  a  dual  objective  optimization  scheme  where 
each  axis  represents  a  performance  index,  J1  or  J2. 


NDSS  or 

efficient  frontier 


Fig.  3.3  Dual  Objective  Optimization  Scheme 
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In  the  two  dimensional  case  the  solution  is  analogous  to  the 
guns  vs.  butter  dichotomy  commonly  cited  in  economic  theory 
However,  PROCES  allows  the  user  to  find  the  actual  values  of 
points  on  the  efficient  frontier. 

Two  problems  with  dual  objective  functions  lend 
themselves  to  vector  optimization  solutions.  The  first  is  a 
variation  of  the  launch  minimization  models  of  Section  3.5. 
However,  in  addition  to  the  number  of  Shuttle  launches  the 
average  ion-propelled  OTV  Time-of-Flight  (TOF)  is  minimized. 
A  second  vector  optimization  application  may  provi  le  a  means 
to  evaluate  an  ABM/ASAT  system.  The  ABM/ASAT  TOF  is 
minimized  subject  to  different  av  and/or  true  anomaly 
constraints . 

3.6.1  Model  VI — Minimized  Shuttle  Launches  vs  Average 
Ion-Propelled  OTV  Time  of  Flight 

Keeping  the  number  of  Shuttle  launches  to  a  minimum  is 
not  the  only  important  consideration  when  designing  an 
Extended  STS.  For  the  case  of  an  ion  propelled  OTV  the 
vehicle's  average  TOF  is  often  equally  important.  A  TOF 
measured  in  years  contravenes  the  occasional  military 
requirement  for  rapid  deployment.  Adoptance  of  an  on-orbit 
spares  philosophy,  whereby  inactive  operational  satellites 
are  deployed  in  advance  of  need,  partially  obviates  the 
military  requirement  for  rapid  deployment.  However,  engin¬ 
eering  common  sense  and  the  satellite  lifetime  dictates  that 
transit  time  be  kept  to  a  minimum. 
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When  the  average  TOF  for  OTV  missions  and  the  number  of 
Shuttle  launches  are  simultaneously  minimized  a  graph 
similar  to  Fig.  3.3  is  generated.  The  graph  depicts  the 
minimum  number  of  Shuttle  launches  as  a  function  of  the 
average  TOF  for  OTV  missions.  The  relationship  between  the 
TOF  and  STS  launches  (called  the  NDSS  or  efficient  frontier) 
provides  the  designer  with  an  important  decision-making 
tool;  namely,  a  trade  off  analysis  between  cost  (STS 
launches)  and  operational  utility  (TOF). 

Application  of  a  vector  optimization  technique  to  ETR 
Model  V  is  a  relatively  straight  forward  procedure.  The 
resulting  dual  objective  minimization  problem  is: 


MIN:  Forb 


n 

i  =  1  ^^sat* ( Ms3t+MF dpl+MF|-tn)  ^ "^MF soc^Mc+NOxMq^v 
”pl  -  MFasc  “  MFdsc 


and 


TIMavg 


n 

i^^TIMdpi  TIMrtn) 
n 


SUBJECT  TO  CONSTRAINTS: 

28.5°  <  igoc  <  570 
6500  km  <  RgQ^,  <  7200  km 
MFoa  +  MFod  <  ^^0830  kg 

Mpi  >  0  kg 
N  >  0 
NO  >  0 

Isoc»  8soc»  ^ 


(Eqs.  3.33) 
ETR  constraint 

Altitude  constraint 
OMS  fuel  constraint 

Payload  constraint 
Number  of  ion  thrusters 
Number  of  OTVs 
Problem  variables 
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Although  most  of  the  expressions  in  the  model  are  standard  a 
few  new  ones  appear  and  several  need  clarification. 
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The  new  expressions  include  TIMg^g  and  NO  which 
respectively  represent  the  average  TOF  for  an  ion  propelled 
OTV  and  the  average  number  of  OTVs  required  to  service  the 
Satellite  Traffic  Model.  The  average  TOF  is  calculated 
above  while  NO  is  given  by: 


NO 


(TIMdpi  TIMrtn)F‘sat 

365.25 


(Eq.  3.34) 


Because  TIM  is  measured  in  days/mission  and  Fgat 
missions/year  the  365.25  days/year  conversion  factor  must  be 
included.  The  quantities  MF^pj^  and  MF^tn  represent  the  fuel 
mass  consumed  by  the  ion-propelled  OTV,  while  MFasc 
^^dsc  represent  the  fuel  mass  consumed  by  the  chemically 
propelled  OTV  when  rendezvousing  with  the  Shuttle. 

The  ion-propelled  OTV  mass  is  again  calculated  by 
Eq.  3.30: 


^otv  =  N  X  (Mthr  ^sp  ^  ^wr  ^ft^  +  f^str 

However,  unlike  the  previous  case  the  number  of  ion 
thrusters,  N,  is  left  as  a  variable  which  the  optimization 
procedure  specifies.  Consequently,  this  problem  has  three 
variables:  SOC  inclination,  SOC  radius  and  the  number  of 

ion  thrusters  on  the  OTV. 
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3.6.2  Model  VII--Minimized  ABM/ASAT  Time  of  Flight 

Optimizing  the  Extended  STS  is  not  the  only  useful 
vector  optimization  application.  The  technique  might  also 
be  applied  to  minimize  the  TOF  to  a  target  satellite  or 
ballistic  missile.  To  simplify  the  problem  only  coplanar 
transfers  are  considered  with  intercept  occurring  at  the 
node  connecting  the  interceptor  and  target  orbits.  It  is 
possible  to  intercept  target  satellites  frequently  by 
deploying  a  string  of  interceptors  at  the  same  inclination 
as  the  SOC.  An  altitude  difference  between  the  SOC  and 
interceptor  orbit  specifies  the  nodal  regression  rate 
between  the  SOC  and  interceptor  longitude  of  ascending 
nodes.  Consequently,  an  OTV  based  at  the  SOC  could  service 
the  interceptor  orbit  at  periodic  intervals  with  only  a 
small  altitude  change.  The  service  interval  for  the  OTV  is 
the  synodic  period  between  the  two  nodal  regression  rates 
and  is  calculated  by  Eqs.  3.18. 

An  ASAT  system  deployed  in  this  manner  has  two  coplanar 
intercept  opportunities  during  each  orbit  of  a  target 
satellite.  The  ABM  scenario  is  much  more  complex  in  that  a 
variety  of  different  interceptor  orbits  are  probably 
required.  The  ABM  interceptors  would  have  to  be  deployed  at 
different  inclinations  in  order  to  intercept  both  ICBMs  and 
SLBMs.  Further,  the  Earth's  rotation  (and  therefore  the 
rotation  of  potential  launch  sites)  requires  that  the 
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interceptors  at  a  given  inclination  be  deployed  at  different 
Longitudes  of  Right  Ascension. 

Although  the  ABM  system  is  considerably  more  complex  than  an 
ASAT  system  it  has  the  added  advantage  of  possibly  being 
useful  against  airborne,  ground  and  naval  targets. 

Minimizing  the  TOF  is  fairly  simple  when  applied  to  a 
specific  case  (Ref.  8,  p.  I8l).  However,  the  general  case 
is  considerably  more  complex.  Applying  a  vector  optimiza¬ 
tion  technique  to  the  problem  might  allow  for  the 
simultaneous  minimization  of  interceptor  TOF,  aV  and  true 
anomaly  at  epoch.  In  essence,  the  solution  would  be  a  NDSS 
with  TOF  given  as  a  function  of  aV  and/or  true  anomaly. 
The  transfer  orbit  parameters  would  be  problem  variables 
specified  by  the  optimization  procedure.  These  values 
indirectly  specify  the  mass  and  quantity  of  ABM/ASAT 
interceptors  required  to  achieve  a  given  TOF.  The  problem 
is  not  simple  but  a  successful  solution  would  allow 
decision-makers  to  quantify  the  cost  and  effectiveness  of 
spaceborne  interceptors  in  future  conflicts.  Such  a 
solution  may  be  well  worth  the  invested  effort. 


4.0  Results 


Numerous  computer  runs  were  completed  on  Models  I 
through  VI  of  Section  3.0.  Model  VII  was  analysed  qualita¬ 
tively  but  not  implemented  on  the  computer.  Appendix  D 
contains  a  listing  of  the  operational  computer  runs 
accomplished  and  selected  outputs.  In  addition  to  the 
operational  runs  of  Appendix  D  several  test  runs  were 
completed  on  the  basic  fuel  and  launch  minimization  models. 
Table  4.1  lists  the  standard  problem  parameters  which  were 
used  in  the  test  and  operational  computer  runs.  Variation 
of  these  parameters  in  the  operational  runs  provided  insight 
into  the  best  means  of  designing  an  Extended  STS.  Table  4.2 
shows  the  Satellite  Traffic  Models  used  for  the  test  and 
operational  runs. 

4.1  Model  Validation 

Validation  of  the  basic  model  was  accomplished  by 
applying  the  ETR  model  to  three  test  missions.  The  missions 
shown  in  Table  4.2  include  orbital  parameters,  masses  and 
launch  frequencies  for  three  satellites.  The  SUMT 
minimization  program  was  run  with  one,  two  and  then  all 
three  test  missions  used  as  data.  In  each  case,  five 
different  performance  indices  were  minimized.  The  Orbiter, 
OTV  and  SOC  orbit  maintenance  fuel  consumption  were 
minimized  individua"' ’  y.  Then  the  combined  (Extended  STS) 
fuel  consumption  and  the  number  of  Shuttle  launches  were 


Table  4.1  Standard  Parameter  Configuration 


SOC  Ballistic  Coefficient  (m^/kg)  =  0.02 
Atmospheric  Scale  Height  (km)  =  30.00 
OTV  Mass  (kg)  =  2270.00 
SOC  and  OTV  Fleet  Mass  (kg)  =  100000.00 
Annual  SOC  Cargo  Mass  (kg/yr)  =  20000.00 
Chemical  OTV  Engine  Specific  Impulse  (sec)  =  455.00 
Chemical  SOC  Engine  Specific  Impulse  (sec)  =  455.00 
ORB  Specific  Impulse  (sec)  =  313.00 
Expended  Mass  at  ETS  (kg)  r  38399.00 
Shuttle  Reference  Mass  (kg)  =  88041.00 
Ion  Thruster  Mass  (kg)  =  51.36 
Ion  Thruster  Thrust  (mN)  =  129.00 
Ion  Thruster  Power  (kw)  =  3.06 
Ion  Thruster  Fuel  Tank  Mass  (kg)  =  12.00 
Ion  Thruster  Specific  Impulse  (sec)  =  2900.00 
Ion  Thruster  Specific  Mass  (kg/kw)  =  10.00 
Ion  OTV  Structural  Mass  (kg)  =  300.00 
Number  of  Ion  Thrusters  (#) 


20.00 


Table  4.2  Satellite  Traffic  Models 


Satellite 

Mission 

^sat 

(km) 

®sat 

^sat 

(deg) 

Mass 

(kg) 

Launch 

Frequency 

(#/YR) 

Test  Missions 

I 

6900 

0.05 

45.0 

20000 

3.0 

II 

40000 

0.00 

0.0 

20000 

3.0 

III 

6500 

0.75 

100.0 

20000 

3.0 

Hypothetical 

Operational 

Missions 

1 

41000 

0.0 

o 

• 

o 

2000 

3.0 

2 

20000 

0.0 

o 

• 

o 

500 

0.5 

3 

6700 

0.0 

28.5 

25000 

0.5 

4 

65000 

0.0 

55.0 

1500 

4.0 

5 

25000 

0.7 

65.0 

1500 

0.5 

6 

12000 

0.0 

90.0 

4500 

3.0 

rainiraized.  In  each  case  SUMT  specified  a  particular 
inclination  and  altitude  for  the  SOC,  which  was  compared 
with  predicted  values  for  that  case.  Table  4.3  shows 
predicted  and  calculated  values  for  all  cases. 

All  of  the  actual  results  either  corresponded  to  the 
predicted  results,  or  would  have  had  no  constraints  been 
present.  As  an  example  the  five  simplified  oases  with  only 
one  satellite  test  mission  (see  Table  4.3)  are  individually 
interpreted.  (All  of  the  predicted  SOC  radii  are 
constrained  within  the  6500  to  7200  km  envelope): 


1.  ETR  Orbiter  Fuel  Minimization--In  this  case 
only  the  Orbiter  fuel  consumption  is  minimized. 
The  predicted  SOC  Radius  is  as  low  as  possible 
(6500  km)  while  inclination  is  as  high  as 
possible.  The  latter  is  a  result  of  the  Orbiter 
being  easiest  to  maneuver  when  its  mass  is  lowest. 
Since  payload  mass  decreases  with  increasing 
inclination,  the  inclination  is  driven  high  sub¬ 
ject  to  the  payload  mass  being  greater  than  zero. 
As  anticipated,  SUMT  drove  the  altitude  to  6500  km 
and  the  inclination  to  115.6°  (payload  mass  infi¬ 
nitesimally  greater  than  zero).  Because  this 
problem  is  independent  of  the  number  of  test  mis¬ 
sions  considered  the  results  are  not  repeated  for 
the  two  and  three  satellite  test  missions  listed 
in  Table  4.3. 

2.  OTV  Fuel  Minimization--In  this  case  only  the 
OTV  fuel  consumption  is  minimized.  Consequently, 
the  predicted  SOC  inclination  and  radius  are  iden¬ 
tical  to  the  satellite  test  missions  inclination 
and  radius.  As  anticipated,  SUMT  drove  the 
inclination  to  45°  and  the  radius  as  high  as 
possible  subject  to  the  Orbiter  fuel  consumption 
being  less  than  10830  kg.  The  latter  constraint 
was  again  infinitesimally  small  per  the 
prediction . 

3.  SOC  Fuel  Minimization — The  SOC  orbit  mainten¬ 
ance  fuel  consumption  is  again  independent  of  the 
number  of  test  missions  considered  and  is  there¬ 
fore  listed  only  once  in  Table  4.3.  The  predicted 
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SOC  radius  is  as  high  above  the  atmosphere  as 
possible  (7200  km)  while  the  inclination  is  as  low 
as  possible.  The  latter  is  due  to  the  atmospheric 
rotation  factor  of  Section  3.2  which  results  in 
slightly  less  drag  at  lower  inclinations.  Again, 
SUMT  drove  the  inclination  to  its  lowest  possible 
value  of  28.5°  and  the  radius  as  high  as  the 
Orbiter  fuel  constraint  would  allow  (an  uncon¬ 
strained  problem  was  also  run  which  drove  the 
radius  to  the  predicted  7200  km). 

4.  ETR  Extended  STS  Fuel  Minimization--In  this 
case  the  Orbiter,  OTV  and  SOC  orbit  maintenance 
fuel  consumption  are  minimized  simultaneously. 
Because  of  the  massive  amount  of  fuel  required  for 
the  OTV  to  change  inclinations,  OTV  fuel 
consumption  is  the  dominant  factor  in  the  minimi¬ 
zation  problem.  Consequently  the  predicted  SOC 
inclination  and  radius  will  be  very  near,  but  not 
identical  to  the  satellite  test  missions  inclina¬ 
tion  and  altitude.  As  anticipated,  SUMT  drives 
the  inclination  to  somewhat  less  than  the 
predicted  45°  inclination,  and  the  radius  towards 
6900  km.  The  actual  radius  was  limited  to  6757  km 
due  to  the  Orbiter  fuel  constraint. 

5.  ETR  Extended  STS  Launch  Minimization--This 
problem  involves  all  the  same  components  as  the 
fuel  minimization  problem  above.  However,  due  to 
the  payload  advantage  of  launching  to  lower  incli¬ 
nations  this  model  will  tend  to  drive  the  SOC 
inclination  somewhat  lower  than  the  previous  case. 
The  radius  will  again  be  driven  to  the  satellite 
test  mission  radius.  As  anticipated,  both  pre¬ 
dictions  were  verified  with  the  Orbiter  fuel  con¬ 
sumption  again  constraining  the  radius. 


In  addition  to  a  single  satellite  test  mission,  two  and 
three  test  missions  were  used  simultaneously  for  data  as 
depicted  in  Table  4.3.  In  both  cases  predictions  can  be 
made  only  within  broad  limits;  nonetheless  the  results  again 
validated  the  accuracy  of  the  model. 

The  test  runs  of  Table  4.3  summarize  the  validation 
methodology  used,  but  were  not  the  only  criterion  used  to 
validate  the  model.  In  addition,  numerous  similar  runs  were 


accomplished,  as  well  as  a  vast  number  of  hand  calculations 
on  every  component  of  the  model.  In  summary,  the  derived 
models  appear  to  be  extremely  accurate  although  an 
independent  verification  by  other  investigators  is  desirable 


4.2  Models  I  and  II — Fuel  4  Launch  Minimization  Solutions 

To  analyze  each  model  several  different  scenarios  were 
compared.  Scenario  A  contains  one  SOC  supplied  by  Shuttle 
launches  from  either  the  ETR  or  WTR,  and  servicing  all  seven 
operational  missions  in  the  traffic  model  of  Table  4.2 
Scenario  B  contains  two  SOCs  with  the  first  supplied  by  ETR 
and  the  second  by  WTR  Shuttle  launches.  The  ETR  SOC  ser¬ 
vices  satellite  missions  1  to  3  while  the  WTR  SOC  services 
missions  4  to  7.  Scenario  C  is  identical  to  B  except  the 
ETR  SOC  services  missions  1  to  5  while  the  WTR  SOC  services 
missions  6  to  7.  Scenario  D  contains  a  total  of  three  SOCs. 
The  first  services  missions  1  to  3  and  is  supplied  by  ETR 
Shuttle  launches.  The  second  services  missions  4  to  5  and 
is  supplied  by  ETR  or  WTR  launches.  The  last  SOC  services 
missions  6  to  7  and  is  supplied  by  WTR  launches.  This 
information  is  summarized  in  Tables  4.4a  and  4.4b  where  the 
combined  fuel  consumption  and  the  number  of  Shuttle  launches 
were  respectively  minimized. 

Both  problem  variables  (SOC  inclination  and  radius), 
fuel  mass,  number  of  Shuttle  launches  and  Orbiter  payload 
mass  are  shown  for  each  run  of  Tables  4.4a  and  4.4b.  The 
total  minimum  fuel  consumption  is  then  calculated  for  each 
scenario  of  Table  4.4a.  Similarly,  the  total  minimum  number 
of  Shuttle  launches  is  calculated  in  Table  4.4b.  All  of  the 
runs  were  made  with  the  standard  parameters  of  Table  4.1. 
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Comparing  the  results  of  the  two  tables  clearly 
demonstrates  the  similarities  between  the  two  performance 
indices.  The  primary  difference  between  minimizing  the 
number  of  launches  instead  of  fuel  consumption  is  that  the 
SOC  inclination  is  driven  to  a  lower  value  as  discussed  in 
Section  4.1.  Due  to  the  similarities  of  the  two  models  and 
the  overwhelming  cost  of  Shuttle  launches  (a  recent  General 
Accounting  Office  report  cited  $56  million  per  launch)  the 
launch  minimization  model  is  used  for  all  future 
calculations . 

To  analyze  the  sensitivity  of  the  ETR  Launch 
Minimization  Model  each  of  the  standard  parameters  of  Table 
4.1  were  varied.  For  the  first  eight  the  parameters  were 
varied  up  and  down  by  50%.  Due  to  their  greater  sensitivity 
the  last  two  (Expended  Mass  at  ETS  and  Shuttle  Reference 
Mass)  were  only  varied  by  10%.  Table  4.5  lists  each  para¬ 
meter  and  its  percent  variation  from  standard.  The  opti¬ 
mized  values  listed  include  the  SOC  inclination,  radius, 
combined  fuel  mass  and  total  number  of  Shuttle  launches. 

The  most  sensitive  parameters  are  the  OTV  specific 
impulse,  OTV  mass,  Orbiter  specific  impulse,  expended  mass 
at  ETS  and  the  Shuttle  reference  mass.  Varying  the  atmos¬ 
pheric  scale  height  also  has  a  significant  effect  on  the 
model.  This  is  especially  noteworthy  since  atmospheric 
density  (and  therefore  scale  height)  varies  drastically  with 
solar  activity.  Parameters  which  appear  to  have  a  minor 
impact  on  the  model  include  the  ballistic  coefficient,  SOC 
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Mass,  annual  SOC  cargo  mass  and  the  specific  impulse  of  the 
SOC  motors  used  for  stationkeeping. 


Three  other  quantities  worth  varying  are  the  two 
problem  variables  (inclination  and  radius)  and  the  launch 
frequencies  in  the  Satellite  Traffic  Model.  Table  4.6  shows 
the  impact  on  the  Extended  STS  fuel  consumption,  Shuttle 
launches  and  Shuttle  payload  mass  when  the  SOC  inclination 
and  radius  are  varied.  Initially  the  SOC  radius  is  left  at 
its  optimum  value  while  the  inclination  is  varied  from  30° 
to  100°.  The  inclination  is  then  set  at  its  optimum  value 
while  the  radius  is  varied  from  6500  to  7100  km.  Both  runs 
are  done  using  ETR  scenario  A.  The  results  clearly  indicate 
that  the  minimum  number  of  Shuttle  launches  is  highly  sensi¬ 
tive  to  the  SOC  inclination  and  altitude.  This  sensitivity 
lends  credibility  to  the  concept  of  designing  the  SOC  such 
that  its  inclination  and  altitude  can  be  changed  over  time. 
For  example  an  altitude  increase  may  be  desirable  during 
peak  solar  activity  while  an  altitude  or  inclination  change 
may  be  needed  as  the  Satellite  Traffic  Model  grows. 

Table  4.7  depicts  optimized  values  for  various  satel¬ 
lite  launch  frequencies.  The  launch  frequencies  of  Table 
4.1  are  varied  up  and  down  by  25%  and  50%.  The  optimized 
values  include  SOC  inclination,  SOC  radius.  Extended  STS 
fuel  mass  and  the  number  of  Shuttle  launches.  Again,  the 
number  of  Shuttle  launches  required  is  especially  sensitive 
to  any  variation  of  the  Satellite  Traffic  Model,  especially 
launch  frequencies. 
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4.3  Model  III--Launch  Minimization  Solutions  with 
OTV/Shuttle  Rendezvous  in  LEO 

Only  the  first  three  scenarios  (A,  B  and  C)  were  used 
to  evaluate  the  utility  of  an  OTV/Shuttle  rendezvous  in  LEO. 
A  previous  study  indicated  that  such  a  strategy  would 
probably  not  be  worthwhile  (Ref  16,  p.  1-5).  In  contrast, 
the  results  of  Table  4.8  indicate  that  significant  savings 
are  realizable  if  the  Shuttle  traffic  is  sufficiently  high. 

Comparing  Table  4.8  to  Table  4.4b  demonstrates  several 
surprising  results.  For  Scenario  A  the  OTV/Shuttle 
rendezvous  strategy  reduces  the  required  number  of  Shuttle 
launches  by  20%  (47.70  to  38.27).  For  Scenarios  B  and  C  the 
corresponding  launch  savings  are  33%  (42.10  to  28.32)  and 
38%  (37.82  to  23.30)  respectively.  In  terms  of  payload  mass 
delivered  to  the  SOC,  these  figures  roughy  correspond  to  a 
20%  payload  increase  for  ETR  launches  and  a  10%  increase  for 
WTR  launches. 

Although  the  OTV/Shuttle  Rendezvous  strategy  appears 
very  promising,  several  mitigating  factors  do  exist.  The 
Shuttle  cargo  and  crew  traveling  to  the  SOC  after  rendezvous 
would  have  to  be  placed  in  modular  containers  that  the  OTV 
could  dock  with.  Such  containers  would  have  to  be  built  and 
rated  for  manned  flight.  Another  potential  problem  lies  in 
the  WTR  MECO  Mass  calculations.  Appendix  A  calculates  only 
the  MECO  mass  above  70°  inclination.  A  third  order  equation 
specifying  the  MECO  mass  as  a  function  of  inclination  is 
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Table  4.8  Model  III  Launch  Minimization — 

OTV/Shuttle  LEO  Rendezvous  Scenarios 
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then  derived  by  linear  regression.  Consequently,  the  data 
for  56°  to  70°  inclinations  is  somewhat  suspect.  Since  one 
of  the  WTR  launches  of  Table  4.8  is  within  this  range  it  is 
also  suspect.  Finally,  the  extent  of  the  savings  will 
depend  on  the  traffic  model  used.  Traffic  models  which 
place  the  optimum  SOC  altitude  in  a  lower  Earth  orbit  are 
not  as  greatly  influenced  by  the  OTV/Shuttle  rendezvous 
strategy  as  higher  altitude  SOCs. 

Several  of  the  parameters  of  Table  4.5  were  again 
varied  within  Model  III.  The  OTV  mass,  OTV  specific 
impulse,  Orbiter  specific  impulse  and  Shuttle  reference  mass 
are  all  varied  up  and  down  by  50t  or  1 Ot  as  indicated  in 
Table  4.9.  As  before  the  OTV  mass,  OTV  specific  impulse  and 
the  Shuttle  reference  mass  were  all  critical  parameters  that 
again  had  a  major  impact  on  the  calculated  number  of  Shuttle 
launches.  Unlike  the  previous  case,  an  increase  in  the 
Orbiter  specific  impulse  resulted  in  only  a  marginal 
decrease  in  the  number  of  Shuttle  launches.  This  outcome 


was  predictable  since  the  Orbiter,  being  deployed  to  a  lower 
altitude,  doesn't  consume  as  much  propellant. 


58 


ft 


I 


to 

x:  u 

\o 

on 

o 

on 

CO  o  >. 

CM 

m 

CM 

oo 

in 

ST 

r— 

in 

H  C  V, 

• 

• 

• 

• 

• 

• 

• 

• 

• 

CO  3  =ite 

CO 

cn 

ON 

CM 

in 

in 

in 

(0  W 

on 

CM 

T— 

in 

m 

CM 

n 

P 

u 

<0 

■p 

(U 

E 

(0 

L. 

(Q 

a. 

P  (0  L. 

CT> 

cy> 

CO 

in 

CM 

oo 

on 

4)  01  >< 

CO 

m 

m 

in 

m 

CM 

on 

1  M 

3  (0  V. 

<3\ 

in 

vO 

CM 

in 

m 

f— 

vO 

C- 

1  3 

tb  Z  00 

in 

rv^ 

cas 

oo 

CM 

on 

=r 

C  O 

oo 

on 

m 

T“ 

vO 

r- 

vO 

o  > 

tn 

o 

m 

f— 

O 

•H  N 

CO 

t— 

<a> 

■P  (U 

(U 

ra  -o 

3 

N  C 

rH 

•H  0) 

(0 

E  OS 

> 

•H 

01 

o 

00 

oo 

in 

CM 

a% 

t- 

in 

CM 

c  « 

T3 

CJ  3^ 

CM 

lO 

0- 

CM 

CM 

CO 

*H  rH 

<U 

O  -H  E 

• 

• 

• 

• 

• 

• 

• 

• 

• 

ac  4-> 

N 

CO  -o 

cy> 

vO 

lO 

CM 

CM 

t- 

VMI 

o 

■p 

•fH 

fl  ^ 

on 

in 

CM 

vO 

vO 

on 

on 

ST 

on 

£i  3 

E 

OS 

t>- 

in 

iM 

t- 

C>- 

u  x: 

•H 

vO 

lO 

vO 

vO 

VO 

VO 

VO 

•o 

VO 

C  Vi 

p 

3  V. 

o. 

«  > 

o 

J  H 

o 

as 

H  Vi 

u  o 

O  P  00 

o 

o 

O 

ca> 

CM 

o 

O 

o 

m 

O  O  41 

o 

o 

O 

CO 

on 

o 

o 

o 

M  C 

CO  z  -a 

« 

• 

• 

• 

• 

• 

• 

• 

• 

l-l  O 

M  w 

VO 

in 

t— 

m 

M  .H 

m 

tn 

in 

in 

in 

in 

in 

in 

in 

rH  (0 

d)  'H 

•o  l. 

O  CO 

•o 

z  > 

L. 

C  (Q 

P  O  T3 

C  -H  C 

4)  P  (0 

41 

• 

u  ra  P 

C 

o 

o 

o 

o 

o 

o 

o 

o 

i-  P  CO 

o 

in 

in 

in 

in 

in 

in 

(U  V 

z 

1 

+ 

1 

+ 

1 

+ 

1 

+ 

0) 

ou  (0  E 

rH 

>  o 

£> 

L. 

CO 

<M 

H 

ci 

Cl 

L. 

•3  t. 

3 

CO 

•H 

in 

•H 

on 

41 

V  4) 

1.  01 

01  --S 

Cm 

in 

Cm 

<M  CQ 

(0  p 

(0  41 

<0  o 

P  41  ^ 

p  41  on 

41  01  »- 

■O  4; 

3  3 

ac 

Cl  01  II 

O  01  II 

OS  <0  sr 

c  E 

C  P 

CM 

41  p  > 

4)  r- 

^  p 

z  o 

(0  CO 

(tl  (Q 

>  CM 

a  3  p 

a  : 

1  L. 

41 

oo 

P  L. 

P  > 

E-I  II 

CO  a  o 

CO  a  o 

P  4100 

CO  m 

CO 

o 

> 

E  Qu 

E  a 

P  c 

II 

(X, 

p 

S»  M  CO 

CQ  M  CO 

P  c 

OS 

o 

H 

M 

OS 

M 

3  41  Z 

a 

■ 

O 

o 

'•Ml' 

ss 

>MI<' 

1 

Vi 

59 


J 


4.4  Model  IV — Launch  Minimization  Solutions  with  an 
Ion-Propelled  OTV 

The  use  of  an  ion  propelled  OTV  to  deploy  and  retrieve 
satellites  can  drastically  reduce  the  number  of  Shuttle 
launches  required  to  service  satellites  within  the  Traffic 
Model.  Indeed,  Tables  4.5  and  4.9  clearly  demonstrate  that 
a  50i  increase  in  the  OTV  specific  impulse  dramatically 
reduces  the  number  of  Shuttle  launches.  With  an  ion- 
propelled  OTV  specific  impulse  can  easily  be  increased  by 
1000?  and  more.  Consequently,  vast  savings  in  fuel  and 
Shuttle  launches  are  possible  if  the  user  is  willing  to 
accept  much  longer  deployment  times. 

Table  4.10  lists  optimized  values  for  a  variety  of 
different  ipecific  impulses.  The  optimized  values  include 
the  SOC  inclination,  SOC  radius,  total  fuel  consumption, 
number  of  Shuttle  launches  and  the  average  mission  duration 
when  deploying  satellites.  The  2900  and  9000  second 
impulses  were  chosen  since  ion  systems  with  those 
capabilities  exist  or  are  being  investigated.  Both  are  30 
cm  mercury  ion  thrusters  being  developed  by  NASA 
(Ref  13,  p.  62).  The  SOC  stationkeeping  thrusters  were 
generally  assigned  the  same  specific  impulse  as  the  OTV 
thrusters.  The  rationale  is  that  if  such  thrusters  are 
available  they  might  as  well  be  used  for  SOC  orbit 
maintenance  and  orbit  changes. 


Table  ^.10  Model  IV  ETR  Launch  Minimizatlon-- 

Variation  of  OTV/SOC  Specific  Impulse 


With  the  use  of  ion  thrusters  the  required  number  ' 
annual  Shuttle  launches  drops  drastically.  As  a  rough  guess 
it  would  probably  require  ten  or  more  annual  launches  for 
the  Shuttle  to  directly  deploy  the  satellites  in  the  Traffic 
Model.  In  the  extreme  cases  where  specific  impulse  is  set 
at  20000  and  30000  seconds,  the  number  of  Shuttle  launches 
is  easily  less  than  half  that  required  for  direct  Shuttle 
insertion.  All  of  these  runs  were  accomplished  with  the  OTV 
mass  set  at  2270  kg. 

Although  such  an  OTV  mass  is  certainly  achievable,  it 
may  be  desirable  to  add  additional  ion  thrusters  that  would 
increase  the  mass,  but  decrease  the  mission  durations.  To 
evaluate  the  more  general  case  where  the  OTV  mass  increases 
with  the  number  of  ion  thrusters  required  the  use  of  Eq. 
3.30: 

^otv  “  ^^^thr  ^sp  ^  ^wr  ^ft^  ^str 

The  parameters  used  in  this  equation  were  for  a  NASA  30  cm 
Mercury  ion  thruster  which  has  been  developed  and  is  being 
refined  (Ref  13,  p.  62).  The  actual  parameter  values  are 
listed  in  Table  4.1.  Also  in  Table  4.1  is  the  structural 
mass  (Mstr^>  number  of  ion  thruster  and  the  specific 

mass  of  the  ion  system.  The  latter  is  simply  the  mass  of 
the  power  system  and  conditioning  equipment  required  to 
generate  a  1  Kw  output.  The  value  of  10  Kg/Kw  reflects  the 
current  level  of  technology.  The  resulting  OTV  mass  is  2179 
Kg,  slightly  less  than  the  previous  case.  For  consistencies 
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sake  the  SOC  orbit  maintenance  thrusters  are  assumed  to  be 
the  same  as  on  the  OTV  with  a  specific  impulse  of  2900  sec. 

Scenarios  A,  B  and  C  were  used  to  evaluate  this  new 
model  with  Table  4.11  summarizing  the  results.  The  ETR 
Scenario  A  required  only  9.08  annual  Shuttle  launches  to 
service  the  Satellite  Traffic  Model.  Scenario  C  with  one 
SOC  inclined  at  28.5°  and  the  other  at  56°  required  about 
2.4  additional  launches  per  year.  However,  it  has  the 
benefit  of  two  operating  stations  (also  the  extra  expense) 
and  somewhat  smaller  mission  durations  when  deploying  polar 
orbiting  satellites.  Scenario  B  is  nearly  identical  to  C 
except  the  average  OTV  mission  duration  is  240  vs  34?  days. 
The  smaller  mission  duration  is  advantageous  since  the 
number  of  OTVs  required  to  service  the  satellite  Traffic 
Model  is  proportionally  smaller. 

As  in  section  4.3  several  of  the  key  parameters  were 
again  varied.  Both  the  Shuttle  reference  mass  and  the  OMS 
specific  impulse  were  varied  as  before.  However,  the  OTV 
mass  and  specific  impulse  were  not  varied  directly. 
Instead,  several  of  the  parameters  used  to  calculate  OTV 
mass  were  varied  up  and  down  by  50%.  Table  4.12  summarizes 
the  results  of  the  sensitivity  analysis.  As  anticipated, 
the  Shuttle  reference  mass  and  Orbiter  specific  impulse  had 
a  pronounced  effect  on  the  number  of  Shuttle  launches. 

The  remaining  ion  engine  parameters  had  a  relatively 
minor  impact  on  the  number  of  Shuttle  launches.  However, 
they  had  a  major  impact  on  the  average  OTV  mission  durations. 
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The  ion  engine  thrust  was  by  far  the  most  significant 
parameter.  A  50%  increase  in  thrust  caused  the  average 
mission  duration  to  decrease  from  415  to  276  days.  A  50% 
variation  in  the  number  of  ion  thrusters,  thruster  mass, 
thruster  power  and  the  specific  mass  had  a  smaller,  but 
still  significant  impact  on  the  mission  durations.  The 
respective  average  durations  were  reduced  to  346,  358,  381 
and  3  81  days. 

Although  the  OTV  specific  impulse  was  not  varied  in 
Table  4.12,  the  runs  of  Table  4.10  adequately  measure  the 
impact  on  Shuttle  launches  and  mission  durations.  Both  are 
significantly  reduced.  For  example,  varying  the  specific 
impulse  from  2900  to  5000  seconds  reduced  the  number  of 
annual  launches  from  9.17  to  6.83  and  the  mission  duration 
from  425  to  362  days. 


a 


4.5  Model  V — Launch  Minimization  Solutions  with  Chemical 
OTV/STS  Rendezvous  in  LEO  and  Ion  OTV  for 
Satellite  Deployment 

The  final  model  minimized  with  SUMT  is  a  combination  of 
Models  III  and  IV.  A  chemically  propelled  OTV  is  used  to 
rendezvous  with  the  Shuttle  and  retrieve  modularized  Shuttle 
payloads.  The  rendezvous  orbit  has  an  altitude  equivalent 
to  the  ballistic  Shuttle  apogee  after  MECO  (160  km).  Conse¬ 
quently,  only  one  Shuttle  OMS  burn  is  required  to  achieve 
orbit  and  one  to  de-orbit.  The  rendezvous  inclination  is 
set  to  the  same  value  as  the  SOC  inclination.  This  simpli¬ 
fication  is  justified  by  the  Model  III  results  illustrated 
in  Table  4.8.  In  that  case  the  SOC  and  rendezvous  inclina¬ 
tions  were  separate  variables  which  the  optimization  proce¬ 
dure  drove  to  nearly  identical  values. 

Another  fleet  of  ion-propelled  OTVs  is  used  to  deploy 
the  satellites  in  the  traffic  model.  Although  the  ion- 
propelled  OTV  is  much  more  fuel  efficient  than  its  chemical 
counterpart,  it  probably  could  not  be  safely  used  for  a  LEO 
Shuttle  rendezvous.  Atmospheric  drag  at  lower  altitudes  may 
well  be  too  great  for  an  ion  thrust  system  to  overcome, 
especially  during  peak  solar  activity.  Further,  a 
chemically  propelled  OTV  would  probably  be  required  for  all 
manned  missions.  The  ion  system  is  too  slow  for  either  a 
manned  rendezvous  with  the  Shuttle  in  LEO  or  any  other 
vehicle  (possibly  to  repair  a  satellitt). 


Another  noteworthy  feature  of  this  model  is  the 
quantity  of  Orbiter  OMS  fuel  consumed.  Even  when  deployed 
to  such  a  low  altitude  the  Orbiter  fuel  consumption  is  over 
4900  Kg  for  ETR  launches  (over  8400  Kg  for  WTR),  still  a 
significant  fraction  of  the  10830  Kg  fuel  tank  capacity.  It 
is  probable  that  OMS  fuel  consumption  could  be  significantly 
reduced  by  using  less  conservative  ascent  and  de-orbit  aVs, 
and  by  altering  the  standard  MECO  conditions  to  better 
utilize  residual  fuel  within  the  External  Tank.  Appendix  A 
delineates  the  crude  mission  planning  equations  used  to 
calculate  ascent  and  de-orbit  aVs.  The  ETR  ascent  equation 
gives  a  aV  roughly  ^0%  higher  than  the  actual  requirement 
while  the  deorbit  equation  appears  even  more  conservative. 
Tailoring  these  equations  would  reduce  fuel  consumption 
slightly  while  launching  to  different  standard  MECO 
conditions  may  result  in  dramatic  fuel  savings.  Both 
options  should  be  further  evaluated  in  the  light  of  data 
from  future  operational  Shuttle  launches. 

Scenarios  A,  B  and  C  were  again  used  to  evaluate 
Model  V.  Table  4.13  summarizes  the  results.  The  ETR  Scena¬ 
rio  A  required  7.8  annual  shuttle  launches,  a  14%  reduction 
from  Model  IV.  Scenario  B  required  9.11  and  Scenario  C 
required  9.08  annual  launches,  both  down  21%  from  Model  IV. 
For  Scernarios  A,  B  and  C  the  average  OTV  mission  TOF  is 
411,  243  and  347  days  respectively.  Scenario  B  is  again 
significantly  faster  than  the  other  two.  However,  at  $56 
million  per  launch  scenario  B  is  also  somewhat  more 
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Ion  Propelled  OTV  and  Shuttle 


e  n 

iH  3  0) 

rascox:^- 
-H  E-*  O  ^ 
O  C  W  c  ^ 

H  -H  3  =Sfe 


U|  0]  c 


OT  C  M  t- 
H  3  <0  >. 
CO  nj  £  ^ — 
J  o  ^ 


CO  ™ 

Om  s_x 


CJ  bO 

o  ^  <u 

CO  ^  T3 


^  'r-i  .  s-/ 
=!*i  rH  -j:  -H 

0  0)”$- 
C/1  4J  -r-t 

«  ac  CO 

CO 


u  c  o 

O  0)  -H 
CO  CJ  u 
CO  (0 


o 

oo  1 

GO 

O  1 

• 

• 

• 

t- 

CT' 

ON 

=T 

OO 

t- 

:T 

=T 

CM 

OO 

o 

CM 

c- 

CM 

MD 

GO 

*“ 

ON 

vO 

t— 

Ov 

*- 

CM 

vO 

CM 

ON 

ON 

CO 

vO 

T— 

CJN 

C~- 

m 

O 

o 

O 

OO 

CM 

CO 

oo 

on 

CTi 

ON 

OO 

o 

oo 

C7V 

CM 

CM 

O 

CM 

oo 

1 

CM 

O 

in 

OO 

oo 

O 

CM 

m 

o 

in 

o 

=r 

t*- 

OO 

CM 

CM 

CM 

O 

oo 

m 

on 

CM 

o 

00 

ro 

cr> 

OO 

iT 

OO 

m 

m 

vO 

CM 

OO 

in 

m 

CM 

CM 

OO 

CM 

o 

c~- 

oo 

GO 

oo 

VO 

VO 

VO 

VO 

vO 

vO 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

O 

o 

• 

• 

* 

• 

• 

o 

vO 

cr> 

VO 

cr> 

vO 

cn 

in 

CM 

m 

CM 

in 

cc 

a: 

os 

cc 

oc 

DC 

E-> 

E- 

E-* 

H 

H 

H 

U] 

3 

u 

3 

UJ 

oo 

in 

1 

I 

1 

1 

1 

t“ 

VO 

V, 

\ 

CM 

CM 

=Ste 

o 

o 

u 

CJ 

o 

o 

o 

O 

o 

CO 

CO 

CO 

CO 

CO 

< 

CQ 

u 

70 


expensive  to  operate,  around  $73  million  in  annual  launch 
costs.  Part  of  this  expense  could  be  recouped  if  the  56° 
inclined  SOC  were  supplied  by  ETR  Shuttle  launches  instead 
of  the  Model's  WTR  launches.  The  ETR  launches  would  have  a 
significantly  greater  payload  capability  than  their  WTR 
counterparts . 


4.6  Model  VI--Model  V  Launch  Minimization  Solutions  versus 
Average  OTV  Flight  Times 

From  the  perspective  of  reducing  launch  costs  the  Model 
V  results  of  the  last  section  appear  very  promising. 
Another  equally  important  performance  parameter  which  must 
be  considered  is  the  ion  propelled  OTV  flight  times.  Long 
flight  times  to  the  satellite's  mission  orbit  have  an 
adverse  impact  on  both  functional  reliability  and 
operational  utility.  Consequently,  from  an  operational 
perspective  it  would  be  ideal  if  both  the  flight  time  and 
number  of  Shuttle  launches  could  be  minimized. 
Unfortunately,  to  a  great  extent  there  is  a  quid  pro  quo 
trade  off  between  the  two  performance  indices.  Flight  times 
can  be  reduced  by  adding  ion  thrusters  to  the  OTV,  but  the 
increased  OTV  and  fuel  mass  requirements  eventually  drives 
up  the  number  of  Shuttle  launches. 

Implementation  of  a  vector  optimization  process  defines 
the  exact  trade  off  between  the  two  indices.  Figure  4.1 
depicts  the  efficient  frontier  on  which  the  Extended  STS 
should  be  operated.  It  was  calculated  using  the  standard 
parameters  of  Table  4.1  with  one  scenario  A  SOC  servicing 
the  entire  Traffic  Model  of  Table  4.2.  Each  point  on  the 
curve  has  associated  with  it  an  optimal  SOC  altitude,  incli¬ 
nation  and  number  of  ion  thrusters  on  the  OTV.  The 
hyperbolic  configuration  demonstrates  that  up  to  the  "knee" 


age  Ion  OTV  Round  Trip  Flight  Times  (days) 

Scenario  A  Launches  vs  OTV  TOF  for  Standard  Parameters 


of  the  curve  average  OTV  flight  times  can  be  dramatically 
reduced  without  substantially  increasing  Shuttle  launches. 
Below  the  knee  annual  Shuttle  launch  requirements  increase 
exponentially.  This  flight  time  reduction  is  primarily  due 
to  a  steady  increase  in  the  number  of  ion  thrusters  on  the 
OTV.  The  exponential  increase  in  Shuttle  launches  for  lower 
flight  times  is  a  function  of  increasing  OTV  mass  and  fuel 
consumption.  It  corresponds  to  a  surprisingly  steep 
exponential  increase  in  fuel  consumption  which  may  be 
partially  due  to  operating  the  ion  engines  at  a  less  than 
optimal  exhaust  velocity.  For  low  thrust  systems  the 
maximum  payload  ratio  for  a  given  AV  is  always  associated 
with  a  particular  exhaust  velocity  (which  corresponds  to 
specific  impulse).  Any  variation  from  the  optimal  exhaust 
velocity  results  in  an  exponential  degradation  of  the 
payload  ratio.  Consequently,  each  mission  in  the  Satellite 
Traffic  Model  has  an  optimal  OTV  specific  impulse  associated 
with  it.  Individually  adjusting  the  ion  engine  thrust  for 
each  mission  may  allow  the  user  to  further  reduce  annual 
Shuttle  launches  or  even  deployment  times. 

Although  better  modeling  may  provide  a  refined 
efficient  operating  frontier,  Fig.  4.1  is  adequate  for  the 
purposes  of  this  study.  Indeed,  its  major  drawback  is  that 
it  is  probably  too  conservative.  To  eliminate  some  of  the 
conservatism  the  standard  parameters  and  Traffic  Model  of 
Tables  4.1  and  ^*.2  were  reviewed  and  assigned  more  realistic 
values  where  appropriate.  Table  4.14  depicts  the  changes  in 
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Traffic  Model 


nine  parameters  for  three  different  levels  of  technology. 
The  first  level  is  currently  achievable  and  varies  only  the 
Shuttle  reference  mass  and  Traffic  Model.  The  former  should 
be  easy  to  achieve  by  fine  tuning  the  STS  flight  profile  and 
mass  configuration  assumptions  in  Tables  3.1  and  3.2.  The 
latter  launch  frequencies  were  reduced  by  25%  (see  the  Low 
Traffic  Model  of  Table  4.7)  to  bring  the  number  of  annual 
Shuttle  launches  in  line  with  more  realistic  projections  of 
military  traffic. 

Level  II  technology  incorporates  only  a  slight 
upgrading  of  the  30  cm  Hg  ion  engine.  The  new  parameters 
were  measured  in  NASA  Lewis  Research  Center  experiments 
(Ref.  17,  Table  3),  and  can  very  probably  be  achieved  in  the 
1990  timeframe.  The  projected  specific  mass  of  5  Kg/Kw  may 
be  slightly  optimistic,  but  it  does  appear  to  be  achievable 
in  the  near  term  (Ref.  18).  The  last  technology  level  is  a 


projection  based  on  current  trends  in  ion  engine  research. 
Reducing  thruster  mass  and  specific  mass  are  both  high 
priority  items.  Ion  engine  thrust,  power  requirements  and 
specific  impulse  are  essentially  linearly  related  over  a 
given  operating  regime  (Ref.  19,  p  5-90).  Consequently, 
these  parameters  are  accordingly  varied  in  Table  4.14.  The 
Shuttle  reference  mass  and  specific  mass  are  again  reduced. 
The  former  increases  the  modeled  Shuttle  payload  mass,  and 
the  introduction  of  light  weight  filament  wound  SRM  casings 
will  essentially  accomplish  the  same.  Again,  the  projected 
specific  mass  may  be  somewhat  optimistic,  but  all  of  the 


6 


projected  parameters  are  realistic  estimates  for  the  2000 
timeframe . 

Applying  vector  optimizaton  to  the  three  technology 
levels  yields  a  consistent  reduction  in  annual  Shuttle 
launches  and  OTV  flight  times.  Figure  4.2  graphically 
illustrates  the  improvement  for  all  three  levels  with  one 
SOC  servicing  the  entire  Low  Traffic  Model  (scenario  A). 
The  three  curves  are  depicted  for  ETR  Shuttle  luanches  while 
the  correspnding  WTR  Shuttle  launch  curve  is  also  shown  for 
Level  I  technology.  The  WTR  supplied  SOC  is  marginally  less 
efficient  than  its  ETR  counterpart.  However,  the  lack  of 
accurate  MECO  mass  data  between  56°  and  70°  casts  some  doubt 
on  the  accuracy  of  this  curve  (the  optimized  SOC  inclina¬ 
tions  are  in  this  region).  In  all  probability  it  should  be 
shifted  upward  somewhat,  thus  increasing  annual  Shuttle 
traffic  and  reflecting  the  decreased  Shuttle  payload 
capability  of  WTR  launches. 

In  an  attempt  to  further  reduce  OTV  flight  times. 
Scenario  B  with  two  SOCs  was  analysed  for  Level  I 
technology.  The  first  SOC  services  satellite  missions  1  to 
3  while  the  second  services  missions  4  to  7.  Because  the 
first  three  missions  are  all  below  28.5°  and  the  last  four 
above  55°  inclinations,  the  variable  SOC  inclinations  were 
logically  set  at  28.5°  and  57.0°  respectively.  These  are 
the  minimum  and  maximum  inclinations  accessible  by  ETR 
Shuttle  launches,  and  specifying  them  as  constants  greatly 
reduced  the  SUMT  optimization  program  convergence  time. 


Average  Ion  OTV  Round  Trip  Flight  Times  (days) 
Fig.  4.2  Annual  Scenario  A  Launches  vs  OTV  TOF 


Figures  4.3  and  4.4  illustrate  the  efficient  operating 
frontiers  for  the  dual  SOC  scenario.  Both  are  supplied  by 
ETR  Shuttle  launches  although  Fig.  4.4  again  includes  the 
corresponding  WTR  Shuttle  resupply  curve.  The  WTR  curve 
reinserts  the  SOC  inclination  as  an  optimized  variable  and 
is  again  probably  somewhat  optimistic. 

If  the  annual  Shuttle  launch  rate  and  the  average  OTV 
mission  duration  are  equally  weighted  performance  indices, 
then  the  ”best"  point  at  which  to  operate  the  Extended  STS 
is  at  the  "knee”  of  the  efficient  operating  frontier.  Table 
4.15  depicts  several  optimized  values  at  the  knee  of  each 
curve  in  Fig.  4.2  through  4.4.  Besides  the  standard  values 
depicted  in  previous  tables,  Table  4.15  includes  the  total 
number  of  ion  propelled  OTVs  required  to  service  the  Traffic 
Model.  This  calculation  assumes  continuous  utilization  of 
every  OTV.  It  should  also  be  remembered  that  the  average 
OTV  flight  times  shown  represent  round  trip  times.  The 
deployment  times  are  of  greater  interest  to  military 
planners  and  are  contained  in  the  computer  runs  cataloged  in 
Appendix  D. 

In  summary,  the  dual  or  single  SOC  scenarios  combined 
with  the  Level  II  technology  available  in  the  near  term 
appears  very  attractive  to  the  military  user.  Shuttle 
launch  rates  can  probably  be  contained  below  cur  ‘ent 
projections  with  reasonable  satellite  deployment  times.  The 
latter  combined  with  the  on-orbit  spares  philosophy  already 
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Average  Ion  OTV  Round  Trip  Flight  Times 
Fig.  4.3  Annual  Scenario  B  Launches  vs  OTV  TOF  for  Missions  1-3 
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Average  Ion  OTV  Round  Trip  Flight  Times 
Fig.  4.4  Annual  Scenario  B  Launches  vs  OYV  TOF  for  Missions  4-7 


Table  4.15  Model  VI  Launch  Minimization  —  Ion  Propelled  OTV 
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III  A  SOC  /M/1-7  ETR  41  .00  7356.41  76.56  40.00 

(2000) 


being  pursued  with  military  programs  can  satisfy  the 
military  requirement  of  rapid  deployment. 


5.0  Summary 


The  ability  of  a  mathematical  model  to  predict  reality 
depends  directly  on  the  accuracy  of  the  underlying 
assumptions.  To  avoid  the  faulty  predictions  of  an  overly 
optimistic  model,  most  of  the  assumptions  in  Model’s  I 
through  V  were  chosen  to  yield  conservatively  high  estimates 
of  annual  Shuttle  traffic.  Conversely,  Model  VI  eliminates 
several  of  the  grosser  assumptions  and  attempts  to  project 
the  impact  of  evolving  technology.  A  brief  review  of  the 
impact  key  assumptions  have  on  all  six  models  follows: 


1.  Satellites  in  the  Traffic  Model  are  currently 
deployed  individually.  In  many  cases  it  will  be 
much  more  efficient  to  deploy  several  satellites 
simultaneously . 

2.  Incorporation  of  satellite  retrieval  and 
repair  missions  requires  additional  OTV  flights  to 
service  satellites  in  the  Traffic  Model.  However, 
the  number  of  new  satellites  deployed  decreases 
proportionally . 

3.  A  constant  atmospheric  scale  height  is  used  to 
calculate  stationkeeping  fuel  consumption.  The 
derived  model  gives  reasonable  results  but  needs 
to  be  verified  by  alternate  methods. 

4.  Altitude  and  inclination  changes  are 
independently  accomplished  with  the  ion  propelled 
OTV.  Combining  the  maneuvers  can  result  in 
significant  fuel  savings  (Models  IV  through  VI 
only ) . 


5.  The  ion  propelled  OTV  fuel  mass  calculation 
assumes  100%  propellant  utilization.  For  compari¬ 
son  the  uprated  30  cm  Hg  engine  of  Table  4.14  has 


a  94.35t  propellant  utilization  rate  (Models  IV 
through  VI  only). 


6.  The  calculated  number  of  ion  propelled  OTVs 
represents  a  lOOt  utilization  rate,  in  reality 
additional  spare  vehicles  would  be  needed  at  the 
SOC  (Model  VI  only) . 


In  addition,  two  key  assumptions  eliminated  by  Model  VI  are 
contained  in  the  first  five  models: 


7.  The  Satellite  Traffic  Model  postulated  repre¬ 
sents  a  high  estimate  of  military  traffic  through 
the  end  of  the  century.  The  Low  Traffic  Model  of 
Table  4.7  and  Model  VI  requires  fewer  Shuttle 
launches  to  support,  but  compares  favorably  with 
more  realistic  projections  of  military  traffic 
requirements . 

8.  The  Orbiter  mass  configuration  and  flight 
profile  assumptions  of  Tables  3.1  and  3.2 
accurately  specify  MECO  mass.  However,  several  of 
the  assumptions  may  be  overly  conservative, 
especially  the  Orbiter  reference  mass  and  the  SRM 
inert  mass.  The  reference  mass  can  probably  be 
substantially  reduced  by  tailoring  resupply 
missions  while  the  use  of  filament  wound  SRM 
casings  will  likewise  reduce  the  inert  mass.  Both 
have  a  dramatic  impact  on  payload  capacity. 


Other  assumptions  made  are  delineated  within  the  text  of 
this  study,  but  are  not  deemed  significant  enough  to  repeat. 

The  consistent  conservatism  exercised  when  making 
assumptions  lends  substantial  confidence  to  the  overall 
results  of  this  study.  Indeed,  extrapolating  the  accuracy 
of  the  MECO  mass  calculations  of  Appendix  A  suggests  that 
the  overall  Shuttle  payload  mass  calculation  is  very  likely 
within  lOt  and  probably  within  5%  of  the  actual  mass. 
Consequently,  the  calculated  number  of  Shuttle  launches 


throughout  this  study  is  probably  also  accurate  to  within 
10%,  or  certainly  calculable  to  that  level  by  fine  tuning 
the  models  and  inputted  data.  In  any  case,  the  accuracy  is 
more  than  sufficient  for  the  gross  requirements  of  this 
transportation  study. 

In  addition  to  the  model  being  accurate,  the 
technological  projections  of  Table  4.14  and  Fig.  4.2  appear 
very  reasonable.  Projecting  technology  into  the  future  is 
always  risky  and  best  avoided.  However,  for  the  sake  of 
those  involved  in  future  planning  it  is  also  essential.  The 
projection  for  1990  timeframe  technology  represents  a  modest 
improvement  in  the  Shuttle  payload  capability  (still  well 
below  the  nominal  65000  lb  figure),  and  an  uprating  of 
current  ion  engine  technology.  The  latter  has  already  been 
demonstrated  in  the  laboratory  and  is  very  probably  feasible 
by  1990.  Projecting  technology  to  2000  is  more  difficult. 
The  projected  Shuttle  capability  at  this  level  is  probably 
close  to  its  maximum  value  without  redesigning  and  refitting 
the  vehicle.  Projected  ion  engine  technology  is  probably 
very  reasonable  with  the  possible  exception  of  the  specific 
mass  which  may  be  slightly  optimistic.  Overall,  both  the 
1990  and  2000  technology  projections  appear  imminently 
feasible  with  current  funding  of  ion  engine  research. 

Comparing  the  Extended  STS  to  current  Shuttle  launch 
costs  is  complicated  by  the  difficulty  of  estimating  the 
number  of  STS  launches  required  to  deploy  the  postulated 
Satellite  Traffic  Model.  The  NASA  STS  Flight  Assignment 


Baseline  projects  that  42  of  the  first  60  operational 
Shuttle  launches  will  be  directly  inserted  to  inclinations 
less  than  57°.  The  total  number  of  listed  payloads  is  106 
thus  requiring  about  2.5  satellite  payloads  per  ETR  Shuttle 
launch  (Ref.  20).  Applying  this  figure  to  the  first  five 
missions  of  the  Traffic  Model  in  Table  4.2  yields  at  least 
three  (3.4)  annual  launches  to  support  those  missions.  The 
more  massive  satellites  at  polar  inclinations  will  probably 
require  one  Shuttle  launch  per  satellite.  Adding  the 
results  gives  a  crude  figure  of  eleven  annual  Shuttle 
launches  to  deploy  the  Satellite  Traffic  Model.  The  corres¬ 
ponding  figure  for  the  Low  Traffic  Model  utilized  in  the 
Model  VI  results  (Section  4.6)  is  eight  annual  Shuttle 
launches. 

Comparing  these  figures  to  those  derived  by  the  various 
models  of  section  4.0  is  enlightening.  The  SOC  utilized  in 
conjunction  with  the  chemical  OTV  of  Model  II  appears 
completely  uneconomical  except  when  used  to  deploy  only 
satellite  missions  one  to  three.  These  satellites  are 
deployed  at  inclinations  ranging  from  equatorial  to  28.5° 
and  altitudes  from  Geosynchronous  to  low  Earth  orbit. 
Although  Model  II  still  requires  a  higher  Shuttle  launch 
rate  (5.75  launches/yr)  it  must  be  remembered  that  this  rate 
includes  deployment  and  resupply  of  the  SOC.  In  Model  III 
the  launch  rate  for  missions  one  to  three  is  reduced  by 
about  21%  (4.56  launches/yr)  by  having  the  OTV  rendezvous 
with  the  Shuttle  at  a  160  km  altitude.  This  compares  even 
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more  favorably  with  direct  Shuttle  insertion  of  satellites. 
These  two  launch  rates  correspond  to  the  Scenario  B  rates  in 
Tables  4.4b  and  4.8.  Both  are  slightly  conservative  esti¬ 
mates  of  annual  Shuttle  launches  due  to  a  minor  sign  error 
in  the  MECO  mass  calculation  (Models  I  to  IV  only--see 
Appendix  C).  The  sign  errors  were  corrected  in  the  computer 
runs  of  Appendix  D,  but  for  consistencies  sake  are  not 
incorporated  in  the  results  of  Section  4.0.  The  corrected 
launch  rates  corresponding  to  the  5.75  and  4.56  values  are 
5.70  and  4.51  launches  per  year. 

Model  IV  utilizes  an  ion  propelled  OTV  to  deploy 
satellites  while  Model  V  combines  the  chemical  OTV/Shuttle 
rendezvous  strategy  of  Model  III  with  an  ion  propelled  OTV 
for  satellite  deployment.  The  use  of  an  ion  propelled  OTV 
drives  the  number  of  Shuttle  launches  to  a  value  comparing 
favorably  with  direct  Shuttle  insertion  of  satellites. 
Indeed,  if  longer  deployment  times  are  acceptable  Model  V 
yields  annual  Shuttle  traffic  levels  (7.80  launches/yr)  that 
are  significantly  less  than  the  11  launches/yr  estimated 
earlier.  Model  VI  accomplishes  a  trade  off  analysis  between 
satellite  deployment  times  and  the  annual  number  of  Shuttle 
launches  in  Model  V.  Figure  4.1  clearly  shows  the  trade 
off  between  deployment  times  versus  annual  Shuttle  traffic. 
By  properly  choosing  the  SOC  inclination,  altitude  and  OTV 
size,  average  deployment  times  can  be  dramatically  reduced 
without  substantially  increasing  Shuttle  traffic. 
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Figures  4.2,  4.3  and  4.4  illustrate  the  same  trade  off 
analysis  with  lower  more  realistic  estimates  of  satellite 
traffic.  In  addition,  several  ion  engine  and  Shuttle  flight 
parameters  have  been  varied  to  eliminate  some  of  the  conser¬ 
vatism  from  the  model.  As  previously  estimated,  direct 
orbital  insertion  of  the  Low  Traffic  Model  would  require 
about  eight  annual  Shuttle  launches.  The  corresponding 
efficient  operating  frontier  for  an  Extended  STS  is  illus¬ 
trated  in  Fig.  4.2  for  three  technology  levels.  As  before, 
annual  Shuttle  launches  are  less  than  eight  with  reasonable 
deployment  times.  The  technology  projections  in  the  figure 
drastically  reduce  deployment  times  and  slightly  reduce 
annual  Shuttle  traffic.  Figures  4.3  and  4.4  illustrate  the 
efficient  operating  frontiers  for  two  SOCs,  one  deploying 
missions  one  to  three  and  the  other  missions  four  to  seven. 
In  this  case  the  combined  number  of  launches  is  slightly 
greater  (still  less  than  eight),  but  the  deployment  times 
are  again  drastically  reduced  even  for  current  technology. 
The  reduced  deployment  times  possible  with  improved 
technology  and/or  the  deployment  of  two  SOCs,  combined  with 
an  on-orbit  spares  philosophy,  can  satisfy  the  military 
requirement  of  rapid  deployment. 

The  efficient  operating  frontier  for  an  Extended  STS 
supplied  by  WTR  Shuttle  launches  is  included  in  Fig.  4.2  and 
4.4.  Both  show  that  ETR  launches  supplying  a  one  or  two  SOC 
scenario  are  marginally  more  economical  than  WTR  launches. 
However,  this  data  is  somewhat  suspect  due  to  the  doubtful 
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quality  of  the  WTR  MECO  mass  calculations  between  56°  and 


70^,  and  the  lack  of  a  realistic  Traffic  Model. 
Nonetheless,  if  future  more  rigorous  analysis  substantiates 
these  results  the  deployment  of  an  Extended  STS  may 
eliminate  the  need  for  a  WTR  Shuttle  launch  site. 

Although  the  Extended  STS  launch  operations  costs  are 
less  than  those  required  to  directly  insert  satellites  into 
orbit  with  the  Shuttle,  other  cost  considerations  also 
exist.  Both  development  and  orbital  operations  costs  will 
contribute  to  the  life  cycle  costs  of  the  Extended  STS.  The 
latter  can  probably  be  accurately  estimated  while  the  former 
may  not  be  as  easy  to  quantify.  However,  by  building  a 
simple  logistics  depot,  development  costs  can  probably  be 
contained  at  a  reasonable  level  without  any  impact  on  the 
stations  operational  utility. 

Again,  although  life  cycle  costs  will  be  pivotal  in  any 
decision  to  deploy  an  Extended  STS,  there  are  additional 
considerations.  The  Extended  STS  offers  a  variety  of 
inherent  advantages  (and  a  few  disadvantages)  over  current 
STS  operations.  Among  the  advantages  are: 

1.  The  Shuttle  can  be  loaded  to  100%  capacity  for 
every  launch.  Current  operations  are  well  below 
this  figure. 

2.  Utilizing  the  Shuttle  as  a  simple  cargo  vessel 
will  reduce  the  complexity  of  payload  integration 
requirements . 

3.  The  current  upper  stages  required  to  deploy 
high  altitude  satellites  will  no  longer  be  needed. 


4.  More  extensive  orbital  test  and  validation 
operations  will  enhance  satellite  reliability  and 
reduce  early  on-orbit  failure  rates.  Ultimately, 
Factory-to-Orbit  testing  may  significantly  reduce 
the  cost  of  system  level  testing  at  the  factory 
and  launch  base. 

5.  Cost  effective  retrieval,  repair  and  servicing 
of  operational  satellites  can  extend  satellite 
lifetimes. 

6.  Orbital  satellite  spares  can  be  stored  at  the 
station  or  at  a  slightly  different  altitude. 
Satellites  deployed  at  different  altitudes  are 
periodically  accessible  with  a  small^V,  yet  are 
geographically  dispersed  and  thus  less  vulnerable 
to  attack  (see  section  3.4). 

7.  An  ABM  or  ASAT  system  can  be  similarly  deployed 
and  periodically  serviced  by  the  OTV. 

8.  The  Extended  STS  reduces  the  annual  number  of 
Shuttle  launches  required  by  current  STS 
operations  and  therefore  significantly  reduces 
costs . 

9.  The  Extended  STS  will  give  the  United  States  a 
manned  presence  in  space  and  a  platform  for  the 
future  expansion  of  manned  activities. 


Disadvantages  of  the  Extended  STS  include  the 
development  and  orbital  operations  costs,  increased  opera¬ 
tional  complexity,  the  risk  of  transporting  LOX  and  Hydrogen 
fuel  into  orbit,  longer  deployment  times  for  the  ion 
propelled  OTV,  and  the  station’s  lack  of  survivability.  The 
latter  consideration  is  often  overstated  since  the  Extended 
STS  is  at  least  as  survivable  as  the  launch  site  and 
considerably  less  prone  to  sabotage.  Further,  destruction 
of  a  SOC  utilized  as  a  logistics  depot  would  not  affect  the 
short  term  viability  of  operational  forces.  Other  advan¬ 
tages  and  disadvantages  can  probably  be  readily  identified, 
but  the  above  summarize  the  key  considerations. 


6.0  Conclusion 


The  current  STS  augmented  by  a  SOC  and  OTV  appears  to 
offer  many  attractive  advantages  over  current  Shuttle 
operations.  Such  an  Extended  STS  will  allow  the  Shuttle  to 
be  utilized  as  a  cargo  vessel  rather  than  as  an  orbital  work 
platform.  The  Shuttle  launched  to  the  SOC  can  consistently 
be  loaded  to  100S  of  its  capacity.  Directly  inserting 
payloads  into  orbit  with  the  Shuttle  potentially  requires 
more  launches  than  transporting  all  satellite  traffic  to  the 
SOC  and  deploying  each  vehicle  to  its  final  orbit  with  an 
OTV.  The  resultant  reduction  of  launch  costs  is  one  of 
three  major  elements  in  the  Extended  STS  life  cycle  costs: 
development  and  orbital  operations  being  the  other  two. 
Only  the  launch  costs  of  the  Extended  STS  were  estimated  in 
this  study. 

The  use  of  both  chemical  and  ion  propelled  OTVs  were 
evaluated  for  use  in  the  Extended  STS.  The  chemically 
propelled  OTV  with  a  specific  impulse  of  455  seconds  was 
uneconomical  except  for  the  limited  mission  of  deploying 
satellites  to  inclinations  below  28.5°.  Without  suffering  a 
massive  payload  penalty  the  Shuttle  is  incapable  of 
launching  to  inclinations  below  the  28.5°  ETR  latitude. 
Consequently,  the  extended  STS  has  a  natural  advantage  over 
direct  Shuttle  insertion  of  low  inclination  satellites.  The 
results  of  Model  II  (Section  4.2)  show  that  over  a  16  year 


period  the  Extended  STS  requires  5.70  annual  Shuttle 
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launches  to  deploy  a  SOC,  OTV  and  the  three  satellite 
missions  in  the  Traffic  Model  with  inclinations  below  28.5°. 
Model  III  (Section  4.3)  offers  additional  savings  by  rendez¬ 
vousing  the  OTV  with  the  Shuttle  in  low  Earth  orbit  and 
returning  the  payload  to  the  SOC.  The  calculated  annual 
Shuttle  traffic  for  that  scenario  was  4.51,  a  21t  savings. 
Both  of  these  figures  can  be  compared  with  a  crude  estimate 
of  about  3.4  annual  Shuttle  launches  (Section  5.0)  for 
direct  insertion  of  Shuttle  payloads.  Although  the  Extended 
STS  with  a  chemically  propelled  OTV  is  marginally 
uneconomical  for  low  inclination  missions,  the  economy  of 
scales  possible  by  including  civilian  and  foreign  satellite 
traffic  will  further  reduce  costs.  Recalling  that  the 
Extended  STS  provides  a  permanent  manned  presence  in  space 
and  obviates  the  need  for  many  expensive  upper  stages. 
Shuttle  launch  rates  actually  compare  favorably  with  direct 
insertion  of  Shuttle  payloads  even  with  the  Satellite 
Traffic  Model  used  in  this  study. 

The  OTV  specific  impulse  is  the  most  sensitive 
parameter  in  the  evaluated  Shuttle  launch  minimization 
models.  A  50%  increase  in  the  Model  III  OTV  specific 
impulse  results  in  a  49%  decrease  in  annual  Shuttle  launches 
with  one  SOC  servicing  all  seven  satellite  missions  in  the 
Traffic  Model  (Table  4.9).  Consequently,  the  use  of  ion 
engines  on  the  OTV  with  their  high  specific  impulses  vastly 
reduces  annual  Shuttle  traffi  •'n  t*-  Extended  STS.  Use  of 
existing  30  cm  diameter  mercury  ion  thrusters  provides 
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reasonable  satellite  deployment  times  which,  when  combined 
with  the  concept  of  orbital  spares,  can  satisfy  the  military 
requirement  of  rapid  deployment.  Figure  4.2  shows  the  trade 
off  between  the  average  OTV  mission  duration  and  annual 
Shuttle  launches  for  one  SOC  servicing  the  entire  Satellite 
Traffic  Model.  Using  the  Low  Traffic  Model  and  deploying 
the  satellites  via  direct  Shuttle  insertion  without  an 
Extended  STS  requires  approximately  8  annual  Shuttle 
launches.  Using  an  Extended  STS  with  30  cm  mercury  ion 
engines  mounted  on  the  OTV  requires  8.16  annual  Shuttle 
launches  per  year  with  an  average  OTV  mission  duration  of 
260  days.  Deploying  two  SOCs  at  28.5°  and  57.0°  inclina¬ 
tions  reduces  annual  Shuttle  traffic  to  7.3  launches  with  a 
170  day  average  OTV  mission  duration.  Alternatively, 
projecting  the  ion  engine  and  Shuttle  technology  into  the 
1990  timeframe  for  a  single  SOC  scenario  reduces  Shuttle 
traffic  to  6.31  annual  launches  with  a  120  day  average  OTV 
mission  duration  (Table  4.15).  In  all  cases  the  annual 
Shuttle  traffic  can  be  further  reduced  if  longer  OTV  mission 
durations  are  deemed  acceptable. 

Combining  a  dual  SOC  scenario  with  the  uprated  ion 
engine  and  Shuttle  technology  available  by  1990  should 
further  reduce  both  Shuttle  traffic  and  the  average  OTV 
mission  duration.  Annual  Shuttle  traffic  should  be  below 
five  to  six  launches  and  the  average  OTV  mission  duration 
well  below  100  days.  Consequently,  using  the  figure  $56 
million  per  Shuttle  launch  estimated  in  a  recent  General 
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Accounting  Office  report,  an  Extended  STS  may  save  as  much 
as  several  hundred  million  dollars  in  annual  Shuttle  launch 


costs.  Accomodating  future  missions  is  simulated  by 
increasing  the  SOC  and  annual  cargo  mass.  Annual  Shuttle 
traffic  is  relatively  insensitive  to  variations  in  both 
parameters  and  only  moderately  increases  launch  costs. 

The  Extended  STS  composed  of  either  a  chemically  or  ion 
propelled  OTV  obviates  the  need  for  the  current  fleet  of 
expensive  upper  stages,  and  provides  an  orbital  platform  for 
the  expansion  of  manned  operations  in  space.  The  launch  and 
orbital  operating  costs  of  the  system  compare  favorably  and 
may  be  less  than  the  cost  of  current  Shuttle  operations. 
Perhaps  more  important  from  a  military  perspective,  the  SOC 
can  be  utilized  as  a  logistics  depot  from  which  the  current 
satellite  fleet  can  be  serviced.  Orbital  test,  repair  and 
service  operations  at  the  SOC  can  enhance  the  reliability 
and  increase  the  lifetime  of  satellites  in  orbit.  The 
evolutionary  development  of  a  Factory-To-Orbit  test  concept 
may  eventually  reduce  expensive  factory  and  launch  base 
testing.  Combined  with  new  satellite  designs  entailing 
modular  maintenance,  system  self  testing  and  simplified 
procedural  testing  military  personnel  in  orbit  may  one  day 
be  able  to  routinely  salvage  billions  of  dollars  in  satellite 
hardware . 

The  models  developed  in  this  study  have  several 
applications.  For  the  strategic  planner  they  offer  a  means 
of  quantifying  the  cost  of  various  Extended  STS 
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configurations.  With  a  few  refinements  and  a  more  rigorous 
analysis  they  could  be  very  useful  in  the  current  national 
debate  over  whether  to  deploy  a  SOC.  For  the  engineer  they 
offer  considerable  insight  into  how  to  efficiently  design  an 
Extended  STS.  They  are  especially  useful  to  the  ion  engine 
designer  trying  to  build  an  efficient  and  effective  OTV.  As 
an  engineering  tool  the  SUMT  program  turned  out  to  be  an 
efficient  numerical  optimization  routine  for  this  type  of 
problem.  Finally,  as  an  academic  exercise  they  have 
radically  altered  my  own  views  on  how  to  best  exploit  the 
military  and  civilian  potential  of  the  space  environment. 

7.0  Research  and  National  Space  Policy  Recommendations 

Although  this  thesis  provides  a  framework  for 
evaluating  the  Extended  STS,  much  more  research  is  needed. 
The  Satellite  Traffic  Model  needs  to  be  refined,  and  more 
realistic  military  and  civilian  satellite  traffic 
incorporated.  In  addition,  the  scenarios  and  assumptions  of 
Sections  3.0  and  4.0  need  to  be  fine  tuned  by  knowledgable 
personnel  involved  in  STS  operations  and  future  planning. 
Finally,  a  means  of  quantifying  the  development  and  orbital 
operations  costs  needs  to  be  established  and  linked  to  the 
overall  life  cycle  costs  of  the  Extended  STS.  Included 
within  this  evaluation  would  be  some  means  of  quantifying 
the  specified  advantages  of  an  Extended  STS,  especially  the 
savings  due  to  satellite  retrieval,  repair  and  service. 
Additionally,  the  impact  of  accomplishing  future  missions 
such  as  deploying  and  maintaining  an  ABM  or  ASAT  system 
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needs  to  be  incorporated  into  the  models,  possibly  utilizing 
some  of  the  ideas  of  Section  3.6.2. 

With  regards  to  a  national  space  policy  the  results  of 
this  research  tentatively  suggests  that  the  military  and 
civilian  utility  of  an  Extended  STS  is  best  pursued  by  an 
evolutionary  approach.  Four  operational  phases  that  fall 
out  of  this  study  are: 


1.  Deploy  a  free-flying  reusable  chemical  OTV  to 
a  28.5°  inclination.  The  vehicle  would  be  used 
to  deploy  low  inclination  and  especially 
Geosynchronous  satellites.  Deployment  of  such 
satellites  eliminates  the  need  for  mission 
peculiar  upper  stages  and  will  prove  the  basic 
concept  of  deploying  satellites  via  a  reusable 
OTV. 


2.  Add  a  low  inclination  SOC  to  the  constellation 
as  soon  as  possible.  The  SOC  will  prove  the 
concept  of  doing  logistics  depot  work  in  orbit. 

3.  Deploy  an  ion  propelled  OTV  to  reduce  fuel 
costs  and  increase  the  number  of  satellite 
missions  serviced.  In  addition,  modify  the 
Chemical  OTV  for  manned  missions  and  retrieval  of 
low  Earth  orbit  Shuttle  payloads. 

4.  Incorporate  satellite  retrieval,  service, 
repair  and  refurbishment  operations.  Consider 
deploying  an  additional  SOC  at  higher  inclinations 
to  reduce  ion  propelled  OTV  deployment  times. 
Other  operational  missions  such  as  Satellite 
Storage  or  an  ABM/ASAT  system  can  also  be 
implemented  during  this  phase. 


Although  these  recommendations  are  somewhat  subjective, 
based  on  this  research  they  appear  to  offer  a  reasonable 
balance  between  designing  a  cost  effective  and  an 
operationally  effective  Extended  STS. 
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Appendix  A 

General  STS  Information 


Modeling  the  Shuttle  payload  mass  as  a  function  of 
inclination  and  altitude  was  difficult  due  to  a  paucity  of 
information  about  Shuttle  flight  characteristics.  Acquisi¬ 
tion  of  the  NASA  document  '’Ascent  Performance  and  Payload 
Estimation  Technique  for  Nominally  Shaped  Operational 
Missions"  (Ref  5)  was  an  invaluable  aid.  The  document 
allowed  Shuttle  MECO  masses  to  be  calculated  for  particular 
Shuttle  configurations.  Tables  3*1  and  3.2  delineate  the 
mass  configuration  and  assumptions  made  in  calculating  the 
Shuttle  MECO  mass.  These  masses  were  tabulated  and  modeled 
as  a  function  of  inclination  with  a  third  order  regression 
analysis.  The  resulting  equations  (Eqs.  3.^)  are  depicted 
for  both  the  ETR  and  WTR. 

As  mentioned  in  section  3.1  the  error  associated  with 
these  MECO  mass  calculations  is  small.  Computer  runs  A-1 
and  A-2  list  all  of  the  "actual"  calculated  MECO  masses  as 
well  as  the  corresponding  "predicted"  masses  from  Eqs.  3.4. 
Reference  5  estimates  that  the  actual  masses  are  accurate  to 
within  136  Kg  (300  lbs)  for  ETR  and  227  Kg  (500  lbs)  for  WTR 
launches.  The  correlation  coefficients  for  Eqs  3.4  are 
greater  than  .99999  for  both  the  ETR  and  WTR  functions, 
while  the  standard  deviations  are  4.1  and  6.1  Kg 
respectively.  Consequently,  if  the  NASA  error  estimates  are 


accurate,  Eqs.  3.4  should  be  able  to  predict  MECO  mass  to 
within  the  maximum  3-sigma  error  of  148  Kg  for  ETR  and  245 
Kg  for  WTR  launches. 

The  MECO  mass  calculation  presumes  that  the  shuttle  is 
launched  to  a  set  of  standard  MECO  conditions.  For  the  ETR 
the  conditions  are  an  altitude  of  57  NM,  a  flight  path  angle 
of  0.65  degrees  and  an  inertial  velocity  of  25680  fps.  The 
WTR  conditions  are  identical  except  the  inertial  velocity  is 
25374  fps.  Although  these  conditions  could  be  varied  for  a 
particular  mission,  the  majority  of  planned  launches  will  be 
targeted  to  achieve  standard  conditions  at  MECO. 

With  the  MECO  mass  known,  the  only  other  variables 
needed  to  calculate  payload  mass  are  the  Shuttle  reference 
mass  and  the  OMS  fuel  consumption.  The  former  depends  on 
the  Orbiter  mass  configuration  detailed  in  Table  3.1  while 
the  latter  can  be  calculated  from  the  velocity  change 
required  to  achieve  a  given  altitude.  This  velocity  change 
is  estimated  by  NASA  mission  planners  (Ref  6)  using  the 
equations : 


3.55  HOT  -  137  fps  for  HOT  <  175  NM 

AVAetr  =  '  3.44  HOT  -  117  fps  for  HOT  >  175  NM 

y  _  281  fps 

(Eqs.  A-1) 

'  276  fps  for  HOT  <  130  NM 

.72  HOT  +  183  fps  for  130  <  HOT  <  170  NM 

aVD  =  <  1.3  HOT  +  84  fps  for  170  <  HGT  <  230 

1.46  HGT  +  48  fps  for  HGT  >  230 
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where  AVAetr»  ^VAv/tr  represent  the  velocities 
required  to  ascend  to  a  given  altitude  (HGT)  and  then  to  de¬ 
orbit.  The  equations  include  a  velocity  reserve  which 
equates  to  an  QMS  propellant  reserve. 

The  altitude  of  Eqs.  A-1  can  easily  be  converted  to 
distance  from  the  Earths  center  (radius).  A  second  order 
linear  regression  can  then  be  run  on  the  radius  and  its 
corresponding  velocity  change.  Computer  runs  A-3  and  A-4 
list  the  "actual”  and  "predicted"  radii  and  aVs  (in  MKS 
units).  The  corresponding  equations  for  AVAetr, 
AVAwtr  listed  in  section  3.1  (Eqs.  3.5)  as  a 
function  of  radius.  The  correlation  coefficients  for 
A  VAetr  0.99999  and  0.99636,  while  the  standard 
deviations  are  0.23649  and  3. 61089  respectively. 

With  the  MECO  mass,  aVA  and  aVD  known,  Eq.  3.8  can  be 
applied  to  calculate  the  payload  mass.  Due  to  the  accuracy 
of  the  initial  MECO  mass  calculation  the  derived  payload 
masses  should  also  be  very  accurate.  Indeed,  because  the 
estimates  for  Shuttle  reference  mass,  ascent  and  deorbit 
AVs,  and  the  RCS  fuel  consumption  are  somewhat 
conservative,  it  is  very  likely  that  the  calculated  payload 
mass  is  somewhat  conservative.  Better  estimates  of  the 
above  variables  by  more  knowledgable  individuals  would  allow 
for  a  more  accurate  determination  of  payload  mass  as  a 
function  of  inclination  and  altitude. 


Computar  Run  A-1  ETR  MECQ  Mass  Versus  Inclination 

Regression  Analysis 


Variable  Name  and  Type: 


Inclination  =  X  ~>  Independent  variable 
Mass  =  Y  =>  Dependent  variable 

Third  Order  Regression  Analysis  Coe-f -f i ci ants: 


Zero  degree  coe-f-f icient  =  +159275.829 
First  degree  coe-f-f i ci ent  =  +  3.3575 

Second  degree  coe-f-f  icient  =  -  3.70994 

Third  degree  coe-ff icient  «  -*•  0.01339 

Regression  Analysis  Parameters: 


Variance  o-f  estimate  =  16.84896 
Standard  error  o-f  estimate  =  4.10475 
Correlation  coe-f-f  icient  =  0.99999 
Degrees  o-f  -freedom  =  12 

Table  o-f  Residuals: 


Actual  ’X' 

Actual  "Y' 

Predicted  ’ 

Y’ 

Residual 

4 

X 

28.4  deg 

156688  kg 

156685.561 

kg 

2.4388 

2 

30 

156399 

156399.083 

-0.0325 

32 

156020 

156022.987 

-2.9366 

4 

34 

155625 

155627.493 

-2. 493 

5 

36 

155212 

155213. 245 

-1.244.6 

6 

38 

154782 

154780.884 

1. 1163 

7 

40 

154335 

154331.053 

3.9469 

8 

42 

153367 

153864.395 

2.6046 

9 

44 

153382 

153381 . 553 

U. 4467 

10 

46 

152882 

152883. 169 

-1 . 1693 

11 

48 

152368 

152369.886 

-1.836 

12 

50 

151840 

151842.346 

-2.3461 

13 

52 

151301 

151301. 192 

-0. 1923 

14 

54 

1 50748 

150747.067 

0 . 9328 

15 

5o 

130182 

150180.613 

1 . 3368 

16 

58 

149602 

149602. 473 

-0. 4731 
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Computer  Run  A-2  WTR  MECO  Mass  Versus  Inclination 

Regression  Analysis 


Variable  Name  and  Type: 


Inclination  =  X  =>  Independent  variable 
Mass  =  Y  =>  Dependent  variable 

Third  Order  Regression  Analysis  Coe-f  +  icients: 


Zero  degree  coe-f-f ici ent  =  +164628.209 
First  degree  coe-f-f i cient  =  -  79.17494 

Second  degree  coe-f-f  ici  ent  =  -  3.01734 

Third  degree  coe-f-f  ici  enc  =  +  0.01256 

Regression  Analysis  Parameters: 


Variance  o-f  estimate  =  36.63971 
Standard  error  of  estimate  =  6.05307 
Correlation  coefficient  =  0.99999 
Degrees  of  freedom  =  17 

Table  of  Residuals: 


Actual  ’X' 

Actual  ’Y' 

Predictec  •’ 

Y" 

S  1  'Ij  LI  ^  1 

1 

70  deg 

148620  kg 

148609. 145 

kg 

1 0 . 3553 

o 

/  ^ 

147974 

147973. 791 

0.2087 

74 

147330 

147336. 003 

-  6.0032 

4 

76 

146683 

146696. 383 

-  8.3833 

5 

78 

146047 

146055.534 

-  8.5345 

6 

30 

145409 

145414.060 

-  5.0597 

7 

82 

144771 

144772.562 

-  1.5613 

8 

84 

144136 

144131.644 

4. 35o3 

9 

86 

143503 

143491.908 

11.0917 

10 

88 

142862 

142353.958 

8.0417 

11 

90 

142224 

142218.397 

5.6031 

12 

92 

141588 

141585.327 

2. 1733 

13 

94 

1 40957 

140956.851 

0. 1491 

14 

96 

140330 

140332.072 

-  2.0721 

15 

98 

139708 

139712.093 

-  4.0933 

16 

100 

139093 

139097.517 

-  4.5:74 

17 

102 

138484 

138488.947 

-  4.9473 

18 

104 

137383 

137886.986 

-  3.9860 

19 

106 

137292 

137292.236 

-  0.2361 

20 

108 

136708 

136705. 301 

2.6993 
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Computer  Run  A-4  ETR  Shuttle  Orbit  Radius  Versus 

Deorbit  Character! st i c  Velocity 
Regression  Analysis 


Variable  Name  and  Typei 


Shuttle  Orbit  Radius  =  X  =>  Independent  variable 
Characteristic  velocity  =  Y  =>  Dependent  variable 

Second  Order  Regression  Analysis  Coe-f f i ci ents: 


Zero  degree  coe-f-f icient  =  +7994.436744 
Rirst  degree  coe-f-f  icient  =  -  2.50173195 

Second  degree  coe-f-f  icient  =  +  0.000197^610171 

Regression  Analysis  Parameters; 


Variance  o-f  estimate  =  13.03351 
Standard  error  o-f  estimate  *  3.61089 
Correlation  coe-ff icient  »  0.99636 
Degrees  o-f  +reedom  =>41 

Table  o-f  Residuals: 


Actual  ’X’ 

Actual  'Y 

64PS  km 

34. 1  m/s 

6503 

84.  1 

3 

6513 

34.  1 

4 

6528 

34.  1 

5 

6538 

84.  1 

6 

6548 

84.  1 

7 

6558 

34.  1 

8 

6568 

84.  1 

9 

6578 

84.  1 

10 

6588 

84.  1 

11 

6593 

84.  1 

12 

6603 

84.  1 

13 

6618 

84.  1 

14 

6628 

85.4 

15 

6638 

86.  6 

16 

6648 

87.8 

17 

6658 

89.0 

18 

6668 

90.  1 

19 

6678 

91.3 

20 

6698 

94.  1 

Predicted  'V'  RssiDual 


73. 1035  m/ s 

5.9915 

78. 7766 

5. 3234 

79.4841 

4.6159 

80.2311 

3 . 3689 

81.0177 

m  O  ^  W 

81.8437 

2.2563 

82.7092 

1 . 3908 

83.6142 

0.4858 

84.5587 

-0. 4597 

85.5427 

-1.4427 

86.5661 

-2.4661 

37.6291 

-3.5291 

83.7316 

-4.6316 

89.8735 

-4. 4735 

9 1 . 0550 

-4.4550 

92.2759 

-4. 4759 

9t».  5w6o 

-4.5363 

94.8363 

—4 . 7363 

96. 1757 

-4.8757 

98.9730 

-4.8730 
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Table  o-f  Residuals  (Continued): 


Actual  ’X' 

Actual  ’Y’ 

Predicted  'V’ 

Resi dual 

21 

6718 

93.3 

101.9232 

-3.6282 

22 

6738 

102.6 

105.0415 

-2.4415 

23 

6758 

106.9 

108.3127 

-1.4127 

24 

6778 

111.2 

111.7419 

-0.5419 

25 

6798 

115.5 

115.3290 

0. 1710 

26 

6813 

120.  4 

119.0741 

1 . 3259 

27 

6838 

125.2 

122.9772 

O  00^*0 

^  m 

23 

6858 

130.0 

127.0383 

2.9617 

29 

6373 

134.8 

131.2573 

3.5427 

30 

6398 

139.6 

135.6343 

3.9657 

31 

6913 

144.4 

140. 1693 

4.2307 

32 

6938 

149.2 

144.3622 

4.3378 

33 

6958 

154.0 

149.7131 

4.2869 

34 

6978 

153.8 

154.7219 

4.0781 

35 

6998 

163.6 

159.3888 

3.7112 

36 

7018 

163.4 

165.2136 

3. 1864 

37 

7038 

173.2 

170.6964 

2.5036 

33 

7058 

173.0 

176.3371 

1 . 6629 

39 

7078 

132.8 

182. 1358 

0.6642 

40 

7098 

187.6 

138.0925 

-0.4925 

41 

7118 

192.4 

194.2071 

-1.3071 

42 

7138 

197.3 

200. 4798 

-3. 1793 

43 

7158 

202.  1 

206.9103 

-4.8103 

7178 

206.9 

213.4939 

-6.5939 
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Scale  Height  Modeled  with  Skylab  Orbit  Maintenance  Data 

The  atmospheric  model  of  section  3.2  is  used  to 
calculate  the  orbit  maintenance  fuel  cost  of  maintaining  the 
SOC  at  constant  altitude.  Fuel  cost  was  calculated  from  the 
average  SOC  acceleration  due  to  air  drag.  Acceleration, 
calculated  by  Eq.  3.11,  is  in  turn  a  function  of  the  SOC 
ballistic  coefficient,  inclination,  radius  and  scale  height. 
All  of  these  parameters  are  known  or  can  be  accurately 
estimated . 

Initially  the  scale  height  was  modeled  by  linear 
regression.  The  mean  COSPAR  international  Reference  Atmos¬ 
phere  (CIRA)  with  the  values  of  scale  height  averaged  over 
the  day/night  cycle  was  used  to  derive  scale  height  as  a 
function  of  radius.  Unfortunately,  the  curve  fit  was  not 
very  good,  and  when  used  in  the  SUMT  minimization  routine  it 
gave  unrealistically  high  estimates  of  fuel  consumption. 
Early  Skylab  data  was  then  used  to  estimate  a  realistic 
value  for  the  scale  height.  To  resolve  the  dilemma  a 
constant  value  of  scale  height  was  decided  on  since  it  could 
be  varied  to  give  realistic  estimates  of  fuel  consumption. 

Skylab  was  a  large  space  station  placed  at  an 
inclination  of  50°  and  a  radius  of  6811  km  (23^NM  altitude). 
Data  from  the  mission  w^s  used  to  calculate  a  realistic 
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value  for  the  scale  height.  Figure  B-1  depicts  the  Skylab 
A  V  maneuvers  required  for  orbit  maintenance  through  the 
first  two  manned  flights.  The  five  av  maneuvers 
accomplished  during  the  144  days  of  flight  correspond  to  an 
average  acceleration  due  to  air  drag  of  1.8  x  10“7  m/sec^. 
Alternatively,  the  three  final  maneuvers  accomplished  during 
the  last  64  days  correspond  to  an  acceleration  of  3.0  x  10-7 
m/sec2.  In  addition  to  estimating  the  Skylab  acceleration 
due  to  air  drag,  mass  and  surface  area  figures  were  used  to 
calculate  a  ballistic  coefficient  in  the  range  .015  to  .02 
(CjS/M) . 

The  above  figures  were  used  in  conjuction  with  Eq.  3.11 
to  calculate  the  scale  height  (H).  Rearranging  Eq.  3.11J 

^^o  -  ^soc^ 

ln[(2AsocxRsoc)/(UxPoxBx( .00175  Isoc  +  .84004))] 

where , 

Rq  =  6498  km 
Pq  =  24.9  kg/km^ 

U  =  398601 .2  km3/sec2 

Isoc  “ 

R 


r> 


=  6811  km 


^Agoc  /B)tnin  =  1 .2  x  10“5  kg/msec2 
(Agoc  /B)max  =  1.5  x  10“5  kg/msec^ 


The  resulting  values  of  scale  height  are  28.6  and  29.2  km. 
These  values  correspond  to  the  average  acceleration 
experienced  by  Skylab  during  its  first  five  months  in  orbit. 
The  specific  value  of  scale  height  will  vary  over  the 
day/night  cycle,  and  even  the  average  value  varies  with 
Solar  activity.  Nonetheless,  the  above  values  inserted  into 
the  atmospheric  model  of  section  3.2  provide  an  average 
estimate  of  the  Skylab  orbit  maintenance  fuel  consumption. 

Consequently,  a  scale  height  of  30  km  was  chosen  for 
use  within  the  orbit  maintenance  fuel  consumption  model. 
Applied  to  the  Skylab  data,  this  figure  of  scale  height 
results  in  a  slightly  conservative  estimate  of  fuel 
consumption.  Higher  values  of  scale  height  would  be  even 
more  conservative.  Although  not  perfect,  a  constant  scale 
height  is  adequate  for  the  purposes  of  this  study. 
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Computer  Programs 

The  computer  models  minimized  by  SUMT  are  contained  in 
this  appendix.  Becuase  of  the  similarity  of  the  six  models 
only  the  program  listings  for  Models  II  and  VI  are  included. 
For  the  sake  of  completeness  Model  II  includes  a  listing  of 
both  the  ETR  and  WTR  problems.  As  with  all  the  models  the 
ETR  and  WTR  problems  are  very  similar.  Wherever  differences 
exist  the  margins  are  marked  with  a  a  or  □  for  ETR  and  WTR 
problems  respectively. 

For  the  sake  of  understandability  the  programs  are 
internally  documented  and  contain  three  different 
subroutines  calculating  OTV,  Orbiter  and  SOC  fuel 
consumption  (subroutines  STSOPT,  ORBFUL  and  SOCFUL 
respectively).  A  fourth  subrouting  called  RESTNT  contains 
the  problem  objective  function  and  constraints.  The  program 
readability  is  enhanced  by  consistent  use  of  the  same 
variable  designations  as  contained  within  the  text  of  the 
thesis.  The  SUMT  library  program  is  documented  in  Ref  9  and 
10.  The  PROCES  program  uses  SUMT  as  a  nested  optimization 
routine  which  is  called  sequentially  and  is  documented  in 
Ref  14,15,21  and  22.  More  detailed  information  on  the 
programs  can  be  obtained  through  the  AFIT  Department  of 
Aeronautics  and  Astronautics,  Capt.  DeWispelare. 


As  with  most  Thesis  efforts  Murphy's  law  proved  itself 
infallible  once  again.  After  completing  more  than  one 
hundred  computer  runs  on  Models  I  through  IV  a  sign  error  in 
the  ETR  MECO  mass  calculation  was  found  (Marked  in  the 
margin  of  Computer  Program  C-1  with  o).  Fortunately,  the 
impact  was  very  minor  and  resulted  only  in  a  slightly 
conservative  estimate  of  MECO  mass  and  therefore  payload 
mass  for  Models  I  through  IV.  The  total  error  was  only  a 
couple  of  hundred  kilograms  or  a  fraction  of  a  Shuttle 
launch,  and  it  was  corrected  in  Models  V  and  VI.  The  errors 
are  appropriately  marked  in  the  attached  programs. 

The  vector  optimization  problem  of  Model  VI  would 
normally  be  accomplished  by  a  program  like  PROCES.  However, 
due  to  the  simplicity  of  a  dual  objective  problem  the 
solution  was  implemented  with  a  simple  variation  of  Model  V. 

A  loop  added  to  the  program  repeatedly  used  SUMT  to  minimize 
Shuttle  launches  subject  to  different  average  OTV  mission 
durations.  The  latter  were  incorporated  as  equality 
constraints  within  the  SUMT  program.  The  resulting 
solutions  were  then  manually  checked  to  get  the  NDSS  or 
efficient  frontier. 

Utilization  of  the  SUMT  program  as  a  numerical  minimi¬ 
zation  routine  was  particularly  efficient  in  the  single 
objective  problems,  models  I  through  V.  The  total  central 
processor  time  required  to  compile  and  execute  the  programs 
was  around  0.5  seconds.  The  vector  optimization  problem 
with  its  multiple  calls  to  SUMT  was  considerably  less  efficient. 


Computer  Program  C-1  ETR  Model  II  Launch 

Minimization  Program 


>=QTS. T35, lOlOO.CMlOCOOO. T820472, SP078 
>=ATTACH,SUMT,SUMT, ID=AFIT. 

)=LIBRARY,SUMT. 

:>=FTN(L=0)  . 

>=LBO. 

)=:(fEOR 

PROBRAM  MA I N  ( I NPUT ,  OUTPUT ,  SAVE ,  TAPE5=  INPUT ,  TAPE6=CUTPUT , 
>=  +  TAPE7=SAVE) 

)=  REAL  ASAT(35) ,ESAT<35) , ISAT<35) ,M3AT(35) . 

+  FSAT(35) ,MFMSN(35) ,PI,U, W,P0,3,NET, 

)=  +  MR,B,NOTV.MSOC,«CSOC,NFdTV.MFSDC,NFORB, 

:=  +  MPL, FORB, iSPOTV, ISPSDC, ISPORB 

>=  INTEBER  i.K 

>=  COMMON / DE VC /N I, NO. NS 

:>=  COMMON/VALUES/ASAT,  ESAT,  ISAT,  MSAT,  FSAT,  MFMSN, 

)=  +  P!,U,W,PO,S,MET.MR,K.B.MQTV.MSOC,MCSOC, 

)=  +  MFOT V , MFSOC , MFOFB , MPL , FORB , I SPOT V , I SP30C , 

>=  +  ISPORB,  H 

)=  CCMMON/SHARE7X (100) , DEL ( 100)  ,  A ( lOO. 100  ,  N  <5) 

:>=  NI=5 

)=  N0=6 

)=  NS=7 

)=C 

)=C  ENTERS  CONSTANT  PROGRAM  DATA 


DATA  PI,U,W,PO,S. MET , MR / 3 . 1 4 1 5926535 .398601.2, 

7 . 2922 1 1 5856E-5 ,24.9,3.7, 38399 . 0 , 8804 1 . 0/ 

ENTERS  PROBLEM  PARAMETER  DATA 

DATA  K . B . MOT V . MSOC , MCSCC / 7 , . 02 , 2270 . , 1 00000 . , 20000 . / 
DATA  H, ISPOTV, ISPSOC, I3P0RB/30. 0, 455. 0, 455. 0, 313. 0/ 

ENTERS  SATELLITE  TRAFFIC  MODEL  DATA 

DATA  <  ASAT  <I),I  =  1.7)/41 000 . . 20000 . , 6700 . , 65000 . , 

25000. .  12000. . 6700. / 

DATA  <ESAT<I) , 1=1 . 7) /O. ,0. ,0. ,0. ,0. 7,0. ,0. / 

DATA  <I3AT(I) , I=l,7)/0. ,0. ,28.5,55. ,65. ,90. ,98. / 

DATA  <MSAT ( I ) , 1=1 , 7) /2000. , 500. , 25000. , 1500. , 

1500. .  4500. . 8000. / 

DATA  f  FSAT  ( I )  .  1  =  1 , 7)  /3.  ,  .  5,  .  5,  4.  ,  .  5,  3.  ,  5.  / 

CALL  SUMT 

FORMAT (//,4X, "ETR  MODEL  II  LAUNCH  ", 

"MINIMIZATION  PROSRAM") 

FORMAT  < 4 X,  "===  =====  ==  ======  , 

PRINT  10 
PRINT  15 

FORMAT (//."  SATELLITE  TRAFFIC  MODEL") 
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600= 

610= 

30 

620= 

630= 

35 

640= 

650= 

36 

660= 

670= 

680= 

690= 

40 

700= 

710= 

720= 

730= 

740= 

50 

750= 

60 

760= 

70 

770= 

30 

780= 

81 

790= 

82 

800= 

83 

310= 

34 

820= 

85 

330= 

86 

340= 

37 

1^50= 

90 

“  860= 

95 

870= 

100 

880= 

101 

890= 

102 

P00= 

103 

910= 

1 10 

920= 

115 

930= 

120 

940= 

121 

950= 

122 

960= 

123 

970= 

980= 

990= 

1000= 

1010= 

1020= 

1030= 

1040= 

1050= 

1060= 

1070= 

1080= 

1090= 

PRINT  20 

FORMAT  (“ - ■',/) 

PRINT  30 

FORMAT (/,5X, "MISSION" , 3X , " ASAT" , 3X , "ESAT" , 3X , "ISAT", 
3X, "MSAT",4X,"FSAT",4X, "OTV  FUEL") 

F0RMAT(7X, " (#) ",5X, " <KM) ", lOX, " (DEG) ",2X, " (KB) ",4X, 
" (LS/YR) ",3X, " (KS) ",/) 

PRINT  35 
PRINT  36 

FORMAT (6X, 13, 5X , F6. 0, 2X , F4. 2, 3X . F5. 1, 1X,F7.0,2X, 
F5.2,2X,Fa.0) 

DO  50  1=1, K 

PRINT  ao, I , ASAT ( I ) , ESAT ( I ) , ISAT ( I ) , 

MSAT ( I ) , FSAT ( I ) , MFMSN ( I ) 

CONTINUE 

FORMAT ( / , "  VAR I ABLE  PARAMETERS " ) 

FORMAT (" - ",/) 

FORMAT ("  BALLISTIC  COEFFICIENT  (MS*2/KS) =" , FIO. 

FORMAT ("  SCALE  HEIGHT  (KM)  =",F10. 

FORMAT ("  OTV  MASS  (KG)  =",F10. 

FORMAT ("  SOC  MASS  (KG)  =",F10. 

FORMAT!"  ANNUAL  SOC  CARGO  MASS  (KG/YR)  =",F10. 

FORMAT!"  OTV  SPECIFIC  IMPULSE  (SEC)  =",F10. 

FORMAT!"  SOC  SPECIFIC  IMPULSE  (SEC)  =",F10. 

FORMAT!"  ORB  SPECIFIC  IMPULSE  (SEC)  =",F10. 

FORMAT!/,"  FUEL  CONSUMPTION") 

FORMAT ("  - ",/) 

FORMAT ("  ORBITER  FUEL  MASS  (KB/MSN)  =",F10. 

FORMAT!"  SOC  FUEL  MASS  (KG/YR)  =",F10. 

FORMAT!"  OTV  FUEL  MASS  (KG/YR)  =",F10. 

FORMAT!"  TOTAL  STS  FUEL  MASS  (KG/YR)  =",F10. 

FORMAT!/,"  SOC  LOCATION  AND  RESUPPLY  PARAMETERS") 

FORMAT!" - ",/) 

FORMAT!"  ORBIT  INCLINATION  (DE6>  =",F10. 

FORMAT!"  ORBIT  RADIUS  (KM)  =",Fu'!. 

FORMAT!"  NO  OF  ORB  LAUNCHES  (LS/YR)  =",F10. 

FORMAT!"  ORB  PAYLOAD  MASS  !K6/MSN)  =",F10. 

DO  130  1=1, 1 
PRINT  60 
PRINT  70 
PRINT  80. B 
PRINT  81, H 
PRINT  82, MOT V 
PRINT  83,MS0C 
PRINT  84,MCS0C 
PRINT  85,ISP0TV 
PRINT  86,ISPS0C 
PRINT  87. ISPORB 
PRINT  90 
PRINT  95 


(MS*2/KS) = 


!KG/YR)  = 


'.FIO. 

",F10. 

'.FIO. 

■,F10. 

•,F10. 

•,F10. 

'.FIO. 

',F10. 


=",F10. 

=",F10. 

=",F10. 

=",F10. 
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*100= 

PRINT  100,MF0RB 

1110= 

PRINT  lOljMFSOC 

1120= 

PRINT  102,MF0TV 

1130= 

PRINT  103, FORB*MFORB+MFSOC+MFOTV 

1140= 

PRINT  110 

1150= 

PRINT  115 

1160= 

PRINT  120, X(l) 

1170= 

PRINT  121, X(2) 

1180= 

PRINT  122,FCRB 

1190= 

PRINT  123, MPL 

1200=  130 

CONTINUE 

1210= 

END 

1220= 

SUBROUTINE  RE3TNT ( IN, VAL) 

1230= 

REAL  ASAT(35) .ESAT<35> , ISAT(35) .MSAT(35> , 

1240= 

FSAT(35) ,MFMSN(35) ,PI,U, W,PO,S,MET, 

1250= 

MR , B , MOTV , MSOC , MCSOC , MFOT V , MFSOC , MFORB , 

1260= 

+ 

MPL.FORB. ISPOTV, ISPSOC, ISPORB 

1270= 

COMMON /VALUES / ASAT, ESAT, ISAT, MSAT, FSAT, MFMSN, 

1280= 

+ 

PI ,  U,  W,  PQ,  S,  MET,  MR,  K,  B.  MOTV,  MSOC,  MCSOC, 

1290= 

+ 

MFOTV, MFSOC, MFORB, MPL. FDRB, ISPOTV, ISPSOC, 

1300= 

+ 

ISPORB. H 

1310= 

COMMON  /  SHARE /XdOO),  DEL  ( 1 00 ) ,  A  ( 1 00 ,  1 00 )  ,  N  ( 5 ) 

1320=C 

1330=C 

X(1)=S0C  INCLINATION  :  X(2)=S0C  ALTITUDE 

^40=C 

(p^;5o= 

IF<IN)  140.140.150 

1360=  140 

CALL  STSOPt<VAL) 

1370= 

RETURN 

13S0=C 

13'^0=C 

INCLINATION,  ALTITUDE  AND  FUEL  CONSTRAINTS 

1 4';!0=c 
1410=  150 

GGT0<151, 152, 153, 154, 155. 156> , IN 

1420=  151 

VAL=X(l)-28.5 

1430= 

RETURN 

1440=  152 

VAL=57.0-X<1) 

1450= 

RETURN 

1460=  153 

VAL=X (2) -6500.0 

1470= 

RETURN 

1480=  154 

VAL=7200.0-X (2; 

1490= 

RETURN 

1500=  155 

CALL  ORBFUL 

1510= 

VAL= 1 0830 . O-MFORB 

1520= 

RETURN 

1530=  156 

CALL  STSOPT<VAL) 

1540= 

VAL=MPL-0.0 

1550= 

RETURN 

1 560= 

END 

1570= 

SUBROUTINE  STSOPT<VAL> 

1580=C 

1 590=C 

CALCULATES  OTV  FUEL  CONSUMPTION  FOR  K  MISSIONS 

A 


A 


I 
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i'600=C 

1610= 

1620= 

1630= 

1640= 

1650= 

1660= 

1670= 

1680= 

1690= 

2  700= 
1710= 
1720= 
1730= 
1740= 
1750= 
1760= 
1770= 
1730= 
1790= 
1800= 
1810= 
1820= 
1830= 

240= 
'J50= 
1860= 
1870= 
1880= 
1890=  160 


1900= 

1910= 

1920=C 

1930=C 

1940=C 

1950= 

1960= 

1970= 

1980= 

1 990=C 
2000=C 
2010=C 


2020= 

2030= 

+ 

2040= 

+ 

2050= 

+ 

2060= 

+ 

2070= 

2080= 

+ 

-2iD90= 

+ 

REAL  I see . RSeC , MCSAT , VOTV , RABAT , M , P , R , S , T , 

ASAT(35) ,ESAT{35) , ISAT (35) . MSAT <35) , 

FSAT  (35)  ,  MFMSN  (35)  ,  PI ,  U,  W,  PO,  6,  MET, 

MR , B , MOT V , Msec , MCSOC , MFDT V , MFSOC , MFORB , 
MPL.FeRB, ISPOTV, ISPSOC, ISPeRB 
INTEGER  i 

CeMMON  /  SHARE  /  X  ( 1 00 )  ,  DEL  ( 1 00 )  ,  A  ( 1 00 ,  1 00 )  ,  N  ( 5 ) 

COMMON/ VALUES/ ASAT , ESAT, I SAT . MSAT , FSAT , MFMSN , 

PI , U, W, PO, G, MET. MR, K, B, MOTV , MSOC, MCSOC, 
MFOTV.MFSOC, MFORB, MPL.FDRB. ISPOTV, ISPSOC, 
ISPGRBjH 
MF0TV=0.0 
MCSAT=0.0 
ISCC=X(1) 

RSOC=X (2) 

DO  160  1=1, K 

RASAT=ASAT ( I ) * ( 1 . O+ESAT ( I ) > 

M=SQRT(ABS(2.0E6*U»(1.0/RS0C-1.0/ (RASAT+RSOC) ) ) ) 
P=SQRT ( U* 1 . 0E6/RS0C  > 

R=1.0E6*U» (4.0/RASAT-i.0/ASAT(I)-2.0/ (RASAT+RSOC) ) 
S=4. 0E6*U*SCRT (ABS ( ( ( 1 - O/RASAT-0. 5/ASAT ( I ) ) 

* (l.O/RASAT-1.0/ (RASAT+RSOC) ) ) ) ) 

T=COS( (ISOC-ISAT(I) )*PI/180.0) 

VOT V= ABS ( M-P ) +SQRT ( ABS ( R-S  *  T ) ) 

MFMSN ( I ) = (MSAT ( I ) +MOTV+MOTV»EXP ( VOTV/ ( ISPOTVCG) ) ) 

* (EXP (VOTV/ ( ISPOTV*G) ) -1 . 0) 

MFOTV=MFOTV+MFMSN ( I ) »FSAT ( I ) 

MCSAT=MCSAT+FSAT ( I ) »MSAT ( I ) 

CONTINUE 
CALL  SOCFUL 
CALL  0R3FUL 

CALCULATES  TOTAL  ANNUAL  FUEL  CONSUMPTION 

FORB= (MCSAT+MFOTV+MFSOC+MCSOC) /MPL 
VAL=FORB 
END 

SUBROUTINE  ORBFUL 

CALCULATES  ORB I TER  FUEL  CONSUMPTION 

REAL  ISOC, RSOC, MM, VORBA. VORBD, MFORBA. MFORBD, 

ASAT(35) ,ESAt(35) , ISAT (35) . MSAT (35) , 

FSAT(35) .MFMSN(35) ,PI.U,W.PO,B.MET, 

MR , B , MOTV , MSOC , MCSOC , MFOT V , MFSOC . MFORB , 

MPL.FORB, ISPOTV. ISPSOC, ISPORB 
COMMON/VALUES/ ASAT. ESAT. ISAT, MSAT, FSAT. MFMSN, 

PI .  U.  W,  PO.  Q,  MET.  MR,  K,  B,  MOTV.  MSOC,  MCSOC, 
MFOTV.MFSOC, MFORB, MPL.FORB, ISPOTV, ISPSOC, 
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^1 00=  + 

2110= 

2120= 

2130= 

2140= 

2.150=  + 

2160= 

2170=  + 

2180= 

2190=  + 

2200= 

2210=  + 

2220= 

2230= 

2240=  + 

2250= 

2260= 

2270= 

2280=C 

2290=C 

2300=C 

2310= 

2320=  + 

2330=  + 

^540=  + 

•350*  + 

2360= 

2370= 

2380=  + 

2390=  + 

240(;}=  + 

2410= 

2420= 

2430= 

2440= 

2450=  + 

2460= 

2470= 

2480= »E0R 
2490=  »DATA 
2500=*E0R 


I SPORE. H 

COMMCN/SHARE/XaOO)  ,DEL(100)  .AdOO,  100)  ,N(5) 
IS0C=X(1)  , 

RS0C=X(2.i  ^ 

MM= 1 59275 . 329-3 . 3575* I SCC-3 . 70994* I SOC*  *2 
-*•0.01339*1300**3 

V0RBA=-2. 023703587E-5*RS0C**2+0. 8474467398*RS0C 
-4623.658232 

VORBD= 1 . 9746 10171 E-4*RSCC*  *2-2 . 50 1 373 1 95*RS0C 
■•■7994.436774 

MPL= (MM-MET) *EXP (-VORBA/ ( ISP0RB*6) ) 
-MR*EXP(VORBD/<ISPORB*G) ) 

MFORBA=  <  MM-MET ) * ( 1 . 0-EXP  < -VORBA/ < 1 3P0RB*S ) ) ) 
MFGRBD= (MM-MET-MFORBA-MPL) 

* ( 1 . 0-EXP (-VORBD/ ( ISPORB*G) ) ) 
MFORB=MFCRBA-*-MFORBD 
END 

SUBROUTINE  SOCFUL 

CALCULATES  SOC  STATIONKEEPING  FUEL  CONSUMPTION 

REAL  ISOC.RSOC. ASOC.F. 

ASAT.'35)  ,ESAT<35)  ,  ISAT  (35)  ,  MSAT  (35)  . 
FSAT<35) .MFMSN(35) .PI.U.W.PO.GjMET, 

MR , B . MOT V , MSOC , MCSOC , MFOT V . MFSOC , MFORB , 
MPL,FORB, ISPOTV. ISPSOC, ISPORB 
COMMCN/SHARE/X  (100)  ,D£L<  100)  ,  AdOO,  100)  ,N  (5) 
COMMON/ VALUES/ ASAT , ESAT , I SAT , MSAT , FSAT , MFMSN , 

PI ,  U,  W,  PO,  G,  MET,  MR,  K,  B.  MOTV ,  MSOC,  MCSDC, 
MFOTO, MFSOC, MFORB, MPL,FORB, ISPOTV. ISPSOC, 
ISPORB, H 
ISOC=X  (1) 

RSCC=X (2) 

F=0 .  00 175*1  SOC-i-0 .  S4004 
ASOC=<U*B*PO*F 

*EXP( <6498.0-RS0C) /H) ) / (2000*RS0C) 

MFSOC=( MSOC* AS0C*31 557600) / ( ISPSOC*G) 

END  I 

y 

N=2,M=6, X=2S. ,6500. ,NT(5)=1.THETA0=1.E-12  *END 


Computer  Program  C-2  WTR  Model  II  Launch 

Minimization  Program 


100=QTS, T35, lOlOO, CMIOOOOO. TS20472, SP078 
1 1 0=ATTACH , SUMT , SUMT , I D=AF I T . 
120=LIBRARY,SUMT. 

130=FTN(L=0) . 


140=Le0. 

150=*E0R 

160= 

170=  + 

180= 

190=  + 

200=  + 

210=  + 

220= 

230= 

240= 

250=  + 

260=  + 

270=  + 

280= 

290= 

300= 

.-..-■310= 

fTri20=C 

^330=C 

340=C 

350= 

360=  + 

370=C 

330=C 

39C-C 

400= 


410= 

420=C 

430=C 

440=C 

450= 


460=  + 

470= 

480= 


490= 


500=  + 

510= 

520= 

530=  10 
540=  + 

550=  15 
560=  + 

--^70= 

■iS0» 

590=  20 


PROGRAM  MA I N (I NPUT , OUTPUT , SAVE , T APE5= I NPUT , T APE6=0UTPUT 
TAPE7=SAVE) 

REAL  ASAT(35> ,ESAT(35) , ISAT<35) ,MSAT(35> , 

FSAT(35) ,MFMSN<35> ,PI,U, W,PO,G,MET, 

MR , B , MOT V , MSOC . MCSOC , MFOT V , MFSOC , MFORB , 

MPL.FORB, I SPOT V, ISPSOC, ISPORB 
INTEGER  i,K 
COMMON/ DEVC/N I , NO, NS 

COMMON/ VALUES/ ASAT, ESAT, ISAT, MSAT, F3AT, MFMSN, 

PI ,  U,  W,  PO,  e,  MET,  MR,  K,  S,  MOTV,  MSOC,  MCSOC, 

MFOTV, MFSOC, MFORB, MPL.FORB, I SPOT V, ISPSOC, 

ISPORB, H 

COMMON/SHARE/XdOO)  ,DEL(100)  ,A(100,  100;  ,N(5) 

NI=5 

N0=6 

NS=7 

ENTERS  CONSTANT  PROGRAM  DATA 

data  PI,U,W,P0,G,MET,MR/3. 1415926535,398601.2, 

7. 29221 15856E-5, 24. 9, 8. 7, 38399.6, 8804i . 0/ 

ENTERS  PROBLEM  PARAMETER  DATA 

DATA  K , B , MOTV , MSOC , MCSOC / 7 , . 02 , 2270 . , 1 00000 . , 20000 . / 
DATA  H, ISPOTV, ISPSOC, ISPCRB/30. 0, 455. 0, 455. 0, 313. 0/ 

ENTERS  SATELLITE  TRAFFIC  MODEL  DATA 

DATA  (ASAT<I) , 1=1 . 7) /41000. , 20000. , 6700. , 65000. , 

25000. .  12000. . 6700. / 

DATA  <ESAT(I) , 1=1 , 7) /O. , 0. , 0. ,0. ,0.7,0. ,0. / 

DATA  <ISAT(I) , 1=1, 7) /O. ,0. ,28.5,55, ,65. ,90. ,98./ 

DATA  <MSAT ( I ) , 1=1 , 7) /2000. , 500. . 25000. , 1500. , 

1500. . 4500. .8000./ 

DATA  (FSAT  ( I)  ,  1  =  1 , 7)  /3.  ,  .  5,  .  5,  4.  ,  .  5,  3.  ,  5.  / 

CALL  SUMT 

FORMAT (//.4X, "WTR  MODEL  II  LAUNCH  ", 

"MINIMIZATION  PROGRAM") 

FORMAT (4X, "===  =====  ==  ======  % 

"ssssssssasas.'s  ==s==s=") 

PRINT  10 
PRINT  15 

FORMAT {//,"  SATELLITE  TRAFFIC  MODEL") 


3! 


[B 

o  •. 
i-'S 


600= 

PRINT  20 

610= 

30 

FORMAT ("  - 

",/) 

620= 

PRINT  30 

630= 

35 

FORMAT (/,5X, "MISSION" , 3X , "ASAT", 

3X,  "ESAT",3X, "ISAT", 

640= 

+ 

3X, "MSAT",4X, "FSAT", 4X, "OTV 

FUEL") 

650= 

36 

F0RMAT(7X, " (#) ",5X, " (KM) ", lOX, " 

(DEG) ",2X, " (KG) ",4X, 

660= 

+ 

(LS/YR)  ",3X.  "  (KG)  ",/) 

670= 

PRINT  35 

680= 

PRINT  36 

690= 

40 

FORMAT (6X, 13, 5X . F6. 0, 2X , F4. 2, 3X , 

F5. 1, 1X,F7.0,2X, 

700= 

710= 

720= 

730= 

740= 

750= 

760= 

770= 

7B0= 

790= 

800= 

810= 

820= 

830= 

|i|^40= 

^-850= 


50 

60 

70 

SO 

81 

82 

83 

84 

85 

86 
87 
90 


/  ; 


F5.2,2X,F8.0> 

DO  50  1=1, K 

PRINT  40. I,ASAT(1) ,ESAT<I) , ISAT(I) , 
MSAT (I ) , FSAT (I) , MFMSN ( I ) 
CONTINUE 

FORMAT</,“  VARIABLE  PARAMETERS") 
P'ORMATC' 

FORMAT  < " 

FORMAT ( " 

FORMAT  <  •• 

FORMAT ( " 

FORMAT ( " 

FORMAT  < " 

FORMAT  < " 

FORMAT  < " 

FORMAT </, 


BALLISTIC  COEFFICIENT 
SCALE  HEISHT  <KM) 

OTV  MASS  (K6) 

SOC  MASS  (KB) 

ANNUAL  SOC  CARGO  MASS 
OTV  SPECIFIC  IMPULSE 
SOC  SPECIFIC  IMPULSE 
ORB  SPECIFIC  IMPULSE 
FUEL  CONSUMPTION") 


(M**2/k;G)=",F10.2) 
=",F10.2) 
=",F10.2) 
=",F10.2) 
(KG/YR)  =",F10.2) 
(SEC)  =",F10.2) 

(SEC)  =",F10.2) 

(SEC)  =",F10.2) 


360= 

95 

FORMAT ("  - 

- - "./) 

370= 

100 

FORMAT ( " 

ORB I TER  FUEL  MASS  (KG /MSN) 

=",F10.2) 

830= 

101 

FORMAT ( " 

SOC  FUEL  MASS  (KB/YR) 

=",F10.2) 

890= 

102 

FORMAT ( " 

OTV  FUEL  MASS  (KG/YR) 

=",F10. 2) 

■ 

900= 

103 

FORMAT ( " 

TOTAL  STS  FUEL  MASS  'KG/VR) 

=  ‘‘,F10.2) 

pi 

910= 

110 

FORMAT (/,"  see  LOCATION  AND  RESUPPLY  PARAME 

TERS") 

r 

920= 

115 

FORMAT ("  - 

—  ■•,/> 

930= 

120 

FORMAT ( " 

ORBIT  INCLINATION  (DEG) 

=",F10.2) 

r. 

940= 

121 

FORMAT ( " 

ORBIT  RADIUS  (KM) 

=".F10.2) 

950= 

122 

FORMAT ( " 

NO  OF  ORB  LAUNCHES  (LS/YR) 

=",F10.2) 

'4 

960= 

123 

FORMAT ( " 

ORB  PAYLOAD  MASS  (KG/MSN) 

=",F10.2) 

:  ■ 

970= 

DO  130  1=1.1 

■' 

980= 

PRINT  60 

r  ' . 

990= 

PRINT  70 

1000= 

PRINT  80, B 

1010= 

PRINT  81,  H 

1020= 

PRINT  82, MOT V 

r. 

1030= 

PRINT  83,MS0C 

1040= 

PRINT  84,MCS0C 

1050= 

PRINT  35. I SPOT V 

» ■  _ 

1060= 

PRINT  86, ISPSOC 

• 

1070= 

PRINT  87, ISPORB 

4 

1080= 

PRINT  90 

—190= 

PRINT  95 

121 


1100= 

1110= 

1120= 

1 130= 
1140= 
1150= 
1160= 
1170= 
1180= 
1190= 
1200=  130 
1210= 
1220= 
1230= 
1240= 
1250= 
1260= 
1270= 
1280= 
1290= 
1300= 
1310= 
1320=C 
7330=C 
•"T^40=C 
^1350= 
1360=  140 
1370= 
1380=C 
1390=C 
1 4C0=C 
1410=  150 
l‘;20=  151 
1430= 
1440=  152 
1450= 
1460=  153 
1470= 
1480=  154 
1490= 
1500=  155 
1510= 
1520= 
1530=  156 
1540= 
1550= 
1560= 
1570= 

1 SS0=C 
*  -^i90=C 


PRINT  lOO.NFORB 
PRINT  101,MFSCC 
PRINT  102,MF0Ty 

PR I NT  1 03 , FORB«NFORB+MFSOC+MFOT V 

PRINT  no 

PRINT  115 
PRINT  120, X<1) 

PRINT  121, X (2) 

PRINT  122,F0RB 
PRINT  123, MPL 
CONTINUE 
END 

SUBROUTINE  RESTNT < IN, VAL) 

REAL  ASAT(35) ,ESAT(35> , ISAT(35) ,MSAT(35) , 
FSAT(35)  ,MFMSN(35)  ,PI,U,  W,P0,S,r1ET, 

MR, B, MOTV, MSOC. MCSOC, MFOTV, MFSOC, MFORB, 
MPL.FQRB, ISPOTV, ISPSOC, ISPGRB 
COMMON/ VALUES/ ASAT , ESAT, I SAT , MSAT , FSAT . MFMSN, 

P I . U , W , PQ , 6 , MET , MR , K , B , MOTV , MSOC , MCSOC . 
MFOTV, MFSbc,MFORB, MPL, FORB, ISPOTV, ISPSOC, 
ISPORB,H 

C0MM0N/SHARE/X(100),DEL<100) ,A<100, 100) ,N<5) 

X(1)=S0C  INCLINATION  :  X(2)=S0C  ALTITUDE 

IF (IN)  140,140,150 
CALL  STSOPt(VAL) 

RETURN 

INCLINATION,  ALTITUDE  AND  FUEL  CONSTRAINTS 

SOTO <151. 152, 153, 154, 155, 156) , IN 

VAL=X  <l)-56.0 

RETURN 

VAL=104.0-X (1) 

RETURN 

VAL=X (2) -6500.0 
RETURN 

VAL=7200.0-X(2) 

RETURN 
CALL  ORBFUL 
VAL= 10830. O-MFORB 
RETURN 

CALL  STSOPTIVAL) 

VAL=MPL-0.0 

RETURN 

END 

SUBROUTINE  STSOPT(VAL) 

CALCULATES  OTV  FUEL  CONSUMPTION  FOR  K  MISSIONS 


122 


1600=C 


1610= 

1620=  + 

1630=  + 

1640=  + 

1650=  + 

1 660= 

1670= 

1680= 

1690=  + 

1700=  + 

1710=  + 


1720= 

1730= 

1740= 

1750= 

1760= 

1770= 

1780= 

1790= 

1800= 

1810= 

1820=  + 

A830= 
*^40= 
^'a850= 

1860=  + 

1870= 

1880= 

1890=  160 
i  900= 

1910= 

1920=C 

1930=C 

1940=C 

1950= 

1960= 

1970= 

1980= 

1990=C 

20C0=C 

2010=C 

2020= 


2030= 

+ 

2040= 

+ 

2050= 

+ 

2060= 

+ 

2070= 

2080= 

+ 

-^90= 

+ 

REAL  I see , RSDC , MCSAT . VOT V , R AS AT , M . P , R , S , T , 

ASAT(35)  ,ESAT(35) , ISAT (35) , MSAT <35) , 

FSAT<35) ,MFMSN<35) ,PI,U, W,PO.S,MET, 

MR, B, MOTV, MSOC, MCSOC, MFOTV. MFSOC. MFORB, 
MPL,F0R8, I SPOT V, ISPSOC, ISPORB 
INTEGER  I 

COMMON/SHARE/X  (100)  .DEL(IOO)  .  AdOO,  100)  ,N  (5) 
COMMON/VALUES/ASAT, ESAT, ISAT, MSAT, FSAT, MFMSN, 

PI , U,  W,  PO, 0, MET, MR. K, B, MOTV, MSOC, MCSOC, 

MFOTV, MFSOC, MFORB, MPL,FORB. ISPOTV. ISPSOC, 
ISPORB, H 
MF0TV=0.0 
MCSAT=0.0 
ISOC=X (1) 

RSOC=X (2) 

DO  160  1=1, K 

RASAT=ASAT  <  I )  *  ( 1 . 0•^ESAT  ( I )  ) 

M=SQRT ( ABS (2. 0E6»U* ( 1 . O/RSOC-1 . 0/ (RASAT+REOC) ) ) ) 
P=SaRT (U* 1 . 0E6/RS0C) 

R=1 . 0E6*U» (4. 0/RASAT-l , 0/ASAT ( I ) -2. 0/ (RASAT+RSOC) ) 
S=4. 0E6*U*SaRT (ABS ( ( ( 1 . O/RASAT-0. 5/ASAT ( I ) ) 

« ( 1 . 0/RASAT-l . 0/ (RASAT+RSOC) ) )  )  ) 

T=COS ( ( I SOC- I SAT ( I ) ) *P I / 1 80 . 0 ) 
VOTV=ABS(M-P)+SQRT(ABS(R“S»T)  ) 

MFMSN ( I ) = (MSAT ( I ) +M0TV+M0TV4EXP ( VOTV/ ( ISP0TV4G) ) ) 
*(EXP(V0TV/(ISP0TV*6) )-l,0) 

MFOTV=MFOTV+MFMSN ( I ) *FSAT ( I ) 

MCSAT=MCSAT+FSAT ( I ) *MSAT ( I ) 

CONTINUE 
CALL  SOCFUL 
ChLL  OREFUL 

CALCULATES  TOTAL  ANNUAL  FUEL  CONSUMPTION 

FORB= (MCSAT+MFOTV+MFSOC+MCSOC) /MPL 
VAL=FORB 
END 

SUBROUTINE  ORBFUL 

CALCULATES  ORB I TER  FUEL  CONSUMPTION 

REAL  I SOC , RSOC , MM , VORBA , VORBD , MFORB A , MFORBD , 

ASAT(35) , ESAT (35) , ISAT(35) , MSAT (3S)  , 

FSAT (35) .MFMSN (35) , PI , U, W, PO, 6 , MET, 

MR, B, MOTV, MSOC, MCSOC, MFOTV, MFSOC, MFORB, 

MPL,FCRB, ISPOTV, ISPSOC, ISPORB 
COMMON/VALUES/ASAT, ESAT, ISAT, MSAT, FSAT, MFMSN, 

PI , U, W, PO, e. MET, MR, K, B, MOTV, MSOC, MCSOC. 

MFOTV, MFSOC, MFORB, MPL, FORB, ISPOTV, ISPSOC. 


2100=  +  ISPORB.H 

2110=  COMMON/SHARE/X  (100)  .DEL(IOO)  ,  AdOO,  100)  ,N  (5) 

2120=  ISDC=X(1) 


2130=  RSOC=X (2) 

2140=  MM= 164628. 209-79. 17494*ISCC-3. 01734*ISaC**2 

2150=  +  +0. 01256*IS0C**3 

2160=  V0RBA=-2. 023703587E-5*RSOC**2+0. 8474467398*RS0C 


2170=  + 

2180= 

2190=  + 

2200= 


2210= 

2220= 

2230= 

2240= 

2250= 

224)0= 

2270= 

2280=C 

2290=C 

2300=C 

2310= 


+ 


+ 


2320= 

+ 

7330= 

+ 

[ir4o= 

“2350= 

+ 

2360= 

2370= 

2380= 

+ 

2390= 

:400= 


2410= 


-4623. 658232+85. 65 

V0RBD=1 . 97461 01 71E-4»RS0C**2-2. 501373195*RS0C 
+7994.436774 

MPL=<MM-MET) *EXP (-VORBA/ ( ISPORB»G> ) 
-MR»EXP(VORBD/ (ISP0RB*6) ) 

MFORB A= ( MM-MET ) » ( 1 . 0-E  X  P ( - VQRB A / ( I SPORB  »  S ) ) ) 
MFORBD=  <MM-MET-MFORBA-MPL) 

*  ( 1 . 0-EXP (-VORBD/ ( ISPORB»G) ) ) 
MFORB=MFORBA+MFGRBD 
END 

SUBROUTINE  SOCFUL 

CALCULATES  SOC  STATIONKEEPING  FUEL  CONSUMPTION 

REAL  ISOC.RSOC.ASOC.F, 

ASAT(35) ,ESAT<35) , ISAT (35) , MSAT (35) , 
FSAT<35) ,MFMSN(35) ,PI,U,W,P0,6,MET, 

MR , B , MOTV , MSOC , MCSOC , MFOT V , MFSOC , MFORB , 
MPLjFORB, I SPOT V, ISPSOC, I SPORB 
COMMON/SHARE/X  ( 100) , DEL  (100)  ,  AdOO,  100)  ,N (5) 
COMMON/ VALUES/ ASAT , ESAT , I SAT , MSAT , FSAT » MFMSN , 
PI ,  U,  W,  PO,  6.  MET,  MR.  K,  B.  MOTV,  MSOC,  MCSOC, 
MFOTV, MFSOC, MFORB, MPL, FORB, I SPOT V, ISPSOC, 
I SPORB, H 
ISOC  =X  d) 


2420= 

2430= 

2440= 

2450=  + 

2460= 

2470= 

2480=*E0R 


RS0C=X(2) 

F=0. 00175* ISOC+0. 84004 
ASOC=(U*B*PO*F 

*EXP ( (6498. O-RSOC) /H) ) / (2000*RS0C) 
MFSOC= (MS0C*AS0C»31 557600) / ( ISPSOC*G) 
END  , 


2490=  HD AT A  N=2, M=6. X=56. , 6500. , NT (5) =1 , THETA0=1 . E-l2  *END 
2500=*E0R 


□ 


□  □  □□ 


Computer  Program  C-3 


ETR  Model  VI  Dual  Objective 
Minimization  Program 


100=QTS. T300, lOlOO, CMIOOOOO. 1820472, SP07S 
1 1 0=ATT ACH , SUMT , SUMT , I D=AF IT. 
120=LIBRARY,SUMT. 

1 3C=FTN ( L=0 , PL= 1 0000 )  . 

140=LS0. 

150=*E0R 

1 60=  PROGRAM  MA I N  <  I NPUT ,  OUTPUT ,  S(: 


170= 

130= 

190= 

200= 

210= 

220= 

230= 

240= 

250= 

260= 

270= 

280= 

290= 

300= 

310= 

320= 

*|^’30= 

^340= 

350=0 

360=0 

370=0 

330=0 

390=0 

40<.!=C 

410= 

420=0 

430=0 

440=0 

450= 

460= 

470=0 

480=0 

490=0 

500= 

510= 

520= 

530= 

540=0 

550=0 

560=0 

570= 

-  -580= 
j90= 


PROGRAM  MAIN  < INPUT, OUTPUT, SAVE, TAPE5= INPUT, 
TAPE6=0UTPUT, TAPE7=SAVE) 

REAL  ASAT(35) ,ESAT(35) , I SAT (35) ,MSAT(35) , ISPCHM, LOWTQF, 
FSAT(35)  ,MFMSN(35>  .PI,U,Pa,G.:lET, 

MR , B , MOT V . MSOC , MCSCC , MFOT V , MFSOC , MFCRB , 

MPL.FORB, isPOTV, ISPSOC, ISPORB, MSP, STEP, 

MTHR . PWR , MFT . THRUST , T I MRTN ( 35 ) , T I MDPL ( 35 ) , 

T I MTOT , NODT V . A VGTOF , R I NT , MFOT V A , MFOT VD , MOMOT V 
INTEGER  I,K,POINT,TOTPTS 
COMMQN/DEVC/NI , NO. NS 

COMMON/ VALUES /ASAT , ESAT , I SAT , MSAT , FS AT , MFMSN , 

P I , U , PO , G , MET , MR , K , S , MOT V , MSOC , MOSOO . I SPOHM . 
MFOTV.MFSdo.MFORB.MPL.FORB, ISPOTV, ISPSOC, MSP, 

I SPORB , H , MTHR , PWR , MFT , THRUST , T I MRTN , T I MDPL , 

T I MTOT , NOOTV , AVGTOF , R I NT , MFOT VA , MFOT VD , MOMOT V 
COMMON/SHARE/X  (100)  ,DEL(100)  ,  AdOO,  100)  ,N(5) 

NI=5 

NC=6 

NS=7 

ENTERS  LOW  VALUE  OF  SECOND  OBJECTIVE  FUNCTION 
(LOWTOF),  STEP  INTERVAL  BY  WHICH  IT  IS 
INCREMENTED  AND  THE  TOTAL  NUMBER  OF  POINTS 
EVALUATED. 

DATA  LOWTOF, STEP, T0TPTS/400. , -20. , 20/ 

ENTERS  CONSTANT  PROGRAM  DATA 

DATA  PI,U,P0,G,MET,MR/3. 1415926535,398601.2. 

24.9, a. 7, 38399.0, 85000.0/ 

ENTERS  PROBLEM  PARAMETER  DATA 

DATA  K , B , MSOC , MOSOO , MSP/ 7 ,.02,1 00000 . , 20000 .,10./ 

DATA  H. ISPOTV, ISPSOC, ISPORB/30. 0, 2900, 0, 2900. 0, 313.  0/ 
DATA  MTHR, PWR, MFT. THRUST/51. 36, 3. 06, 12.0. . 129/ 

DATA  MCMOtv, I SPOHM/ 2270. 0, 455. 6/ 

ENTERS  SATELLITE  TRAFFIC  MODEL  DATA 

DATA  ( ASAT ( I ) , I = 1 , 7 ) / 4 1 000 . , 20000 . , 6700 . . 65000 . , 

25000. , 12000. ,6700. / 

DATA  (ESAT (I) , 1=1 . 7) /O. , 0. , 0. ,0. ,0. 7,0. ,0. / 


600= 
610= 
620= 
630= 
640=C 
a50=C 
660=C 
670=C 
680= 
690= 
700=  2 
710= 
720= 
730=  10 
740= 
750=  15 
760= 
770= 
730= 
790=  20 
800= 
310=  30 
920= 
830=  35 
*^40= 
^-^50=  36 
860=  37 
870= 
880= 
390= 
900= 


DATA  <ISAT (I) ,1=1,7) /O. ,0. ,28.5,55. ,65. ,90. ,98. / 

DATA  ( MSAT ( I ) , I = 1 , 7 ) / 2000 . , 500 . , 25000 . , 1 500 . , 

1500. , 4500. , 8000. / 

DATA  (FSAT(I)  ,  1  =  1 , 7) /2.  25,  .  375,  .  375,  3.  0,  .  375, 2.  25, 3.  75/ 

LOOP  TO  INCREMENT  2ND  OBJ.  FUNCTION,  THE  ION  OTV 
AVERAGE  TIME  OF  FLIGHT  (AV6T0F) 

AV8T0F=L0WTGF 

DO  130  P0INT=1,T0TPTS 

FORMAT  <//,  lOX.  "ION  OTV  TIME  OF  FLIGHT  IS  ■•,F10.2) 
PRINT  2,AVGTCF 
CALL  SUMT 

FORMAT (//. 4X, "ETR  MODEL  VI  DUAL  OBJECTIVE  ", 
"MINIMIZATION  PROGRAM") 

FORMAT  ( 4 X,  "===  =====  ==  ====  =========  , 

PRINT  10 
PRINT  15 

FORMAT</,"  SATELLITE  TRAFFIC  MODEL") 

PRINT  20 

FORMAT (" - ") 

PRINT  30 

FORMATS "MSN", 3X, "ASAT",2X, "ESAT",2X, "ISAT",2X, 

" MSAT " , 2X , " FSAT " , 3X , " OTV " , 3X , " DEPLOY " , 2  X , " RETURN " ) 
F0RMAT(37X, "FUEL",3X, "TIME" , 4X , "TIME" ) 

FORMATdX,  "#",5X,  "KM",9X,  "DEG",4X,  "KG" ,  3X ,  "#/YR" , 

4X, "KG",4X, "DAYS",4X, "DAYS",/) 

PRINT  35 
PRINT  36 
PRINT  37 


910= 
920= 
930= 
940= 
950= 
960= 
970= 
980= 
990= 
1000= 
1010= 
1020* 
1030= 
1040* 
1050= 
1060= 
1070= 
1080= 
>  -090= 


40  FORMAT < 12, 2X,F7,0, 1X.F4.3, IX, F5. I, 1X.F6.0, 

+  1X,F4.2,  1X,F7.0, IX, F6. 1,2X,F6. 1) 

DO  50  1=1, K 

PRINT  40, I,ASAT(I) ,ESAT(I) , ) SAT (I )  , 

+  MSAT ( I ) , FSAT  < I ) , MFMSN ( I ) , 

+  TIMDPL  < I ) /86400. , TIMRTN  < I ) /86400. 


50 

CONTINUE 

60 

FORMAT ( / , " 

CONSTANT  PARAMETERS") 

70 

FORMAT ("  - 

- "  ) 

72 

FORMAT ( " 

ION  THRUSTER  MASS  (KG) 

—  ll 

,F10.2) 

73 

FORMAT  < " 

THRUSTER  FUEL  TANK  MASS  (KS) 

=  •* 

,F10.2) 

74 

FORMAT ( " 

SPECIFIC  MASS  (KG/KW) 

,F10.2) 

75 

FORMAT ( " 

THRUSTER  POWER  (KW/ THRUSTER) 

=  ** 

,F10.2) 

76 

FORMAT ( " 

ION  THRUSTER  THRUST  (MN) 

.F10.2) 

77 

FORMAT ( " 

NUMBER  OF  OTV  ION  THRUSTERS  (#) 

=  *• 

,F10.2> 

78 

FORMAT ( " 

EXPENDED  MASS  AT  ETS  (KG) 

,F10.2) 

79 

FORMAT ( " 

shuttle  reference  MASS  (KG) 

,F10.2) 

80 

FORMAT ( " 

BALLISTIC  COEFFICIENT  (MiK*2/KG) 

—  II 

,F10. 2) 

81 

FORMAT ( " 

SCALE  HEIGHT  (KM) 

=  ** 

,F10.2) 

▲ 


126 


A. 


1100= 

82 

1110= 

89 

1120= 

83 

1130= 

84 

1140= 

88 

1150= 

85 

1160= 

86 

1170= 

37 

1130= 

90 

1190= 

95 

1200= 

100 

1210= 

101 

1220= 

102 

1230= 

103 

1240= 

110 

1250= 

115 

1260= 

120 

1270= 

121 

1280= 

122 

1290= 

123 

1300= 

124 

1310= 

125 

1320= 

126 

1330= 

^^350= 

1360= 


1370= 
1380= 
1 390= 


1  400= 
1410= 
1 420= 
1430= 
1440= 
1450= 
1460= 
1470= 
1480= 
1490= 
1500= 
1510= 
1520= 
1530= 
1540= 
1550= 
1560= 
1570= 
1580= 
^90= 


FORMAT  ( *■ 
FORMAT  ( •• 
FORMAT  { 
FORMAT  ( 
FORMAT  < " 
FORMAT  < “ 
FORMAT  ( 
FORMAT  <  ■■ 
FORMAT  <  •• 
FORMAT  <  '• 
FORMAT ( " 
FORMAT  <  •• 
FORMAT  ( ■■ 
FORMAT { " 
FORMAT  < " 
FORMAT ( “ 
FORMAT ( “ 
FORMAT ( " 
FORMAT ( " 
FORMAT  ( 
FORMAT  ( •• 
FORMAT  ( '• 
FORMAT  <  •• 
PRINT  60 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


ION  OTV  MASS  (KS) 

CHEMICAL  OTV  MASS  <KS) 

SOC  MASS  (KS) 

ANNUAL  SOC  CARGO  MASS  <KS/YR) 
CHEMICAL  OTV  SPECIFIC  IMPULSE 
ION  OTV  SPECIFIC  IMPULSE  (SEC) 
SOC  SPECIFIC  IMPULSE  (SEC) 

ORB  SPECIFIC  IMPULSE  (SEC) 

FUEL  CONSUMPTION") 


> 


ORB I TER  FUEL  MASS  (KG/ MSN) 

SOC  FUEL  MASS  (KG/YR) 

OTV  FUEL  MASS  (KG/YR)  = 

TOTAL  STS  FUEL  MASS  (KG/YR) 

SOC  LOCATION  AND  RESUPPLY  PARAMETERS 


70 

72, 

73, 

74, 

75, 

76, 
78. 


ORBIT  INCLINATION  (DEG) 

ORBIT  RADIUS  (KM) 

NO  OF  ORB  LAUNCHES  (LS/YR) 

ORB  PAYLOAD  MASS  (KG /MSN) 

NO  OF  ION  PROPELLED  OTV’S  (#)  = 

AVG  OTV  TIME  OF  FLIGHT  (DAYS)  = 
PAYLOAD  MASS  TO  SOC  (KG)  = 


MTHR 

MFT 

MSP 

PWR 

THRUST* 1000.0 
MET 


■,F10, 

',F10. 

■,F10. 

*,F10. 

■,F10. 

’,F10. 

•,F10. 

■,F10. 


,F10. 

',F10. 

■jFlO. 

•.FIO. 

•) 

■) 

■ , FIO. 
'.FIO. 
•.FIO. 
'.FIO. 
•,F10. 
.FIO. 
.FIO. 


PRINT  79. MR 
PRINT  80,8 
PRINT  81, H 
PRINT  89,MCM0TV 
PRINT  83.MS0C 
PRINT  84,MCS0C 
PRINT  88, ISPCHM 
PRINT  85, I SPOT V 
PRINT  86, ISPSOC 
PRINT  87,ISP0RB 
PRINT  90 
PRINT  95 
PRINT  lOOjMFORB 
PRINT  101,MFS0C 

PRINT  102, MFOT V+FORB* ( MFOT V A+MFOT VD ) 

PRINT  103, FORB* (MFORB+MFOTVA+NFOTVD) +MFSOC+MFOTV 
PRINT  110 
PRINT  115 
PRINT  120, X(l) 


127 


w  w  M  w  r  j  ro  ro  nj  io  hj  m  w  ro  m  m  to  to 


1600= 

PRINT  121, X (2) 

1610= 

PRINT  77, X (3) 

1620= 

PRINT  122,F0RB 

1630= 

PRINT  125,TIMT0T/K 

1640= 

PRINT  82, MOT V 

1650= 

PRINT  124,N0aTV 

1660= 

PRINT  123, MPL 

1670= 

PRINT  126,MPL-MF0TVA 

1680= 

AVBTOF=AVQTOF+STEP 

1690=  130 

CONTINUE 

1700= 

END 

1710= 

SUSROUT I NE  RESTNT (IN, VAL ) 

1720= 

REAL  ASAT<35) ,ESAT<35) , ISAT (35) , MSAT (35) , ISPCHM, 

1730= 

+ 

FSAT(35) ,MFMSN(35) ,PI,U,PO,S,MET, 

1740= 

+ 

MR, B, MOTV, MSOC, MCSOC, MFOTO, MFSOC, MFORB, 

1750= 

MPL,FORB, ISPOTV, ISPSOC, ISPDRB,MSP, 

1760= 

+ 

MTHR , PWR , MFT , THRUST , T I MRTN ( 35  > , T I MDPL ( 35  > , 

1770= 

+ 

TIMTOT, NOOTV, AV6T0F, RINT, MFOTOA, MFOTOD, MCMOTO 

1780= 

COMMON/ VALUES / AS AT , ESAT , I SAT , MSAT . FS AT , MFMSN , 

1790= 

PI , U, PO. e, MET, MR. K. B, MOTV, MSOC, MCSOC, ISPCHM, 

1800= 

+ 

MFOTV, MFSOC, MFORB, MPL, FORB, ISPOTV, ISPSOC, MSP, 

1810= 

+ 

I SPCRB , H , MTHR , PWR , MFT , THRUST , T I MRTN , T I MDPL , 

1820= 

+ 

T I MTOT , NOOTV , A VSTOF .RINT, MFOTV A , MFOtvD , MCMOTV 

1830= 

COMMON/SHARE/XdOO)  ,DEL(i00)  ,  AdOO,  100)  ,N  (5) 

WdSO-C 

X(1)=S0C  INCLINATION  :  X(2)=S0C  RADIUS 

1860=C 

X(3)=NUMBER  OF  OTV  ION  THRUSTERS 

1370=C 

1380= 

IF(IN)  140,140.150 

1890=  140 

CALL  STSOPt(VAL) 

l‘?00= 

RETURN 

1910=C 

1920=C 

INCLINATION,  ALTITUDE  AND  FUEL  CONSTRAINTS 

1930=C 
1940=  150 

6CT0d51 , 152,  153, 154,  155,  156,  157,  158)  .  IN 

1950=  151 

VAL=X (l)-28.5 

1960= 

RETURN 

1970=  152 

VAL=57.0-X (1) 

1980= 

RETURN 

1990=  153 

VAL=X( 2) -6500.0 

2000= 

RETURN 

2010=  154 

VAL=X (3) -0.0 

2020= 

RETURN 

2030-  155 

CALL  ORBFUL 

2040= 

VAL= 10830. O-MFORB 

2050= 

RETURN 

2060=  156 

CALL  STSOPT(VAL) 

2070= 

VAL=MPL-MF0TVA-0. 0 

2080* 

RETURN 

*  -^»90=  157 

CALL  STSOPT(VAL) 
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2100= 
2110= 
2120=  158 
2130= 
2140= 
2150= 
2160= 
2170=C 
21S0=C 
2190=C 
2200= 
2210= 
2220= 
2230= 
2240= 
2250= 
2260= 
2270= 
2230= 
2290= 
2300= 
2310= 
2320= 
2330= 

2360= 

2370= 

2380= 

2390= 

2400= 

2410= 


2420= 
2430= 
2440= 
2450= 
2460= 
2470= 
2480= 
2490= 
2500= 
2510= 
2520= 
2530= 
2540= 
2550=  160 
2560= 
2570= 
2580* 
.'1590=0 
•y00=C 
2610=C 


+ 


+ 


VAL=N00TV-0.0 

RETURN 

CALL  STS0PT<VAL) 

VAL=T I MTOT/ K-A V6T0F 
RETURN 
END 

SUBROUT I NE  STSOPT  <  VAL ) 

CALCULATES  OTV  FUEL  CONSUMPTION  FOR  K  MISSIONS 

REAL  I SOC . RSOC , MCSAT , VOTV . RASAT , TEMP , 

ASAT.:35)  ,ESAT<35)  .  ISAT<35)  ,MSAT(35)  ,  ISPChM, 
FSAT<35) ,MFMSN(35) ,PI,U,PO,S.MET, 

MR , B . MOT V , MSOC , MCSCC , MFOT V . MFSOC , MFORB , 

MPL.FORB. I SPOT V, ISPSOC, ISPORBjMSP, 

MTHR , P WR , MFT , THRUST . MFRTN , MFDPL , 

TIMRTN(35) ,TIMDPL(35)  , 

TIMTOT, NQCtv, AVGTOF, RINT, MFOTVA, MFOTVD, MCMOTV 
INTEGER  I 

COMMON/SHARE/X (100) ,DEL(100) . A<100. 100) .N<5) 

COMMON / VALUES/ AS AT , ESAT , I SAT . MS AT . FS AT , MFMSN , 

P I , U , PO , G , MET , MR , K , B , MOT V , MSOC , MCSOC , 1 SPCHM , 
MFOTV, MFSOC, MFORB, MPLjFORB, I SPOT V, ISPSOC, MSP, 

I SPORB , H , MTHR , PWR , MFT , THRUST , T I MRTN , T I MDPL , 

T I MTOT , NOOTV , AVGTOF , R I NT , MFOT VA , MFOTVD , MCMOT V 
TIMT0T=0.0 
MF0TV=0.0 
MCSAT=0.0 
IS0C=X<1) 

R30C=X (2) 

MOTV=X  (3)  *  fMTHR+MSPCPWR-t-MFT)  +300. 

DO  160  1  =  1.  K 

VOTV=(ABS( <ASAT(I)/6378. 145) ** (-.5; 

-<RS0C/6378. 145) »«<-.5) ) 

+PI«ABS( (ISAT(I)-ISOC) CPI/ISO.O) 

/ <2.0»SQRT(ASAT(I) /637a. 145) ) >  *7905.36828 
TEMP=EXP(VOTV/ <ISP0TV*6) )-1.0 
MFRTN=MOTV*TEMP 

T I MRTN  < I )  =  ( MFRTN* I SPOT V  *6)/(X(3)*THRUST) 

MFDPL=  <MSAT  < I ) +MOTV+MFRTN) »TEMP 
TIMDPL<I)=(MFDPL»ISPOTV«G) / (X (3) «THRUST) 

MFMSN  < I ) =MFRTN+MFDPL 

TIMTOT=TIMTOT+(TIMDPL<I>+TIMRTN(I) ) /86400. 
MFOTV=MFOTV+MFMSN  < I ) *FSAT ( I > 

MCSAT=MCSAT+FSAT ( I ) *MSAT ( I ) 

CONTINUE 

NOOT V=T I MTOT / 365 . 25 
CALL  SOCFUL 
CALL  ORBFUL 

CALCULATES  OTV  FUEL  NEEDED  TO  RENDEZVOUS  WITH  ORB I TER 
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2620= 

2630= 

2640= 

2o50= 

2660= 

2670= 

2680=  + 

2690= 

2700= 

2710= 

2720= 

2730=C 

2740=C 

2750=C 

2760= 

2770=  + 

2780= 

2790= 

2800= 

2aio=c 

2820=C 

2S30=C 


,2840= 

1^50= 

”d60= 

+ 

+ 

2870= 

+ 

2880= 

+ 

28‘='0= 

+ 

2900= 

2^0  J 

2920= 

2930= 

2940= 

+ 

2950= 

2960= 

2970= 

2980= 

2990= 

+ 

3000= 

3010= 

+ 

3020= 

3030= 

3040= 

3050= 

+ 

3060= 

3070= 

3080= 

+ 

3090= 

IINT=X  <l; 

RINT=6538. 1 

M=SQRT<ABS<2.0E6»U»<1.0/ftINT-1.0/(RI^4T+RSCC)  )  )  ) 
P=3QRT  <U»1 . 0E6/RINT) 

R=1 . 0E6*U* (3. O/RSOC-2. 0/ (RSOC+RINT) > 

8=4. 0E6*U*SQRT ( ABS  < ( 1 . 0/ (2. 0*R20C) ) 

» (l.O/RSOC-1.0/ (RSOC+RINT) ) ) ) 

T=COS( (IINT-ISCC) *PI/180.0) 

VOTV=ABS  <  M-P ) +SQRT  <  ABS ( R-S* T ) ) 

MFOTVD=MCMOTV* (EXP (VOTV/ ( ISPCHM*G) ) -1 . 0) 

MFOTOA= (MCMOTV+MPL) * ( 1 . 0-EXP (-VOTO/ ( ISPCHM^B) ) ) 

CALCULATES  TOTAL  ANNUAL  FUEL  CONSUMPTION 

FORB= (MCSAT+MFOTV+MFSOC+MCSOC+ (MOTV»NOOTV) / 16. 0) / 
( MPL-MFOTOA-MFOTVD ) 

VAL=FORB 

END 

SUBROUTINE  ORBFUL 

CALCULATES  ORB I TER  FUEL  CONSUMPTION 

REAL  I I NT , R I NT , MM , VORB A , VORBD , MFORB A , MFORBD , 

ASAT (35) , ESAT (35) , ISAT (35) , MSAT (35) , ISPCHM, 
FSAT(35) ,MFMSN(35) ,PI,U.PO.e,MET, 

MR, B, MOTV, MSOC, MCSOC, MFOTV, MFSOC, MFORB, 
MPL,FORB, ISPOTV, ISPSOC, ISPORB.MSP, 
MTHR,PWR.MFT, THRUST, TIMRTN (35) , TIMDPL (35)  , 
TIMTOT, NCOTV, AVSTOF, MFOTVA, MFOTVD, MCMOTV 
CCMMON/VALUES/ ASAT, ESAT, iSAT, MSAT, FSAT, mFMSN, 

P I , U , PO , 3 , MET , MR , K . 3 . ^  OT V , MSGC . MCSOC , 1 SPCHM , 
MFOTV, MFSOC, MFORB, MPL.FQRB, ISPCTV, ISPSCC,MSP, 

I SPORB , H - MTHR , PWR , MFT , THRUST , T I MRTN .TIM  OPL , 
TIMTOT, NOOTV, AVSTOF, RINT. MFOTVA, MFCTVO. MCMQTV 
COMMON/SHARE/X (100) ,DEL(100) ,A(iOO, 100) ,N(5) 
IINT=X<1) 

RINT=653a. 1 

MM=159275. 829+3. 3575»I INT-3. 70994* I INT»»2 
+0.01339»IINT«*3 
V0RBA=51.7 
V0RBD=a4. 1 

MPL= (MM-MET) *EXP (-VORBA/ ( ISPCRB*S) ) 

-MR*EXP (VORBD/ ( ISPORB*S> ) 

MFORBA=( MM-MET) * ( 1 . 0-EXP (-VORBA/ (ISPORB*S) ) ) 
MFORBD= (MM-MET-MFORBA-MPL) 

♦( 1.0-EXP (-VORBD/ (ISPCRB*G) ) ) 
MFCRB=MFORBA+MFORBD 
END 
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»> 


i 


3100= 
3110= 
3120= 
3130= 
3140= 
3150= 
3160= 
3170= 
3130= 
3190= 
3200= 
3210= 
3220= 
3230= 
3240= 
%i^50= 
3260= 
3270= 
3280= 
3290= 
3300= 
3310= 
3320= 
3330= 
ni<40= 
^350= 
3360= 
3370= 
3380= 
3390= 
"400= 
34  !.0= 
3420= 
3430= 
3440= 
3450= 
34to0= 
3470= 
3480= 
3490= 
3500= 
3510= 
3520= 
3530= 
3540= 
3550= 


«E0R 

$DATA 

*DATA 

*DATA 

5DATA 

5DATA 

$DATA 

>DATA 

»DATA 

fDATA 

«DATA 

*DATA 

«0ATA 

♦DATA 

♦DATA 

♦DATA 

♦DATA 

♦DATA 

♦DATA 

♦DATA 


SUBROUTINE  SOCFUL 

CALCULATES  SOC  STATIONKEEPING  FUEL  CONSUMPTION 

REAL  ISOC.RSOC,ASOC,F,. 

ASAT(35) ,ESAT(35) , ISAT<35) ,MSAT<35) , ISPCHM, 
FSAT<35) ,MFM3N(35) , PI , U, PO, S, MET, 

MR , B , MOT V , MSOC , MCSOC , MFOT V , MFSOC , MFDRB , 

MPL, FORB, ISFOTV, ISPSOC, ISPORS,  MSP, 
MTHR,PWR,MFT,THRUST,TIMRTN(35) ,TIMDPL(35) , 
TIMTOT, NOOTV. A06T0F, RINT, MFOTOA, MFOTVD, MCMOTV 
COMMON/SHARE/X  <100)  ,DEL(i00)  ,  AdOO,  100)  ,N  (5) 
COMMON/ VALUES/ ASAT, ESAT , 1 3AT , MSAT , FSA  f , MFMSN , 

P I , U , PO , Q , MET , MR , K , B , MOT V , MSGC , MCSOC , I SPCHM . 
MFOTV, MFSOC, MFORB, MPL, FORB, ISFOTV, ISPSOC, MSP, 
I SPORB ,  H ,  MTHR , PWR , MFT , THRUST , T I MR t N , T I MDPL , 

T I MTOT . NOOT V , AVGTOF , R I NT , MFOTV A , MFOTVD . MCMOT V 
IS0C=X(1) 

RSOC=X (2) 

F=0 .00175*1 SOC+0 . 84004 
ASOC=<U*B*PO*F 

*EXP< (6498. O-RSOC) /H) ) / (2000»RS0C) 

MFSOC=»  (MS0C*AS0C*31557600)  /  ( ISPS0C*6) 

END  ^ 

N*3 . M»7 , MZ = 1 , X=2a . , 6500 . , 1 . , NT ( 5 ) = 1 , THET A0» 1 . E-3  ♦END 
N=3,M=7,MZ=l,X=2a. ,6500. , 1 . , NT (5) *1 , THETAO*! . E-3  ♦END 
Ns3,M*7.MZ=1,X=28. ,6500. , 1 . , NT (5) =1 , THETA0=1 . E-3  ♦END 
N4=3,  M*7,  MZ=1  .  X=23.  ,  6500.  ,  1 .  ,  NT  (5)  =1 ,  THETA0=1 .  E-3  ♦END 
N=3,  M=7.  MZ=1  ,  X=2S.  ,  6500.  ,  1 .  ,  NT  (5)  -1 ,  thETA0=1  .  E-3  ♦£,%£■ 
■M=3,M=7,  MZ  =  1,X=2S.  ,6500.  ,  1,  ,  NT  <5)  =  1 ,  THETA0=1  .  E-3  SEND 
,\=3.  i  MZ=1 .  X=2S,  o500.  ,  I ,  ,  NT  <5)  -i  ,  .'HETAu^i  .  E-3  -END 
N=3,  M=7.  MZ=1 ,  X=28.  ,  o500.  ,  1 .  ,  NT  (5)  -=1 ,  "HETAO=l .  E-3  iEND 
N=3,M=7,MZ=1, X=28. ,6500. , 1 ., NT (5) =1 . ThETAO=l . E-3  SEND 

♦END 


N=3,M=7,,MZ  =  1,  X=28.  ,6500. 
N=3,M=7,MZ=1,X=23. ,6500. 
N=3,M*7,MZ»1, X=28. ,6500. 


,  NT (5) =1 , thETA0=1 . E-: 


,  NT  ( 5 )  a  1 ,  THETAOa  1 .  E-3  «E:ND 


N=3,Ma7,MZ=l,X=28. ,6500. , 1 . , NT (5) =1 , THETA0=1 . E-3  ♦END 
N=3, Ma7. MZ=1 , X=28. , 6500. . 1 . , NT  <5) =1 , THETA0=1 . E-3  ♦END 
N=3,M»7,MZ=1,X=28. ,6500. . 1 . , NT (5) =1 , THETA0=1 . E-3  ♦END 
N=3,M=7.MZ=1, X=28. ,6500. , 1 . , NT (5) =1 , THETA0=1 . E-3  SEND 
N=3, M=7, MZ=1 , X=28. . 6500. , 1 . , NT (5) =1 , THETA0=1 . E-3  SEND 
N=3.M=7,MZ=1, X=28. ,6500. . 1 . . NT (5) =1 , THETA0=1 . E-3  SEND 
N=3,M=7,MZ=1,X=28. ,6500. . 1 , , NT (5) *1 , THETA0=1 . E-3  ♦END 
N=3,M=7,MZ=1,X=28. .6500. , 1 . , NT (5) =1 . THETA0=1 .  E-3  ♦END 
N=3,Ma7,MZ=l,X=28. ,6500. , 1 . , NT (5) =1 , THETA0=1 . E-3  ♦END 
N=3,M=7,MZ=1,X=28. ,6500. , 1 . , NT (5) =1 . THETA0=1 . E-3  ♦END 


131 


APPENDIX  D 


Operational  Computer  Runs 

Over  a  hundred  computer  runs  were  accomplished  for  this 
thesis  effort.  Table  D-1  lists  all  of  the  runs  and  the 
configuration  of  the  Model  during  the  run.  The 
configuration  includes  the  particular  scenario  considered, 
the  number  of  satellite  missions  serviced  and  any  parameters 
that  are  varied.  In  addition,  the  output  of  selected 
computer  runs  are  included  for  those  requiring  a  more 
detailed  understanding  of  the  different  models  and 
scenarios.  The  minor  MECO  mass  sign  error  discussed  in 
Appendix  C  was  corrected  in  all  of  the  attached  computer 
runs.  Thus,  the  attached  Model  II  and  III  results  are  more 
accurate  than  the  results  of  Section  4.0  although  the 
differences  are  very  small. 


Table  D-1  Computer  Runs  Accomplished 


9 


3 

i  tr 

■ 


Run 

# 

Model 

# 

Launch 

Range 

soc 

Scen¬ 

ario 

Satellite 

Traffic 

Model 

Varied 

Para¬ 

meters 

1 

I 

ETR 

C 

1-5 

None 

2 

rt 

II 

A 

1-7 

II 

3 

II 

If 

D 

4-5 

If 

4 

If 

ETR 

B,D 

1-3 

11 

5 

I 

WTR 

C,D 

6-7 

It 

6 

It 

II 

B 

4-7 

It 

7 

If 

If 

A 

1-7 

It 

8 

11 

If 

D 

4-5 

II 

9 

II 

ETR 

A 

1-7 

II 

10 

tt 

If 

B,D 

1-3 

11 

11 

II 

If 

D 

4-5 

II 

12 

w 

If 

C 

1-5 

It 

13 

14 

15 

16 

II 

II 

If 

If 

WTR 

If 

If 

II 

B 

D 

C,D 

A 

4-7 

4-5 

6-7 

1-7 

11 

II 

ft 

II 

17 

II 

ETR 

A 

1-7 

B 

=  .01 

18 

II 

If 

If 

II 

If 

=  .03 

19 

II 

It 

It 

II 

H 

=  15 

20 

11 

ft 

II 

It 

It 

s  45 

21 

II 

ETR 

A 

1-7 

Motv 

S  1135 

22 

11 

It 

ft 

II 

II 

=  3405 

23 

If 

If 

II 

11 

^soc 

=  50000 

24 

It 

If 

II 

II 

It 

=  100000 

25 

II 

ETR 

A 

1-7 

SOC 

=  10000 

26 

II 

II 

II 

n 

It 

=  30000 
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Table  D-1  Computer  Runs  Accomplished  (Cont.) 


Run 

Model 

Launch 

soc 

Satellite 

Varied 

# 

# 

Range 

Seen- 

Traffic 

Para- 

ario 

Model 

meters 

53 

II 

ETR 

A 

1-7 

Mr 

79236.9 

54 

II 

II 

tl 

II 

II 

= 

96845.1 

55 

III 

ETR 

A 

1-7 

None 

56 

II 

II 

C 

1-5 

It 

57 

II 

WTR 

A 

1-7 

II 

58 

II 

II 

C 

6-7 

11 

59 

III 

ETR 

A 

1-7 

”otv 

1135 

60 

II 

II 

11 

It 

II 

5 

3405 

61 

II 

II 

II 

If 

ISPotv 

s 

227.5 

62 

n 

11 

tl 

If 

= 

682.5 

63 

III 

ETR 

A 

1-7 

ISPotv 

• 

2900 

64 

It 

II 

11 

II 

It 

s 

9000 

65 

II 

tl 

It 

If 

ISPorb 

- 

156.5 

66 

II 

II 

II 

II 

n 

r 

469.5 

67 

III 

ETR 

A 

1-7 

Mr 

79236.9 

68 

II 

II 

II 

II 

tl 

= 

96845.1 

69 

IV 

ETR 

A 

1-7 

ISP  Qtv 

• 

9000 

70 

71 

II 

II 

II 

II 

II 

If 

If 

tl 

It 

II 

* 

29o8°°  ' 

5000 

72 

11 

II 

If 

II 

II 

ISPsoc  = 

9000 

73 

IV 

ETR 

A 

1-7 

ISPotv 

ISPsoc  “ 

2900 

74 

11 

II 

If 

ft 

II 

= 

II  _ 

10000 

75 

II 

n 

If 

It 

II 

It  - 

30000 

76 

II 

n 

tl 

II 

II 

If  - 

20000 

77 

IV 

ETR 

A 

1-7 

None 

78 

II 

11 

B 

1-3 

II 
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Table  D-1  Computer  Runs  Accomplished  (Cont.) 


Run 

# 

Model 

# 

Launch 

Range 

SOC 

Scen¬ 

ario 

Satellite 

Traffic 

Model 

Varied 

Para¬ 

meters 

79 

IV 

ETR 

C 

1-5 

None 

80 

II 

WTR 

A 

1-7 

II 

81 

II 

It 

B 

4-7 

11 

82 

II 

If 

C 

6-7 

II 

83 

III 

ETR 

B 

1-3 

II 

84 

11 

WTR 

B 

4-7 

It 

85 

86 

11 

III 

II 

II 

11 

ETR 

ff 

It 

ff 

A 

ff 

ff 

If 

1-7 

If 

It 

If 

Mr 

n 

ISP  orb 

79236.9 

96845.1 

156.5 

469.5 

89 

III 

ETR 

A 

1-7 

Thrust 

.1935 

90 

If 

ff 

If 

ft 

ti 

.0645 

91 

If 

ff 

If 

ff 

N 

10 

92 

If 

11 

IT 

If 

M 

30 

93 

94 

95 

96 

III 

If 

II 

If 

ETR 

If 

If 

If 

A 

ft 

II 

ft 

1-7 

11 

II 

If 

Mthr 

11 

Pwr 

‘^wr 

=  25.68 

=  77.04 
=  1.53 
=  4.59 

97 

III 

ETR 

A 

1-7 

”SD 

=  5 

98 

tf 

ft 

It 

ft 

sp 

ft 

=  15 

99 

V 

ETR 

A 

1-7 

None 

100 

It 

If 

B 

1-3 

II 

101 

It 

ft 

C 

1-5 

11 

102 

If 

WTR 

A 

1-7 

It 

103 

If 

tf 

B 

4-7 

If 

104 

If 

If 

C 

6-7 

If 
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Table 

D-l  Computer  Runs 

Accomplished 

( Cont . ) 

Run 

# 

Model 

# 

Launch 

Range 

soc 

Scen¬ 

ario 

Satellite 

Traffic 

Model 

Varied 

Para¬ 

meters 

105 

VI 

ETR 

A 

1-7 

Standard  Parameters 
T0F=300,280,  . . . 

106 

II 

It 

ti 

II 

Standard  Parameters 
TOF=300,320,  .  .  . 

107 

II 

WTR 

II 

II 

Level  I  Technology 
TOF=400,380,  . . . 

108 

II 

ETR 

II 

II 

Level  I  Technology 
TOF=400,380, .  . . 

109 

II 

II 

B 

1-3 

Level  I  Technology 
Isoc=28.50 

TOF=400,380,  . . . 

110 

fl 

11 

A 

1-7 

Level  II  Technology 
TOF=400,380, .  . . 

in 

n 

II 

A 

II 

Level  III  Technologj 
TOFr400,380, .  . . 

112 

II 

It 

B 

4-7 

Level  I  Technology 
Isoc=57.0O 

T0F  =  400,380,  . . . 

113 

II 

WTR 

B 

4-7 

Level  I  Technology 
TOF=400,380, . . . 

Computer  Run  D-1  Run  Number  9 


1040= 

1050= 

1060= 

1070= 

1080= 

1090= 

1100= 

1110= 


ETR  MODEL  II  LAUNCH  MINIMIZATION  PROGRAM 


SATELLITE  TRAFFIC  MODEL 


1120= 

MISSION 

A3AT 

ESAT 

ISAT 

MSAT 

FSAT 

□TV  FUEL 

1130= 

1140= 

<#) 

(KM) 

(DEG) 

(KG) 

(LS/YR) 

(KG) 

1150= 

1 

41000. 

0.  00 

0.  0 

2000. 

3 . 00 

31577. 

1160= 

2 

20000. 

0.00 

0.0 

500. 

.  SO 

36754. 

1170= 

3 

6700. 

0.00 

28.5 

25000. 

.50 

53592. 

1130= 

4 

65000. 

0.00 

55.  0 

1500. 

4.00 

1 8370. 

1190= 

5 

25000. 

.  70 

65.0 

1500. 

.  50 

7953- 

1  200= 

6 

12000. 

0.00 

90.0 

4500. 

si .  00 

24841. 

1210= 

1220= 

7 

6700. 

0.00 

98.0 

8000. 

5.  00 

55980. 

1230= 

1240= 


VARIABLE  PARAMETERS 


1250= 

(^60= 

BALLISTIC  COEFFICIENT  (M»»2/KG) 

s 

.02 

L^70» 

SCALE  HEIGHT  (KM) 

s 

30.00 

1280= 

OTV  MASS  (KG) 

s 

2270.00 

1290= 

SOC  MASS  (KG) 

s 

100000.00 

1300= 

ANNUAL  SOC  CARGO  MASS  (KG/YR) 

= 

20000. 00 

1310= 

OTV  SPECIFIC  IMPULSE  (SEC) 

= 

455.00 

1320= 

SOC  SPECIFIC  IMPULSE  (SEC) 

zs 

455.00 

1330= 

ORB  SPECIFIC  IMPULSE  (SEC) 

=' 

w*  1  ■  00 

1340= 

1350= 

FUEL  CONSUMPTION 

1360= 

1370= 

1380= 

ORBITER  FUEL  MASS  (KG/MSN) 

8277.73 

1390= 

SOC  FUEL  MASS  (KG/YR) 

= 

35629.92 

1400= 

OTV  FUEL  MASS  (KG/YR) 

s 

571779.56 

1410= 

1420= 

1430= 

1440= 

1450= 

1460= 

1470= 

1480= 

1490= 


TOTAL  STS  FUEL  MASS  <KG/YR) 


SOC  LOCATION  AND  RESUPPLY  PARAMETERS 


992732. 79 


ORBIT  INCLINATION  <DE0) 
ORBIT  RADIUS  <KM) 

NO  OF  ORB  LAUNCHES  (LS/YR) 
ORB  PAYLOAD  MASS  (KG/MSN) 


=  57.00 
=  6670.37 
=  46.35 
=  15175.57 


I 
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Computer  Run  D- 2  Run  Number  10 


1500=  ETR  MODEL  II  LAUNCH  MINIMIZATION  PROGRAM 
1510=  ===  =====  ==  ======  ============  ======= 

1520= 

1530= 

1540=  SATELLITE  TRAFFIC  MODEL 

1550=  - 

1560= 

1570= 


1580= 

MISSION 

ASAT 

ESAT 

I  SAT 

MSAT  FSAT  OTV  FUEL 

1590= 

(4)) 

(KM) 

(DEG) 

(KG)  (LS/YR)  (KG) 

1600= 

1610= 

1 

41000. 

0.00 

0.0 

2000.  3.00  20820. 

1620= 

20000. 

0.  00 

0.0 

500.  .50  15542. 

1630= 

6700. 

0.00 

28.5 

25000.  .50  126. 

1640= 

1650= 

VARIABLE  PARAMETERS 

1660= 

1670= 

1680= 

BALLISTIC 

:  COEFFICIENT 

(M**2/KS) 

=  .02 

1690= 

SCALE  HEIGHT  (KM) 

=  30.00 

j--^00= 

OTV  MASS 

(KG) 

=  2270.00 

SOC  MASS 

(KG) 

=  100000.00 

1720= 

ANNUAL  SOC  CARGO 

MASS 

(KG/YR) 

=  20000.00 

1730= 

OTV  SPECIFIC  IMPULSE 

(SEC) 

=  455.00 

1740= 

SOC  SPECIFIC  IMPULSE 

(SEC) 

=  455.00 

1 750= 

ORB  SPECIFIC  IMPULSE 

(SEC) 

=  313.00 

1760= 

1770=  FUEL  CONSUMPTION 

1780=  - 

1 790= 

1300=  ORBITER  FUEL  MASS  (KG/MSN)  =  10278.08 

1810=  SOC  FUEL  MASS  <KS/YR)  =  4669.75 

1820=  OTV  FUEL  MASS  (KG/YR)  =  70294.70 

1830=  TOTAL  STS  FUEL  MASS  (KG/YR)  =  133549. 18 

1840= 

1850=  SOC  LOCATION  AND  RESUPPLY  PARAMETERS 

1860=  - 

1870= 

1880=  ORBIT  INCLINATION  (DEG)  =  28.50 

1890=  ORBIT  RADIUS  (KM)  =  6729.44 

1900=  NO  OF  ORB  LAUNCHES  (LS/YR)  =  5.70 

1910-  ORB  PAYLOAD  MASS  (KQ/MSN)  =  19950.01 


Computer  Run  D-3  Run  Number  55 


100=  ETR  MODEL  III  LAUNCH  MINIMIZATION  PROGRAM 

120= 

130= 

140=  SATELLITE  TRAFFIC  MODEL 

150= - 

160= 

170= 


180= 

MISSION 

ASAT 

ESAT 

I  SAT 

MSAT 

FSAT 

OTV  FUEL 

190= 

(#) 

(KM) 

(DEG) 

(KG) 

(LS/YR) 

(KG) 

200= 

210= 

1 

41000. 

0.00 

0.0 

2000. 

3 . 00 

31279. 

220= 

o 

20000. 

0.00 

0.0 

500. 

.  50 

36417. 

230= 

s 

W 

6700. 

0.00 

28.5 

25000. 

.50 

53733. 

240= 

4 

65000. 

0.  00 

55.0 

1500. 

4 . 00 

18148. 

250= 

5 

25000. 

.70 

65.  0 

1500. 

.  50 

7851. 

260= 

6 

12000. 

0.00 

90.0 

4500. 

3.00 

24593. 

270= 

7 

6700. 

0.00 

98.0 

8000. 

5.00 

56177. 

280= 

290=  VARIABLE  PARAMETERS 

300= - 

310= 

ii|^20=  BALLISTIC  COEFFICIENT  <M»*2/KG)=  .02 

"-^30=  SCALE  HEIGHT  <KM)  =  30.00 

340=  OTV  MASS  (KG)  =  2270.00 

350=  SOC  MASS  (KG)  =  100000.00 

360=  ANNUAL  SOC  CARGO  MASS  (KG/YR>  =  20000.00 

370=  OTV  SPECIFIC  IMPULSE  (SEC)  =  455.00 

380=  see  SPECIFIC  IMPULSE  (SEC)  =  455.00 

390=  OPB  SPECIFIC  IMPULSE  (SEC)  =  313.00 

400= 

410=  FUEL  CONSUMPTION 

420= - 

430= 

440=  ORBITER  FUEL  MASS  (KG/MSN)  =  4859.20 

450=  SOC  FUEL  MASS  (KQ/YR)  =  4189.47 

460=  OTV  FUEL  MASS  <KG/YR)  =  595837.88 

470=  TOTAL  STS  FUEL  MASS  (KG/YR)  =  782239.98 


480= 

490=  SOC  LOCATION  AND  RESUPPLY  PARAMETERS 

500= - 

510= 

520=  ORBIT  INCLINATION  (DEG)  =  57.00 

530=  ORBIT  RADIUS  (KM)  =  6734.31 

540=  NO  OF  ORB  LAUNCHES  (LS/YR)  =  37.50 

550=  ORB  PAYLOAD  MASS  (KG/ MSN)  =  18641.47 

560=  RENDEZVOUS  INCLINATION  (DEG)  =  56.84 

570=  RENDEZVOUS  RADIUS  (KM)  =  6538.10 

580=  PAYLOAD  MASS  TO  SOC  (KG)  =  18024.02 


Computer  Run  D-4  Run  Number  83 


::90=  ETR  MODEL  III  LAUNCH  MINIMIZATION  PROGRAM 

610= 

620= 

630=  SATELLITE  TRAFFIC  MODEL 

640= - 

650= 


660= 


670= 

MISSION 

ASAT 

ESAT 

ISAT 

MSAT 

FSAT 

□TV  FUEL 

680= 

<#) 

(KM) 

(DEG) 

(KG) 

(LS/YR) 

(KG) 

690= 

700= 

1 

41000. 

0.00 

0.  0 

2000 . 

3.00 

20443. 

710= 

20000. 

0.00 

0.0 

500. 

.50 

15229. 

720= 

3 

6700. 

0.00 

28.5 

25000. 

.50 

550. 

730= 

740=  VARIABLE  PARAMETERS 

750= - 

760= 


770= 

BALLISTIC  COEFFICIENT  (M*»2/KG) 

• 

780= 

SCALE  HEIGHT  (KM) 

= 

30.00 

790= 

OTV  MASS  (KG) 

=s 

2270.00 

800= 

SOC  MASS  (KG) 

100000.00 

810= 

ANNUAL  SOC  CARGO  MASS  (KG/YR) 

= 

20000.00 

320= 

OTV  SPECIFIC  IMPULSE  (SEC) 

455.00 

830= 

SOC  SPECIFIC  IMPULSE  (SEC) 

= 

455. 00 

ORB  SPECIFIC  IMPULSE  (SEC) 

313.00 

860= 
0*7 1^— 

FUEL  CONSUMPTION 

8SC» 

890= 

ORBITEP  FUEL  MASS  (KG/MSN) 

4985.72 

9o<:  = 

SOC  FUEL  MASS  (KG/YR) 

•= 

169.53 

910= 

OTV  FUEL  MASS  (KG/YR) 

= 

74306.78 

920= 

930= 

TOTAL  STS  FUEL  MASS  (KG/YR) 

97433.88 

940= 

SOC  LOCATION  AND  RESUPPLY  PARAMETERS 

950= 

960= 

970= 

ORBIT  INCLINATION  (DEG) 

28.50 

980= 

ORBIT  RADIUS  (KM) 

=5 

6823.47 

990= 

NO  OF  ORB  LAUNCHES  (LS/YR> 

= 

4.51 

1000= 

ORB  PAYLOAD  MASS  (KG/MSN) 

S 

25242. 37 

1010= 

RENDEZVOUS  INCLINATION  (DEG) 

s 

23.50 

1020= 

RENDEZVOUS  RADIUS  (KM) 

s 

6538. 10 

1030= 

PAYLOAD  MASS  TO  SOC  (KG) 

= 

24100.54 

p 


h 


5 


Computer  Run  D-5  Table  4.15  Efficient  Operating  Point 


-.00= 

lio-- 

120= 

130= 

140= 

ETR  MODEL  VI 

DUAL  OBJEcflvi’\?:^lVlSP^ 

SATELLITE 

TRAFFIC  MODEL 

150=MSN 

160= 

ASAT 

ESAT 

I  SAT 

MSAT 

FSAT 

OTV 

FUEL 

DEPLOY 

TIME 

RETURN 

TIME 

170= 

180= 

# 

KM 

DEG 

KG 

#/YR 

KG 

DAYS 

DAYS 

190= 

1 

41000. 

.  000 

0.0 

2000. 

2.25 

6051. 

125.7 

76.5 

200= 

2 

20000. 

.000 

0.0 

500. 

.  38 

5415. 

105.9 

75.0 

210= 

3 

6700. 

.000 

28.5 

25000. 

.38 

1559. 

43.  1 

9.0 

220= 

4 

65000. 

.000 

55.0 

1500. 

3.00 

5171. 

104.5 

68.3 

230= 

5 

25000. 

.700 

65.0 

1500. 

.38 

5658. 

115.  1 

73.9 

240= 

6 

12000. 

.000 

90.0 

4500. 

2.25 

11673. 

268.9 

1  04  O 

4  ^  X  a  ^ 

250= 

260= 

-7 

0 

6700. 

.000 

98.0 

8000. 

3.75 

18938. 

471.6 

161.3 

270 

280 

290 

300 

310 

320 

330 

340 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

vT/ilO 


=  CONSTANT  PARAMETERS 


(KS) 


ION  THRUSTER  MASS  (KG) 

THRUSTER  FUEL  TANK  MASS 
SPECIFIC  MASS  <KG/KW> 

THRUSTER  POWER  (KW/THRUSTER)  = 
ION  THRUSTER  THRUST  <MN) 

EXPENDED  MASS  AT  ETS  (KG) 

SHUTTLE  REFERENCE  MASS  (KG) 
BALLISTIC  COEFFICIENT  (M*»2/KB>= 
SCALE  HEIGHT  (KM) 

CHEMICAL  OTV  MASS  (KG)  * 

SOC  MASS  (KG) 

ANNUAL  SOC  CARGO  MASS  (KG/YR)  = 
CHEMICAL  OTV  SPECIFIC  IMPULSE  = 
ION  O’V  SPECIFIC  IMPULSE  (SEC)  = 
SOC  SPECIFIC  IMPULSE  (SEC) 

CRB  SPECIFIC  IMPULSE  <'SEC) 

FUEL  CONSUMPTION 


ORBITER  FUEL  MASS  (KG/MSN) 

SOC  FUEL  MASS  (KG/YR) 

OTV  FUEL  MASS  (KG/YR) 

TOTAL  STS  FUEL  MASS  (KG/YR) 

SOC  LOCATION  AND  RESUPPLY  PARAMETERS 

ORBIT  INCLINATION  (DEG) 

ORBIT  RADIUS  (KM) 

NUMBER  OF  OTV  ION  THRUSTERS  (#)» 
NO  OF  ORB  LAUNCHES  (LS/YR)  = 

AVQ  OTV  TIME  OF  FLIGHT  (DAYS)  = 

ION  OTV  MASS  (KB) 

NO  OF  ION  PROPELLED  OTV»S  (#)  = 

ORB  PAYLOAD  MASS  (KG/MSN)  s 

PAYLOAD  MASS  TO  SOC  (KG)  = 
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51.36 
12.00 
10.00 
3.06 
129.00 
38399.00 
85000.00 
.02 
30.00 
2270.00 
100000.00 
20000.00 
455.00 
2^0i) .  '>0 
2900 . 00 
313.00 


4874.09 
7258.64 
136096.27 
183132. 19 


33.  12 
6661.22 
67.74 
8.  16 
260.00 
6665.01 
4.98 
27531.12 
26968.35 


n 


Computer  Run  D-6  Table  4.15  Efficient  Operating  Point  for 

WTR  Scenario  A,  Missions  1-7,  Level  I 
Technology,  Run  Ntjmber  107 


■-X>0 

"1 10 
120 
130 


WTR  MODEL  VI  DUAL  OBJECTIVE  MI^4IMI ZATION  PROGRAM 


=  SATELLITE  TRAFFIC  MODEL 


140= 

— 

- — 

— - 

— 

150=MSN 

160= 

ASAT 

ESAT 

I  SAT 

MSAT 

FSAT 

OTV 

FUEL 

DEPLOY 

TIME 

RETURN 

TIME 

170= 

180= 

# 

KM 

DEG 

KG 

#/YR 

KG 

DAYS 

DAYS 

190= 

1 

41000. 

.000 

0.0 

2000. 

2.25 

7451. 

139.7 

104.4 

200= 

2 

20000. 

.000 

0.0 

500. 

.38 

7610. 

184.9 

115.4 

210= 

3 

6700. 

.  000 

28.5 

25000. 

.33 

11088. 

372.2 

65.5 

220= 

4 

65000. 

.000 

55.0 

1500. 

3.00 

3477. 

32.5 

54.7 

230= 

5 

25000. 

.  700 

65.0 

1500. 

.38 

2641. 

61.7 

42.5 

240= 

6 

12000. 

.000 

90.0 

4500. 

2.25 

5500. 

147.  1 

70.0 

250= 

260= 

7 

6700. 

.000 

98.0 

8000. 

3.75 

8349. 

242.6 

36.9 

)=  CONSTANT  PARAMETERS 


ION  THRUSTER  MASS  (KG) 

THRUSTER  FUEL  TANK  MASS  (KG) 
SPECIFIC  MASS  (KG/KW) 

THRUSTER  POWER  (KW/THRUSTER) 

ION  THRUSTER  THRUST  (MN) 

EXPENDED  MASS  AT  ETS  (KG) 

SHUTTLE  REFERENCE  MASS  (KG) 
BALLISTIC  COEFFICIENT  (M»»2/KB>= 
SCALE  HEIGHT  (KM) 

CHEMICAL  OTV  MASS  (KB) 
sac  MASS  (KG) 

ANNUAL  SOC  CARGO  MASS  (KG/YR)  = 
CHEMICAL  OTV  SPECIFIC  IMPULSE  = 
ION  OTV  SPECIFIC  IMPULSE  (SEC)  = 
SOC  SPECIFIC  IMPULSE  (SEC) 

CRB  SPECIFIC  IMPULSE  (SEC) 

FUEL  CONSUMPTION 


ORB I TER  FUEL  MASS  (KG/MSN) 

SOC  FUEL  MASS  (KG/YR) 

OTV  FUEL  MASS  (KG/YR)  = 

TOTAL  STS  FUEL  MASS  (KG/YR) 

SOC  LOCATION  AND  RESUPPLY  PARAMETERS 

ORBIT  INCLINATION  (DEG) 

ORBIT  RADIUS  (KM) 

NUMBER  OF  OTV  ION  THRUSTERS  (#)= 
NO  OF  ORB  LAUNCHES  (LS/YR> 

AVG  OTV  TIME  OF  FLIGHT  (DAYS)  = 

ION  OTV  MASS  (KG) 

NO  OF  ION  PROPELLED  OTV'S  (#)  = 

ORB  PAYLOAD  MASS  (KG/MSN) 

PAYLOAD  MASS  TO  SOC  (KG)  > 


51.36 
12.00 
10.  00 
3.06 
129.00 
38399.00 
35000.00 
.02 
30.00 
2270.00 
100000.00 
20000 . 00 
455.00 
2900. 00 
29C.)0 . 00 
313.00 


8255.52 
1588.41 
83950.63 
153014. 12 


59.  (.)0 
6708.07 
57.36 
8.  17 
260.00 
5689.28 
4.98 
20378.57 
19816.11 
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Computer  Run  D-7  Table  4.15  Efficient  Operating  Point 

for  ETR  Scenario  B,  Missions  1-3,  Level  I 
Technology,  Run  Number  109 


00=  ETR  MODEL  VI  DUAL  OBJECTIVE  MIh4IMI ZATIDN  PR06RAM 
120= 

130=  SATELLITE  TRAFFIC  MODEL 


140= 


150=MSN 

160= 

ASAT 

ESAT 

I  SAT 

MSAT 

FSAT 

□TV 

FUEL 

DEPLOY 

TIME 

RETURN 

TIME 

1 70= 
130= 

# 

KM 

DEG 

K3 

#/'yR 

KG 

BAYS 

DAYS 

190= 

1 

41000. 

.  000 

0.0 

2000. 

2.25 

4638. 

121.7 

72.0 

200= 

20000. 

.  000 

0.0 

500. 

.38 

3934. 

95.3 

68.  6 

210= 

W 

£j700. 

.  000 

28.5 

25000. 

.33 

46. 

1 . 6 

220= 


230=  CONSTANT  PARAMETERS 


240= 

250= 

ION  THRUSTER  MASS  <K3) 

S 

51.36 

260= 

THRUSTER  FUEL  TANK  MASS  (KS> 

s 

12.00 

270= 

SPECIFIC  MASS  (KS/KW) 

= 

10. 00 

230= 

THRUSTER  POWER  (KW/ THRUSTER;- 

3.06 

290= 

ION  THRUSTER  THRUST  (MN) 

129.00 

300= 

EXPENDED  MASS  AT  ETS  (KG) 

38399.00 

310= 

SHUTTLE  REFERENCE  MASS  <KS) 

= 

85000. 00 

320= 

BALLISTIC  COEFFICIENT  <M*#2/KG) 

.02 

330= 

SCALE  HEIGHT  (KM) 

= 

30.00 

.  340= 

CHEMICAL  OTV  MASS  (KG) 

s 

2270.00 

1^50= 

^360= 

SOC  MASS  (KG) 

s: 

100000.00 

ANNUAL  SOC  CARGO  MASS  (KG/YR> 

= 

20000.00 

370= 

CHEMICAL  OTV  SPECIFIC  IMPULSE 

sz 

455.00 

380= 

ION  OTV  SPECIFIC  IMPULSE  (SEC) 

S 

2900.00 

390= 

SOC  SPECIFIC  IMPULSE  (SEC) 

= 

2900.00 

4(.>0= 

CPB  SPECIFIC  IMPULSE  (SEC) 

= 

313. 00 

4 1  (.■>= 
420= 
430= 

-UEL  CONSUMPTION 

ORB I  TER  FUEL  MASS  (KG/ MSN) 

s 

4390.33 

440= 

SOC  FUEL  MASS  <K3/YR) 

= 

293.96 

450= 

OTV  FUEL  MASS  (KG/YR) 

= 

14122.37 

460= 

TOTAL  STS  FUEL  MASS  (KG/YR) 

s 

22328.21 

470= 

SOC  LOCATION  AND  RESUPPLY  PARAMETERS 

480= 

490= 

ORBIT  INCLINATION  (DEG) 

28.50 

500= 

ORBIT  RADIUS  (KM) 

6756.71 

510= 

NUMBER  OF  OTV  ION  THRUSTERS  (#) 

= 

54.  19 

520= 

NO  OF  ORB  LAUNCHES  (LS/YR) 

1.72 

530= 

AVG  OTV  TIME  OF  FLIGHT  (DAYS) 

120.00 

540= 

ION  OTV  MASS  (KG) 

5392. 16 

550= 

NO  OF  ION  PROPELLED  OTV’S  (#) 

= 

.99 

560= 

ORB  PAYLOAD  MASS  (KG /MSN) 

s 

28378.75 

570= 

PAYLOAD  MASS  TO  SOC  (KG) 

= 

27194.41 
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Computer  Run  D-8  Table  4.15  Efficient  Operatine  Point 

for  ETR  Scenario  B,  Missions  4-7, 
Level  I  Technology,  Run  Number  112 


I 

I 


100=  ETR  MODEL  VI  DUAL  OBJECTIVE  MINIMIZATION  PROORAM 


120= 

130=  SATELLITE  TRAFFIC  MODEL 

140= - 


150=MSN 

160= 

ASAT 

ESAT 

I  SAT 

MSAT 

FSAT 

OTV 

FUEL 

DEPLOY 

TIME 

RETURN 

TIME 

170= 

180= 

# 

KM 

DEG 

KG 

#/YR 

KG 

DAYS 

DAYS 

190= 

4 

65000. 

.000 

55.0 

1500.  3.00 

2842. 

83.7 

53. 8 

200= 

5 

25000. 

.700 

65.0 

1500. 

.  38 

2357. 

68.7 

45.3 

210= 

6 

12000. 

.000 

90.0 

4500.  2.25 

5058. 

170.2 

74.5 

220= 

7 

6700. 

.000 

98.0 

8000.  3. 75 

7938. 

290.6 

93.3 

230= 

240=  CONSTANT  PARAMETERS 
250= - 


260= 
270= 
230= 
290= 
300= 
310= 
320= 
330= 
340= 
350= 
360= 
370= 
380= 
390= 
400= 
410= 
A2C  = 


ION  THRUSTER  MASS  (KG) 

THRUSTER  FUEL  TANK  MASS  (KG)  = 
SPECIFIC  MASS  <KG/KW) 

THRUSTER  PQUjER  <KW/ THRUSTER) 

ION  THRUSTER  THRUST  (MN)  = 

EXPENDED  MASS  AT  ETS  (KG) 

SHUTTLE  REFERENCE  MASS  (KG) 
BALLISTIC  COEFFICIENT  (M»*2/K6)= 
SCALE  HEIGHT  (KM)  = 

CHEMICAL  OTV  MASS  (KG)  = 

SOC  MASS  (KG)  = 

ANNUAL  SOC  CARGO  MASS  (KS/YR>  = 
CHEMICAL  OTV  SPECIFIC  IMPULSE  = 
ION  OTV  SPECIFIC  IMPULSE  (SEC)  = 
see  SPECIFIC  IMPULSE  (SEC) 

ORB  SPECIFIC  IMPULSE  (SEC)  = 

FUEL  CONSUMPTION 


51.36 
12.00 
10.  00 
3.  06 
129. CO 
38399.00 
85000. 00 
.02 
30.  00 
2270.00 
100000.00 
20000.00 
455.00 
2900.00 
2900.00 
313.00 


430= - 

440=  ORB I TER  FUEL  MASS  (KG/MSN) 

450=  SOC  FUEL  MASS  <KG/YR) 

460=  OTV  FUEL  MASS  (KQ/YR) 

470=  TOTAL  STS  FUEL  MASS  (KQ.'YR)  = 

480=  SOC  LOCATION  AND  RESUPPLY  PARAMETERS 

490= - 

500=  ORBIT  INCLINATION  (DEG) 

510=  ORBIT  RADIUS  (KM) 

520=  NUMBER  OF  OTV  ION  THRUSTERS  (#)= 

530=  NO  OF  ORB  LAUNCHES  (LS/YR)  = 

540=  AVG  OTV  TIME  OF  FLIGHT  (DAYS)  = 

550=  ION  OTV  MASS  (KG)  = 

560=  NO  OF  ION  PROPELLED  OTV’S  (#)  = 

570=  ORB  PAYLOAD  MASS  (KG/MSN)  » 

580=  PAYLOAD  MASS  TO  SOC  (KG)  = 


4762.92 
995.30 
54462. 15 
82054.50 


57.00 
6721.92 
46.30 
5.58 
220. 00 
4697.76 
2.41 
21731.42 
21094.33 
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Computer  Run  D-9 


Table  4.15  Efficient  Operating  Point  for 
WTR  Scenario  B,  Missions  4-7,  Level  I 
Technology,  Run  Number  113 


■00= 

110= 

120= 

130= 


WTR  MODEL  VI  DUAL  OBJECTIVE  MINIMIZATION  PROGRAM 


SATELLITE  TRAFFIC  MODEL 


140= 

— 

- - 

— 

150=MSN 

160= 

ASAT 

ESAT 

I  SAT 

MSAT 

FSAT 

OTV 

FUEL 

DEPLOY 

TIME 

RETURN 

TIME 

170= 

180= 

« 

KM 

DEG 

KG 

#/YR 

KB 

DAYS 

DAYS 

190= 

4 

65000. 

.000 

55.0 

1500. 

3.00 

2190. 

99.4 

58.2 

200= 

5 

25000. 

.700 

65.0 

1500. 

.38 

1592. 

71.3 

43.3 

210= 

6 

12000. 

.000 

90.0 

4500. 

2.25 

3650. 

191.4 

71.3 

220= 

230= 

7 

6700. 

.  000 

98.0 

8000. 

3.75 

5907. 

337.  1 

88.0 

240= 

250= 

260= 

270= 

280= 

290= 

300= 

310= 

320= 

330= 

340- 

350= 

CoO* 

370= 

380= 

390= 

400= 

410= 

42'.)= 

430= 

440= 

450= 

460= 

470- 

480= 

490= 

500= 

510= 

520= 

530= 

540= 

550= 

560= 

570- 

580= 


CONSTANT  PARAMETERS 

ION  THRUSTER  MASS  <KG) 

THRUSTER  FUEL  TANK  MASS  (KG) 
SPECIFIC  MASS  (KG/KW) 

THRUSTER  POWER  <KW/ THRUSTER) 

ION  THRUSTER  THRUST  (MN) 
EXPENDED  MASS  AT  ETS  (KG) 
SHUTTLE  REFERENCE  MASS  (KB) 
BALLISTIC  COEFFICIENT  (M**2/KB) 
SCALE  HEIGHT  (KM) 

CHEMICAL  OTV  MASS  (KB) 

SOC  MASS  (KB) 

ANNUAL  SOC  CARSO  MASS  (KS/YR) 
CHEMICAL  OTV  SPECIFIC  IMPULSE 
ION  OTV  SPECIFIC  IMPULSE  (SEC) 
SOC  SPECIFIC  IMPULSE  (SEC) 

CRB  SPECIFIC  IMPULSE  (SEC) 

FUEL  CONSUMPTION 


ORBITER  FUEL  MASS  (KG/MSN) 

SOC  FUEL  MASS  (K6/YR) 

OTV  FUEL  MASS  (KG/YR)  = 

TOTAL  STS  FUEL  MASS  (KG/YR) 

SOC  LOCATION  AND  RESUPPLY  PARAMETERS 

ORBIT  INCLINATION  fDEG) 

ORBIT  RADIUS  (KM)  = 

NUMBER  OF  OTV  ION  THRUSTERS  (#)= 
NO  OF  ORB  LAUNCHES  (LS/YR) 

AVQ  OTV  TIME  OF  fiklQHT  (DAYS)  = 
ION  OTV  MASS  (KG) 2  = 

NO  OF  ION  PROPELl  WI^TV'S  (#)  = 

ORB  PAYLOAD  MASS  '  tKB^N) 

PAYLOAD  MASS  TO  SOC  (KG)  = 


51.36 

12.00 

10.00 

3.06 

129.00 

38399.00 

85000.00 

.02 

30.00 

2270.00 

100000.00 

20000.00 

455.00 

2900.00 

2900.00 

313.00 


8239.82 

5648.36 

40213.27 

91660.95 


60.06 
6670.23 
31.43 
5.56 
240.00 
3255. 46 
2.63 
20075.22 
19638. 13 
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Computer  Run  D-10  Table  4.15  Efficient  Operating  Point 

for  ETR  Scenario  A,  Missions  1-7,  Level  II 
Technology,  Run  Number  110 


ETR  MODEL  VI  DUAL  OBJECTIVE  MI^4IMI ZATION  PROORAM 


)=  SATELLITE  TRAFFIC  MODEL 


140= - 

150-MSN 

ASAT 

ESAT 

I  SAT 

MSAT 

FSAT 

OTV 

DEPLOY 

RETURN 

160= 

170= 

# 

KM 

DEG 

KG 

#/YR 

FUEL 

KG 

TIME 

DAYS 

TIME 

DAYS 

180= 

190= 

1 

41000. 

.000 

0.  0 

2000. 

2.25 

3653. 

61.9 

35.7 

200= 

2 

20000, 

.  000 

0.0 

500. 

.38 

a 

51.0 

36.2 

210= 

o 

6700. 

.  000 

28.5 

25000. 

'*>000 

a 

67.8 

10.3 

220= 

4 

65000. 

.000 

55.0 

1500. 

O  a  00 

2629. 

42.8 

27.5 

230= 

5 

25000. 

.700 

65.0 

1500. 

.38 

2756. 

44.9 

28.7 

240= 

6 

120C0. 

.000 

90.0 

4500. 

2.25 

5934. 

111.6 

47.0 

250= 

7 

6700. 

.000 

98.0 

8000. 

3.75 

10272. 

210.4 

64.  1 

260= 
270= 
280= 
290= 
300= 
310= 
320= 
330= 
.  .340= 
^S0= 
^^60= 
370= 
380= 
390= 
400= 
41.0  = 
420= 
430= 
440= 
450= 
460= 
470= 
480= 
490= 
500= 
510= 
520= 
530= 
540- 
550= 
560= 
370= 
580= 
590= 
--*-600= 
r.lvjlO- 


=  CONSTANT  PARAMETERS 


ION  THRUSTER  MASS  (K6) 

THRUSTER  FUEL  TANK  MASS  (KS) 
SPECIFIC  MASS  (KG/KW) 

THRUSTER  POWER  (KW/THRUSTER) 

ION  THRUSTER  THRUST  (MN) 

EXPENDED  MASS  AT  ETS  (KB) 

SHUTTLE  REFERENCE  MASS  (KG) 
BALLISTIC  COEFFICIENT  (M**2/Ke>= 
SCALE  HEIGHT  (KM) 

CHEMICAL  OTV  MASS  (KB) 

SOC  MASS  (KG) 

ANNUAL  see  CARGO  MASS  (KG/YR)  = 
CHEMICAL  OTV  SPECIFIC  IMPULSE  = 
ICN  OTV  SPECIFIC  IMPULSE  (SEC)  = 
SOC  SPECIFIC  IMPULSE  (SEC) 

ORB  SPECIFIC  IMPULSE  (SEC) 

FUEL  CONSUMPTION 


ORB I TER  FUEL  MASS  (KB/MSN) 

SOC  FUEL  MASS  (KB/YR) 

OTV  FUEL  MASS  (KG/YR) 

TOTAL  STS  FUEL  MASS  (KB/YR) 

SOC  LOCATION  AND  RESUPPLY  PARAMETERS 

ORBIT  INCLINATION  (DEG) 

ORBIT  RADIUS  (KM) 

NUMBER  OF  OTV  ION  THRUSTERS  (#)= 
NO  OF  CRB  LAUNCHES  (LS/YR) 

A VS  OTV  TIME  OF  FLIGHT  (DAYS)  - 
ION  OTV  MASS  (KS) 

NO  OF  ION  PROPELLED  OTV’S  (#)  = 

ORB  PAYLOAD  MASS  (KS/MSN) 

PAYLOAD  MASS  TO  SOC  (KG) 


51.36 
12.00 
5.  00 
6.86 
298.00 
38399.00 
85000.00 
.02 
30.  00 
2270.00 
1 00000 . 00 
20000 . 00 
455. 00 
.”448.  00 
2900. 00 
1  .  00 


4850.57 
.00 
85966. 65 
116554.06 


39.02 
7087.77 
43.59 
6.  ..>1 
120.00 
4557.63 
2.30 
26304.03 
24166.57 
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L 


Computer  Run  D-11 


Table  4.15  Efficient  Operating  Point  for 
ETR  Scenario  A,  Missions  1-7,  Level  III 
Technology,  Run  Number  111 


P 


1 00= 
110= 
120= 

ETR  MODEL  VI 

DUAL  OBJECTIVE  MINIMI ZATI 

ON  PROGRAM 

130= 

SATELLITE 

TRAFFIC  MODEL 

140= 

— 

- - 

- 

150=MSN 

ASAT 

ESAT 

I  SAT 

MSAT 

FSAT 

OTV 

DEPLOY 

RETURN 

160= 

FUEL 

TIME 

TIME 

170= 

180= 

# 

KM 

DEG 

KG 

#/YR 

KG 

DAYS 

DAYS 

190= 

1 

41000. 

.  000 

0.0 

2000. 

2.25 

2549. 

20.6 

12.9 

200= 

2 

20000. 

.000 

0.0 

500. 

m  >^8 

2315. 

17.2 

13.2 

210= 

3 

6700. 

.000 

28.5 

25000. 

.38 

2527, 

28.5 

4.7 

220= 

4 

65000. 

.  000 

55.0 

1500. 

3.00 

1768. 

13.  7 

9.5 

230= 

5 

25000. 

.700 

65.0 

1500. 

.38 

1323. 

14.2 

9.8 

240= 

6 

12000. 

.  000 

90.0 

4500. 

2.25 

3735. 

33.  9 

15.9 

250= 

260= 

7 

6700. 

.  000 

98.0 

8000. 

3.75 

6520. 

63.9 

21.9 

270= 

CONSTANT  PARAMETERS 

280= 

— 

- - 

— 

290= 

ION  THRUSTER 

MASS 

(KG) 

= 

40.00 

300= 

THRUSTER  FUEL  TANK 

MASS 

(KG) 

S 

12.00 

310= 

SPECIFIC  MASS  (KG/KW) 

= 

1.00 

320= 

THRUSTER  POWER  (KW/THRUSTER) 

10.00 

330= 

ION  THRUSTER 

THRUST  (MN> 

=  500.00 

.  .  340= 
1^30= 
360= 
370= 
380= 
390= 
400= 
410= 
420= 
430= 
440= 
450= 
460= 
470= 
480= 
490= 
500= 
510= 
520= 
530= 
540= 
550= 
560= 
570= 
580= 
590= 
-■">00= 
■  *10= 


EXPENDED  MASS  AT  ETS  (KS) 

SHUTTLE  REFERENCE  MASS  <KS) 
BALLISTIC  COEFFICIENT  <Mt*2/KB)= 
SCALE  HEISHT  <KM) 

CHEMICAL  OTV  MASS  (KG) 

SOC  MASS  (KG)  = 

ANNUAL  SOC  CARGO  MASS  <KG/YR>  = 
CHEMICAL  OTV  SPECIFIC  IMPULSE  = 
ION  OTV  SPECIFIC  IMPULSE  (SEC)  = 
SOC  SPECIFIC  IMPULSE  (SEC) 

0R3  SPECIFIC  IMPULSE  (SEC) 

FUEL  CONSUMPTION 


ORB I TER  FUEL  MASS  (KG /MSN) 

SOC  FUEL  MASS  (KG/YR) 

OTV  FUEL  MASS  (KG/YR) 

TOTAL  STS  FUEL  MASS  (KG/YR)  = 
SOC  LOCATION  AND  RESUPPLY  PARAMETERS 

ORBIT  INCLINATION  (DEG) 

ORBIT  RADIUS  (KM) 

NUMBER  OF  OTV  ION  THRUSTERS  (#)= 
NO  OF  ORB  LAUNCHES  (LS/YR) 

AV0  OTV  TIME  OF  FLIGHT  (DAYS)  « 
ION  OTV  MASS  (KG)  = 

NO  OF  ION  PROPELLED  OTV'S  <#)  = 

ORB  PAYLOAD  MASS  (KB/MSN) 

PAYLOAD  MASS  TO  SOC  (KG) 


38399.00 
80000.00 
.02 
30.00 
2270.00 
200000.00 
26250.  t:,-0 
453 .  0<;' 
5'.)0«  > .  00 
2900. 00 
313-00 


4685.21 

.00 

65068.40 

87849.65 


41.00 
7356.41 
76 . 56 
4.86 
40.00 
5046.57 
.77 
31015.72 
27469. 14 
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